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ABSTRACT 

Using  Van  Hove’s  density  correlation-function 
technique  we  have  calculated  the  intensity  spectrum 
of  light  scattered  by  a  relaxing  liquid.  Our  treatment 
differs  from  the  previous  work  on  the  subject,  due  to 
Mountain,  in  two  respects:  (a)  the  liquid  is  assumed 
to  have  a  relaxing  instead  of  a  static  shear  viscosity, 
and  (b)  instead  of  using  the  time  integral  representation 
for  the  relaxing  bulk  viscosity,  we  consider  it  to  arise 
from  the  relaxation  of  an  ordering  parameter,  and  apply 
the  thermodynamic  theory  of  relaxation  processes  to 
calculate  the  correlation  function.  The  first  modifi¬ 
cation  is  necessary  if  the  theory  is  to  explain  the 
occurrence  of  Brillouin  components  in  highly  viscous 
liquids  like  glycerine.  The  second  modification  gives 
an  insight  into  the  nature  of  approximations  involved 
in  using  the  integral  representation  for  the  bulk 
viscosity . 

We  first  solve  the  scattering  problem  for  a  model 
in  which  the  shear  and  the  bulk  relaxation  processes 
are  each  characterized  by  a  single  relaxation  time.  Our 
analysis  shows  that,  provided  the  sound  absorption  per 
cycle  is  small  compared  with  unity,  the  spectral  intensity 
may  be  approximately  written  as  a  sum  of  five  Lorentzians. 
Two  of  these  represent  the  Brillouin  components,  and  the 
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other  three  form  the  central  component.  It  is  also  shown 
that  the  time  integral  representation  of  the  bulk  vis¬ 
cosity  cannot ' in  general  describe  a  temperature  induced 
relaxation  process,  but  is  appropriate  for  a  density 
induced  relaxation  process. 

We  next  solve  the  scattering  problem  for  the  case 
when  the  relaxing  viscosities  are  each  described  by  a 
multiple  or  a  distribution  of  relaxation  times,  instead 
of  a  single  relaxation  time.  Both  exact' and  approximate 
expressions  for  the  intensity  spectrum  are  obtained  for 
this  case. 

Further  it  is  pointed  out  that  one  should  distin¬ 
guish  between  the  sound  velocity  v  for  a  temporally 
absorbed  hypersonic  wave  and  the  velocity  vg  as  determined 
from  the  Doppler-shifted  Brillouin  peaks.  The  dispersion 
relations  for  v  and  Vg  for  the  case  of  a  single  relaxation 
time  are  discussed  in  detail,  and  the  conditions  for 
negative  dispersion  (decrease  in  v  or  Vg  with  increasing 
wave  number)  are  obtained.  The  difference  between  v  and 
vD  is  significant  whenever  the  sound  absorption  per  cycle 

D 

is  not  much  smaller  than  unity. 
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GENERAL  INTRODUCTION 

t 

The  subject  of  light  scattering  is  nearly  a  hundred 
years  old,  dating  back  to  the  work  of  Rayleigh  [l]  and 
Einstein  [2j .  It  has  been  a  useful  tool  for  the  study  of 
a  variety  of  molecular  phenomena.  In  the  last  decade, 
with  the  development  of  laser  techniques  and  high-precision 
detection  devices,  light  scattering' has  opened  up  new 
possibilities  for  the  study  of  molecular  kinetics  and  for 
the  investigation  of  acoustic  properties  of  fluids  in  the 
ultra  high  frequency  (hypersonic)  region  (frequency 
~  101(^  c/s),  which  is  otherwise  difficult  to  study  by 
the  usual  ultrasonic  methods. 

It  is  well  known  that  the  scattering  of  light  by 
a  homogeneous  medium  mainly  arises "from  the  local  fluctua¬ 
tions  in  the  density  p  of  the  medium  due  to  thermal 
agitations.  These  density  fluctuations  give  rise  to 
fluctuations  in  the  dielectric  constant  e  of  the  medium 
and  hence  to  scattering  of  light.  A  Fourier  analysis  of 
the  density  fluctuations  over  the  volume  of  the  scattering 
medium  results  in  a  very  large  number  of  hypersonic  waves 
each  specified  by  a  wave  vector  it.  These  sound  waves 
act  substantially  as  diffraction  gratings  for  the  scatter¬ 
ing  of  light.  The  spectrum  then  consists  of  two  Doppler 

shifted  or  Brillouin  components  corresponding  to  sound 
propagation  in  opposite  directions  (for  a  given  it).  In 
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addition  there  is  a  central  or  Rayleigh  component  which 
has  a  peak  at  the  incident  frequency.  This  component 
arises  from  a  nonpropagating  (diffusive)  mode  which  decays 
as  entropy  fluctuations. 

Both  the  Brillouin  components  and  the  Rayleigh 
component ■ have  certain  line  widths  which  depend  on  the 
rates  at  which  the  propagating  (sound)  and  the  diffusive 
modes  decay.  It  is  only  in  the  last  few  years,  with  the 
advent  of  laser  techniques  and  the  availability  of 
intense  monochromatic  light  beams,  that  it  is  becoming 
increasingly  possible  to  measure  the  widths  and  line 
shapes  of  the  Rayleigh  and  Brillouin  components.  This 
in  turn  has  revived  interest  in  theoretical  calculations 
of  the  spectral  distribution  of  the  scattered  light. 

In  Chapter  1  of  this  thesis  we  briefly  review  some 
of  the  important  early  work  on  the  subject.  Some  aspects 
of  sound  absorption  in  fluids  and  of  the  thermodynamic 
theory  of  relaxation  processes  are  reviewed  respectively 
in  Chapters  2  and  3.  Chapter  4  describes  the  Van  Hove’s 
density-correlation  function  technique,  as  adapted  by 
Kamarov  and  Fisher  [20]  to  calculate  the  spectral  distri¬ 
bution  of  scattered  light.  These  chapters  form  the 
prerequisite  to  our  work  which  is  described  in  subsequent 
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In  Chapter  5  we  derive  the  spectral  distribu¬ 
tion  for  a  fluid  in  which  sound  is  absorbed  owing  to 
its  thermal  conductivity  and  to  relaxing  shear  and 
bulk  viscosities,  each  having  a  single  relaxation  time. 
The  relationship  of  our  work  to  those  of  previous 
authors  is  described  in  the  Introduction  to  this 
Chapter.  Chapter  6  points  out  the  necessity  of  dis¬ 
tinguishing  between  the  sound  velocity  as  measured  by 
Brillouin  scattering  and  that  by  other  (ultrasonic) 
techniques.  Chapter  7  is  devoted  to  the  extension  of 
the  theory  of  Chapter  5  to  the  case  of  a  fluid  in  which 
the  two  viscosities  are  each  characterized  by  a  set 
(discrete  and  continuum)  of  relaxation  times.  A  summary 
and  discussion  of  the  main  results  is  given  in  Chapter  8. 
The  calculation  of  the  spectral  intensities  involves 
expressions  for  the  mean  square  fluctuations  in  the 
various  thermodynamic  quantities.  These  are  derived 
in  the  Appendix. 
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Chapter  1 


EARLY  WORK 


§1.1  Rayleigh  Scattering 

In  an  experiment  carried  out  by  Tyndall  in  1869 
£3] ,  in  which  he  investigated  the  scattering  of  light 
by  fine  particles,  he  observed  a  bluish  hue  in  the 
scattered  spectrum  and  that  the  scattered  light  was 
polarized  at  an  angle  of  90°  to  the  direction  of  the 
incident  light.  Theoretical  investigation  of  the 
problem  was  later  taken  by  Rayleigh  in  1871  [l] .  He 
applied  the  electromagnetic  theory  to  a  system  of 
optically  isotropic  particles  whose  size  is  small 
compared  with  the  optical  wavelength.  The  incident 
field  is  then  assumed  to  be  practically  uniform  through¬ 
out  each  of  the  particles  and  induces  oscillatory  electric 
dipoles  in  them.  These  serve  as  secondary  sources  of 
light  and  radiate  light  in  the  form  of  scattering.  For 
observation  at  a  large  distance  from  the  scattering 
medium,  it  is  enough  to  consider  only  dipole  scattering. 
The  total  scattered  intensity  should  depend  on  the  super¬ 
position  of  the  individual  wavelets  with  the  proper  phases 
from  each  of  the  molecular  dipoles.  But  Rayleigh  assumed 
that  these  phases  are  at  random  so  that  the  total  inten¬ 
sity  is  the  sum  of  the  individual  intensities.  This  is 
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equivalent  to  assuming  that  there  is  no  correlation 
between  the  molecules,  which  is  justified  for  very 
dilute  systems" like  an  ideal  gas.  The  expression  for 
the  scattered  intensity  per  unit  volume  per  unit 
incident  intensity  is 


I 


^  4  2  0 

l6?  g01  sin2 
R  A 


(1) 


where  a  is  the  average  polarizability  of  the  particles, 

A  the  wavelength  of  incident  light,  n  the  number  density 
of  the  particles,  R  is  the  point  of  observation  assumed 
to  be  far  away  from  the  scattering  medium,  and  <!>  the 
angle  between  R  and  the  direction  of  the  incident  electric 
vector . 


If  the  incident  light  is  not  polarized,  Eq„  (1) 
must  be  averaged  over  all  directions  of  the  electric 
vector  in  a  plane  perpendicular  to  the  direction  of 
propagation  of  the  incident  light »  One  then  has 


I 


o  4  „  2  0 

ott  na  / -.  , 

— p - r—  ( 1  +  cos  0  )  , 

R  A 


(2) 


where  0  is  the  scattering  angle. 


§1.2  The  work  of  Einstein  and  Smolukowski 

As  we  stated  before,  the  assumption  which  one  makes 
in  Rayleigh’s  theory  of  light  scattering  is  that  the 
particles  are  random  scatterers.  For  dense  fluid  media. 
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especially  liquids,  one  should  expect ' that'  the  molecules 
will  not  move  independently  of  each" other,  and  therefore 
the  above  assumption  will  not  be  justified.  Indeed, 
experiments  show  that  the  scattering  formula  given  by 
Eq.  (1)  or  (2)  does  not  predict  correct  results  for 
liquids , 

Let  us  now  consider  a  homogeneous  liquid 
as  the  scattering  medium.  A  simple " application  of  Huygen’s 
principle  shows  that  such  a  medium" can ' scatter  no  light 
because  of  the  destructive  interference" from  neighbouring 
small  volumes,  except,  of  course,  in  the- forward  direction. 
This  is  contrary  to  observations.  One  can  therefore 
conclude  from  this  that  there  must  exist  in  the  medium 
some  inhomogeneities  which  are  responsible  for  the 
observed  scattering  of  light,  although  the  liquid  is  on 
the  average  homogeneous.  This  problem  was  solved  by 
Einstein  and  Smolukowski  [2,4^  who  first  used  the  concept 
of  thermodynamic  fluctuations  to  explain  the  large 
intensity  of  scattering  observed  for  gases  near  their 
critical  temperatures,  a  phenomenon  known  as  critical 
opalescence.  According  to  them,  the  molecules  in  a 
liquid  are  not  distributed  uniformly,  but  rather  their 
number  fluctuates  from  their  mean  value  in  any  small 
volume  element  of  the  liquid.  These  fluctuations  arise 
from  thermal  agitations  of  the  molecules  and  they  in 
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turn  produce  fluctuations  in  the  dielectric  constant 
+ 

e  .  Thus  the  scattering  system  is  really  not  homo¬ 
geneous  locally.  So  unlike  the  single  particle 
scattering  theory  of  Rayleigh,  the  theory  of  Smolukowski 
and  Einstein  considers  scattering  by  small  volumes  in 
the  fluid,  so  that  though  there  may  be  molecular  corre¬ 
lations,  the  small  volumes  are  large  enough  to  be 
statistically  independent  of  each  other,  and  scatter 
light  incoherently. 

Let  us  divide  the  scattering  medium  into  small 
regions  or  cells  whose  linear  dimensions  are  small 
compared  with  the  optical  wavelength  X  but  are  large 
enough  to  contain  sufficiently  many  molecules  to  allow 
one  to  use  thermodynamic  description.  Consider  a  par¬ 
ticular  cell  of  volume  6V . .  Let  the  dielectric  constant 

J 

at  some  instant  of  time  be  e.  =  e  +  6e.,  where  e  is 

J  o  j  o 

the  mean  dielectric  constant  of  the  medium  and  5e>is 


the  fluctuation. 

The  relevant  electromagnetic  equations 
Gaussian  units)  for  our  system  are 


J 

(in 


The  fluctuations  in  the  dielectric  constant  can 
also  arise  from  the  fluctuations  in  the  orientation 
of  the  molecules,  if  they  are  anisotropic.  This  effect 


will  not  be  considered  here. 
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curl 


3  =  —  D  a  div  D  =  0 


curl  3  =  -  H  ,  dlv  3=0 

c  5 

6  =  e  £  =  i  +  4*  ?  ,  ?  = 


(3) 


Hence  a  fluctuation  6en.  will  result  in  an 

J 

additional  polarization  AP,  giving  rise  to  the  scattered 
-*■  ->  -*■ 

fields  E’j  D’  and  H’ .  We  are  only  interested  in  linear 
effects.  Assuming  that  |6e.  |<<E0,  |  E  *  |  << ( E |  etc.,  we 
have 


cur 


1  Hf  =  -  D’  ,  div  6'  =  0 


curl  E’  = 


14-  -* 

-  —  H  3  div  H’  =  0 
c  5 


3’  =  e  3’  +  6e3  =  e  3’  +  4tt AP 

O  0  0 


A?  =  -rf^-  E 


Ae 

TF  o 


(4) 


where  E  is  the  incident  electric  vector.  Eliminating  H’ 
o 

and  3’  from  (4),  one  gets 


v  3 ’ - 1  S’  = 

c" 


-  curl  curl  (6e.E  ) 

J  ° 


(5) 


We  now  define  t  =  t,  so  that  (5)  is  transformed 


into 


V2D*  -  —r  S’  =  -  curl  curl  (6e. 3  )  . 
2  ^  2  jo 

C  0  T 


(6) 


The  solution  of  this  equation  is  readily  obtained 


in  terms  of  retarded  potential 


3  6  +  ’ 
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D’  =  curl  curl  n 


s  = 


/ 

6  V 


A?  ( ryi ) 


d3? 


-*• 


“  i 

R-r 


T  - 


±  -> 
R-r 


57  ' 


(6e.lV 


<SV.  R-r 


Q-V 

d3r 

 ll=2 


(7) 


As  the  volume  6V .  has  dimensions  much  smaller  than  A, 

J 

(7)  can  be  written  as 


n 


1  (5£.iE0)6V1 

57  - 


-v 

R-r 


zr  —r 

R-r 


(8) 


T  - 


Here  R  is  the  point  of  observation  and  r  .  is  the  position 

of  the  volume  6V . , 

J 

We  shall  be  interested  only  in  observations  at 

large  distances  from  the  scattering  medium,  so  that 

S-r  .  I  =  R, 

J  1 

The  incident  electric  vector  E  can  be  written 

o 

as 


*  3  — i  (  k  . r.-ft  t) 

E0  =  E1  e  o  j  c 


(9) 


where  is  the  angular  frequency  of  the  incident  light 

and  k  is  its  wave  vector .  Since  the  wave  number  of  the 
o 

^  o/ 

scattered  wave  k’  is  equal  to  — - —  ,  (8)  can  be  written 
as 


r  ift  t 

'le  °  -Ik’cR^  rir  N  i ( k  1  —  k  ).r. 

n  =  57r  <6ej  6vj5  e  0  J 


(10) 
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Now  regarding  e  as  a  function  of  p  and  T,  one  has 


9  e  •  9  e  . 

6ej  =  (9  p  '  T  6p  +  ^9~T“^p  6T  ° 


(11) 


We  assume,  as  customary,  that  the  second  term  can  be 

•j* 

neglected  compared  with  the  first  .  Hence, 


6ej  ~  (3  p}T  6pj  In'T  6nj 


(11'  ) 


and  (10)  becomes 


j  _  g  „ift„t  _-i<’  oR  1  /9  e 


1 


e  o  e 


.  /  \  -*■ 

,-o  tx  1  (.  K  T  -K  )  .  r  . 

( ^  „  )  m  o  N  j  e  o  j 


^ttR  9  n '  T  ,i 


(12) 


where  6N .  (E  6n . 6V . )  Is  the  number  fluctuation  In  the  volume 
J  J  J 

6V.c  To  get  E’,  one  must  evaluate  curl  curl  ff .  For  large 
J 


R, 


curl  curl 


E  e1^  R 
E  e 

R 


-]  i  ic  t  p 

^  curl  curl  E-^1 


=  —  {grad ( dlv  S1e1K?oR)  -  V2 (S-Le1K '  0 R)  } 


1  ’  **  *  R 

^  {- |  E-jJ  k  ’  k  ’  cose})  e1K  °  +  «,fcE^e 


,2g  e1K’ oR}  e  5  a  gt 


(13) 


This  neglect  Implies  that  (^-&) 2  <5TP>  .  «  (|“^0 2  <Sp2>  9 

dip  dpi 

2  2 

where  <6T  >  and  <6p  >  are  the  mean  square  fluctuations 
of  temperature  and  density  respectively.  There  is  some 
evidence,  though  not  conclusive,  that  for  water  this 
assumption  may  not  be  a  very  good  approximation.  See  f5] • 


n  ' 


11 


where  <p  is  the  angle  between  k  f  and  From  (13),  it 

can  be  verified  easily  that  2' . k  t  =  0,  and  that 


sintj) 


(14) 


Now  far  away  from  the  scattering  medium,  we  have 
A?  =  0,  so  that  (4)  and  (14)  give 


gt  =  i_  g i  =  g  (|-|-)m  6N,  el(K'  Ko)orj 


I  7T  £ 


O 


d  n  T 


(15) 


Hence 


I®'  I  w-  2  2  2  K’  sin  ^  ^n'T  6Nj  0  ('l6) 

1 6  7T  Re  d 

o 

The  intensity  per  unit  incident  intensity  is  then  given 
by 


I . 
J 


Q 

_ o 

IV  2^2  4 

lOTT  R  C 


sin  <t> 


M)2 

Kd  n  V 


2 

6N  .  > 

J 


(17) 


We  can  now  make  the  assumption  that  for'  a  liquid 
not  close  to  its  critical  point,  fluctuations  in  different 
cells  are  statistically  independent,  since  then  molecular 
correlation  extends " only  to  molecular  distances.  There¬ 
fore,  the  total  intensity  will  be  the  sum  of  the  contri¬ 
butions  from  all  the  cells,  i.e.  1=21..  From  the  thermo- 

j  J 

dynamic  theory  of  fluctuations  [see  Appendix  [|  , 


«5N .  >  = 
J 


kDT  n  6V . 
Boo  ,i 


B 


T 


3 


(18) 
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where  nQ  and  Tq  are  the  equilibrium  number  density  and 
temperature  respectively.  Therefore  the'  total  intensity 
is 


I  = 


V  knT 
B  o 

16tt2BtR2 


^4  2  0 

tp(|f)Tlp=Po  sin  * 


(19) 


where  V  is  the  total  scattering ' volume .  The  expression 
(19)  is  the  well-known  result  of  Einstein. 

For  an  ideal  gas,  PV=N  ku  T  ,  so  that  B=n  k^T. 
Further  the  polarizability  a= ( e-1 )/4im ,  so  that  one  may 
write 


—  1  /8  £  \ 

°  =  ^  °  T  ‘ 

If  one  substitutes  (20)  into  (19),  one  sees' that 
becomes  identical  with  the  Rayleigh  formula,  eq. 
above . 


(20) 

(19) 

(1) 


Now  in  order  to  make  (19)  useful  for  practical 
purposes,  one  needs  to  know  how  e  and  p  depend  on  each 
other.  Einstein  used  the  Clausius-Mossott i  or  Lorentz- 
Lorenz  relation 


e-1  1 

—To  —  =  constant 

£  +  2  p 


Therefore 


/9e  \ 
po  3  p  T 


(e0-  1)Uo+  2) 

3 


(21) 


(22) 


If  (22)  is  substituted  into  (19) ,  the  resulting  formula 
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is  known  to  hold  well  for  dilute  gases  (e  -  1)  but  not 
as  well  in  dense  fluids  and  rather  poorly  for  polar 
molecules.  Ramanathan  and  others  [6,7^  have  ascribed 
this  to  the  inadequacy  of  using  the  Lorentz  local  field 
approximation.  According  to  these  authors,  a  more 
satisfactory  relation  to  use,  for  light  scattering 
problems,  is^* 

e  -  1 

p  =  constant  (23) 

If  we  substitute  (23)  into  (19),  we  have 
V  knT  ft  i  o  o 

I  =  - p - j  (— )  (e  -  1)  sin  (p  .  (24) 

16tt  BtR^  c  ° 

Expression  (24)  is  known  to  give  better  agreement  with 
experiments  for  most  liquids  than  the  use  of  (21)  in 
(19)  [6,7]. 

§1.3  Fine  Structure  of  the  Rayleigh  Line 

A .  Brillouin  Mandelstam  Doublets 

The  fine  structure  of  the  Rayleigh  line  was  first 
predicted  by  Brillouin  in  1922  [8J  ,  and  later  also  by 
Mandelstam  £9].  In  analogy  to  the  Debye  waves  in  solids. 


Mathematical  justification  for  this  step  is  lacking 


at  the  present  time. 


. 

B  1'  •* 
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Brlllouin  considered  that  the  description  of  the  thermal 

motions  of  the  molecules  in  a  liquid  is  equivalent  to  a 

superposition  of  a  large  number  of  elastic  sound  waves 

of  different  wave  vectors  k.  Provided  such  a  Fourier 

analysis  is  carried  out  over  a  volume  of  the  scattering 

medium  whose  linear  dimensions  are  all  large  compared 

with  the  wavelength  of  light,  each' Fourier  component 

substantially  acts  as  a  diffraction  grating;  In  other 

words  the  wave-vectors  k  and  k’  of  the  incident  and 

o 

scattered  light  beams  respectively  are  related' to  the 

wave-vector  k  of  the  sound  wave  by  the  well  known  Bragg 

condition  k ' -k  =  k.  One  can  arrive  at  this  formula 

o 

most  easily  if  in  (7),  one  considers  the  contribution 
to  n  from  the  whole  scattering  volume,  so  that 


4- 

n  = 


Then 


(15) 


(6eS) 

|  =4  -4  | 

R-r 


T 


becomes 


—4 

G 


(— ) 
v9  p'T 


r  x 
J  6pe 

V 


-4 


-k  ) . r  ,  3+ 
o  d  ur 


(25) 


(26) 


Fourier  analysing  6p,  noting  that  the  elastic 
sound  waves  are  standing  waves,  we  write 


6p 


v  -4  v  “4-4 

£  p ( k )  cos  k . r 

t 

i  v  / -i  \  r  ik.r  .  -ik . r-i 
^  £  P(k)  [e  +  e  J 

k 


(27) 


r ol flhnoo 
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Substituting  this  into  (26),  and  making  use  of 

the  fact  that  the  integral  /elc*'r  dpr  differs  from  zero 

V 

for  large  volumes  of  V  and  is  equal  to  V  only  under 
the  condition  q=0,  we  obtain 

=  Str  5  (ff}T  *  v&(k)  +  p(-k)]  ,  (28) 

o  p 

with  the  condition 

±k  =  k1  -  kq  .  (29) 


Expression  (28)  will  lead  to  Einstein 1 s ■ formula 
(19)  if  one  makes  use  of  the  following  results  from  the 
theory  of  fluctuations  [see  Appendix^ 


<  p(k) p(-k)>  = 

po  B  o 

(30) 

v  bt 

<  p ( k) p (k)  >  = 

0  . 

(31) 

If  0  denotes  the 

angle  of  scattering,  then  the 

Bragg  condition 

(29)  may 

be  written  in  the  more  familiar 

form  as 

2^ot 

k  =  -  sin  % 

c 

0 

(32) 

or 

2  A  sin  e=  A 

3 

(33) 

where  y  is  the  refractive  index  of  the  medium,  A  is  the 
wavelength  of  the  sound  wave  and  A  is  the  optical  wave¬ 
length  in  the  medium. 


* 

Jbni  V  la  esniii'.ov  io 


.  c  , 
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Now  if  the  homogeneities  were  static,  the  scattered 
light  would  have  the  same  frequency  as  the  incident  beam. 
Since  the  sound  waves  are  moving,  the  time  varying  inhomo¬ 
geneities  will  modulate  the  scattered  light.  Let  the 
latter  be  denoted  by 

S'  a  einot  <(.(t )  ,  (3*0 

where  cf>(t)  is  the  modulation  function.  We  shall  see  later 
that  <J>(t)  arises  from  fluctuations  in  pressure  or  density 
at  constant  entropy,  which  propagate  through  the  medium 
with  the  sound  velocity  v.  Hence  to  a  first  approximation 
<J>(t)  may  be  expected  to  obey  the  wave  equation 

<J>  ( t )  -  v2V2<f>  =  0  .  (  35  ) 

Since  the  standing  wave  solution  of  (35)  is  rele¬ 
vant  ,  we  have 

4>  ( t )  a  cos  ajt  cos  it.r  ,  (36) 

where  u)=vk.  Therefore,  we  have 

f  ift  t  r  iujt  ,  — io)t-r 

E  a  e  o  e  +  e 
p  L  J 

=  eiD^o+aj^  +  e1Pfio"a)'t^  .  (37) 

Thus  in  the  scattered  light  there  are  two  components. 
One  has  angular  frequency  ft  -w,  and  is  called  the  stokes 
line.  The  other  has  angular  frequency  ftQ+a),  and  is  called 
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the  anti-Stokes  line.  From  (32),  the  shift  in  frequency 
is  given  by 

Aft  =  w(ic)  =  2  — y  ft  sin  %  6  .  (38) 

This  result  is  consistent  with  the  fact  that  two 

opposite  going  sound  waves  give  rise  to  a  Doppler  shift 

effect  of  precisely  the  same" quantity  as  (38). 

The  Brillouin  component s *  obtained  from  (37)  are 

infinitely  sharp.  In  practice  they  have  finite  line 

widths.  We  shall  see  later  that  this  is  because  sound 

waves  are  absorbed  in  the  medium. 

From  (33)  it  can  be  seen  that' for  incident  light 

in  the  visible  spectrum,  A  is  of  the  order  of  5  *  10 

cm.  Since  the  velocity  of  sound  is  of  the  order 

4 

10  cm/sec,  the  frequency  of  sound  wave  involved  is 

q  io 

of  the  order  of  10^  or  10  cycles  per  second.  Sound 
waves  of  such  high  frequencies  are  called  hypersonic 
or  ultra  high  frequency  waves. 

In  liquids  since  shear  waves ' are ' highly  damped, 
only  longitudinal  hypersonic  waves  are  propagated. 
Therefore  only  a  single  pair  of  Brillouin  lines  are 
observed . 

B .  The  Central  or  Rayleigh  Component 

In  the  first  experimental  determination  of  the 
frequency  distribution  of  the  light  scattered  by  a  liquid. 


. 


i 
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Gross  QlcQ  observed  that  In  addition  to  the  doublet 
predicted  by  Brillouin’s  theory,  a  central  line  was 
also  present.  The  explanation  of  this  was  given  by 
Landau  and  Placzek  QllJ  who  considered  that  the  spon¬ 
taneous  density  fluctuations' in  a  liquid,  besides  having 
an  adiabatic  pressure  fluctuation  part  to  give  rise  to 
Brillouin’s  doublets,  also  has  a  part  due  to  isobaric 
entropy  fluctuation.  The  latter  will  also  scatter  light. 
Let  this  be  described  by  another  modulation  function  i|;(t). 
The  contribution  to  the  scattered  light  will  be 


E’  a  e^o1"  i|;(t)  . 
s 


(39) 


The  heat  transport  equation  is  normally  given  by 
[see  Sec  2.3  3 


(40) 


where  x  is  the  thermal  conductivity,  and  all  quantities 
are  referred  to  unit  mass.  For  isobaric  fluctuation, 
(40)  can  be  written  as 

If  =  V2  T  .  (41) 

3  t  pc 

o  p 

If  we  assume  that  ip ( t )  satisfies  a  similar  equa¬ 
tion,  then 


(42) 
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Eq.  (42)  can  be  solved  easily  and  (39)becomes 


E_  a  e  o  e 


P  c 

Ko  p 


t 


(43) 


This  means  that  the  entropy  wave  is  nonpropaga- 
tlve  and  damped  out  temporally.  This  mode  therefore 
represents  an  unshlfted  line,  and  is  conventionally 
called  the  Rayleigh  component. 


C .  Temporal  Absorption  and  Linewidths 

A  damping  factor  such  as  e-^  ^/P0cp  ^3)  is 

equivalent  to  temporal  absorption.  The  spectral  dis¬ 
tribution  of  the  intensity  can  be  obtained  by  a  Fourier 
transform  of  (43).  Writing 

oo 

e'  (t)  =  /  S'(io’)  e1"^  doo'  (44) 

o  S 

—  oo 

+ 

we  have 

oo  •  M  t 

e’  (a)")  =  -  f  £'(tf)  e"la)  dt'  .  (45) 

S  TTqS 

Putting  (43)  into  (45),  and  defining  the  frequency 
of  the  scattered  light  in  terms  of  the  frequency  shift, 

a)1  =  w"  -  ft  (46) 

o 


We  have  restricted  to  positive  t.  This  is  equivalent 

to  assuming  Ef(t)  =  £ ' ( | t | ) .  See  for  example  Ref.  29 

s  s 

or  (5  -  14). 


•• 


v  vr.  ;  Jf;  -- 


q°cc  s  fious  ‘toJ'obI:  ^niqmBb  A 
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we  get 


I  (wf )  a 

O 


aj 


-2  +  (X 


k‘ 


(47) 


p  c 
Ko  p 


This  is  a  Lorentzian  distribution  with  maximum 
at  the  unshifted  frequency  i.e.  mt=0  and  with  half  width 

xiiL 


equal  to 


P  c 
Ko  p 


Now  consider  the  Brillouin  components.  Temporal 
absorption  can  be  formally  described  by "letting  the 
angular  frequency  go  in  (37)  be  complex.  Denoting  the 
real  part  of  the  angular  frequency  to  be  go  again,  we 
can  write  (37)  as 


S' (t )  a  e 
P 


-Tt  r  i(fi  +aj)t  ,  i(ft  -a>)t 
|_e  o  +  e  o 


■J  . 


(48) 


where  ai=vk,  v  being  the  sound  velocity,  and  r  is  the 

temporal  absorption  and  arises  from  the • imaginary  part 

of  the  complex  frequency.  The  result  of  (48)  can  also 

be  obtained  if  one  adds  to  (35)  a  loss  term  -  TV  4> . 

If  we  now  do  similar  thing  as  before  and  evaluate 

S'(go’)  from  (48),  using  (46),  one  obtains,  for  small 
P 

sound  absorption  (i.e.  f/vk  <<  1), 


± » 


I '  -  const.  {— tt 
P  -2 


1 


f  +(mf--vk) 


r2  +  (a)f  +vk) 2 


} 


(49) 


Thus  we  obtain  a  doublets  which  are  symmetrically 
centred  at  w’  =  ±vk  with  half  width  f  for  each  of  the  two 


component  s . 
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D 


Ratio  of  the  Intensities  of  the  Fine  Structure 


The  ratio  of  the  integrated  intensity  of  the 

central  Rayleigh  component,  denoted  by  I  ,  to  the 

c 

intensities  of  the  displaced  Brillouin  doublets, 
denoted  by  2  Ig,  was  first  given  by  Landau  and  Placzek 
En].  Their  result  is 


I  C  —  C-TJ- 

c  =  P  V 
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B 


'V 


=  y  -  1  , 


(50) 


where  c  and  c,T  are  the  specific  heats  at  constant 
p  V 

pressure  and  volume  respectively,  and  y  is  the  specific 

heats  ratio.  (50)  is  then  known  as  the  Landau-Placzek 

ratio.  This  can  be  derived  as- follows: 

The  total  intensity  of  the  scattered  light  is 

2 

proportional  to  <  6e  >  ,  where  6e  is  the  dielectric 

fluctuation.  If  one  neglects  (^-ft)  »  then  one  can  write 

dip 

<  6e2>  =  (|^-)2  «5p2>. 

Now  6p=  (3_|)p  6S  +  0S  6p  . 

From  thermodynamic  calculations  [see  Chapter  3l  , 


P)  =  _o  /9  Pn 
8  p  s  Bs  ;  v9  s'p 


P0eTo 

— 5 


Q2 

3  =  p. 


c  -c,r 

_P _ Z 

TBm 


(51) 


From  fluctuation  theory  Qsee  Appendix^ , 


c  6s2>  =  kBc  ,  <  6p2>  =pokB  T^B, 


<  6s  6p>  =  0 


(52) 


s  * 


q  ;  ■: 
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Making  use  of  (51)  and  (52),  one  obtains 


*  2 

6  p  >  = 


p°knT 
Ko  B  o 

Bm 


{(1  _  i)  +  i 
Y  Y 


Hence 
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B 


=  (1  -  7V7  =  Y  -  1 


(53) 


(50) 


A  good  account  of  what  we  have  just  described  in 
this  chapter  can  be  found  in  "Molecular  Scattering  of 
Light"  by  Pabelinskii  [37]  • 


Chapter  2 


SOUftfD  ABSORPTIONS  IN  VISCOUS  LIQUIDS 
§2.1  Preliminary  Remarks 

We  have  seen  in  the  last  section  that  the  absorp¬ 
tions  of  sound  waves  in  a  liquid  modulate  the  scattered 
light.  Landau  and  Placzek  [llj  first  related  the  width 
of  the  central  line  of  the  Rayleigh  fire  structure  to 
the  heat  conductivity  and  the  width  of ' the* doublets  to 
the  absorption  of  sound  in  the  liquid.  Later  works  [l2j  have 
confirmed  this  theory.  In  this  section  we  shall  discuss 
some  of  the  major  causes  of  sound  absorption,  namely  the 
shear  and  bulk  viscosities  and  the  thermal  conductivity 
\j-3]  •  We  have  said  that  a  calculation  of  density  fluctua¬ 
tions  by  thermodynamic  theory  alone  cannot  give  the  line 
widths.  In  the  most  general  case,  the  kinetic  motions  of 
the  molecules  in  the  liquid  should  be  understood  from  a 
formalism  based  on  a  microscopic  theory.  But  this  is 
extremely  difficult  to  do.  However,  for  the  problem  of 
light  scattering  or  ultrasonic  absorptions  in  liquids, 
since  the  wavelengths  of  disturbances  are  much  larger 
than  the  int ermolecular  distances,  one  can  use  macroscopic 
treatment  of  the  problem,  and  the  dissipative  mechanisms 
can  be  conveniently  discribed  by  macroscopic  equations  of 
hydrodynamics  and  thermodynamics  of  irreversible  processes. 


<  £  ^  i 1  c  «,;■  >r* -j  b.i  jci  rl  ;n 
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§2,2  Hydrodynamical  Equations 

In  a  fluid,  conservation  of  matter  implies  that 
if  p  is  the  density  of  the  fluid  and  u,  the  velocity  at 
a  given  point  in  the  fluid,  we  have 

/(pu)  .d§  =  -  gY  /  pdV  .  (1) 

This  holds  for  any  arbitrary  volume  .  Hence  we 

have 


|f  +  div  (pu)  -  0  ,  (2) 

which  is  the  continuity  equation.  For  small  disturbances, 
(2)  can  be  linearized,  which  then  becomes 

|f  +  PQ  div  u  =  0  ,  (3) 

For  a  force-free  fluid,  the 'force  required  to  move 
an  element  of  fluid  is  only  to  overcome  the  hydrostatic 
pressure,  which  acts  normal  to  the  surface  of  the  element. 
Thus 

/  p|£  dV  =  -/pd§  =  -/VpdV  .  (4) 


Hence 


oat 


=  -  Vp  , 


(5) 


which  is  the  linearized  Euler’s  equation. 


. 


39tnoo9d  t  r  o  riw  tt>ds^B9ni  !  6  o  (!  ) 
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In  a  real  fluid,  however,  the  stress  on  a  surface 
element  contains  both  the  normal  and  the  tangential  com¬ 
ponents,  so  that  one  has  to  describe  it  by  a  symmetric 
tensor  a..  .  of  rank  2,  Hence  the  force  density  equation 
of  Euler  becomes 


3  u 


3  a 


o  3  t 


3  x  . 
J 


JA  = 


-  E 


n.  3  7T  .  . 

+  Z 

J  3  xj  13  j 3  xj 


(6) 


where 


°ij  _  p6ij  +  "ij  ’ 


(7) 


and  7T .  .  ( — it  .  .  )  is  called  the  viscosity  stress  tensor, 
ij  Ji  J 

We  assume  that  the  rates  of  strain  3  u^/3  x  .  are  small 
so  that  the  stress  can  be  taken  to  depend  linearly  on 
them,  and  that  the  fluid  is  isotropic,  so  that  the 
properties  along  any  axis  are  the  same.  Therefore, 


3  u 


7T  .  . 
li 


=  r) '  div  u  +  2 n 


3  x . 
i 


7T 


IJ 


3  u .  3  u  . 

n  L— -  +  .-e1) 


3  x  . 
J 


3  x. 


i  /  J 


where  n  and  r\ '  are  constants 


Putting  (8)  into  (6),  we  obtain 


3_u 
*o  3  t 


=  -  Vp  +  b  grad  (div  u) 


(8) 


+  ^  n  grad  (div  u)  +  nV  u  , 


(9) 


where 


? 


b  =  n'  +  -r  n  . 


(10) 


9rT9flW 


vi  ;  11.30  it  ,  6 

'  •  0  1  »(J  £12?  i  8  )  -  «  j  'v-ri  O  3 


. 


27 


Stokes  assumed  the  special  case  b  =  0,  so  that  (9) 
is  reduced  to  the  Navier-Stokes  Equation, 

n  is  called  the  shear  viscosity  and  b  the  bulk 
viscosity . 


§2.3  Bulk  Viscosity  and  Sound  Absorption 


It  can  be  shown  that  for  ideal  monatomic  gas  b=0 
[14]  .  In  general  b^O  [13,15]. 

For  a  nonzero  b,  the  sound  absorption  a  (per  wave¬ 
length)  due  to  viscosities  alone  can  be  calculated  as 
follows : 

We  assume  that  sound  waves  are  of  weak  intensity 
so  that  linearized  hydrodynamical  equations  may  be  used. 
This  is  often  referred  to  as  thq  acoustic  approximation. 
Combining  (9)  with  (3),  We  obtain  for  one  dimensional 
flow, 

9  2u  6p  9  2u  ,  ,  4  N  9  3  2u  , 

— o  TT  — ~  +  Cb  +  5-r) )  jwr  — 2  *  (ID 

p  L  3x' 

solved  for  a  solution  of 


9  t 


u  op  9  x 


Eq.  (11)  can  be 

i(o)t-kx)  .  . 

u  a  e  ,  giving 


a)2  (4^-)-1  {1  +  iu)(b 
op 

2  [o_  1 

“  Bs  1  +  Wu>v 


4  v  1  /9 

+  3n  3  P 


J 


where 


03  .  = 


B 

s 

b  +  jn 


(12) 

(13) 


(14) 
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For  k=k^-ik,-)  ,  u 


-k0x  i(o)t-knx) 
a  e  2  e  1 

the  sound  attenuation  per  wavelength  is  defined  as 

27Tk0 


a 


k. 


(15) 


Substituting  (15)  into  (13),  we  have 


2 

7TV  p 


U)/a) 


a 


2 

v  = 


V 


B 


2B 


1+  (  a)/ 03 v  ) 


1+  (co/coy) 

l+{l+(u/(0v)2}^ 


For  the  low  frequency  case. 


(16) 


(17) 


CO  <  <  0) 


V 


(18) 


Then 


a  = 


Trk 


p  v 
Ko  o 


(b  +  5-  n) 


v  =  v  =  (B  /p  ) 2 
o  s  Mo; 


(19) 

(20) 


For  b  =  0,  Eq .  (19)  together  with  the  contribution 
to  sound  absorption  from  thermal  conductivity  is  fre¬ 
quently  called  the  classical  absorption.  It  has  been 
established  QL3,  16^  that  the  classical  absorption  is 
too  small  to  account  for  the  observed  absorption  in  many 
liquids  and  that  a  nonzero  b  is  necessary  to  explain  the 
additional  absorption.  Eq.  (19)  is  about  the  only  means 


of  determining  b  experimentally 
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To  understand  the  meaning  of  b,  let  us  calculate 
the  mean  pressure  from  (7),  (8)  and  (3). 


1  v  „ 

3  ?  °ii  =  P 
1 


+  L.  Le. 


P0  3t 


(21) 


We  see  that  b  is  a  measure  of  the  resistance  to  the 
rate  of  change  of  density  of  the  fluid,  which  may  be 
caused  by  different  mechanisms.  In  a  monatomic  fluid 
a  nonzero  b  can  arise  only  from  the  local  spatial 
rearrangement  of  the  atoms  i.e.  structural  relaxation 
[see  Chapter  3] ,  and  b=0  for  a  very  dilute  system  of 
such  constituents.  For  more  complicated  molecules,  b 
can  arise  also  from  other  types  of  relaxation  processes. 

In  some  liquids  [13,  16] ,  the  observed  sound 
attenuation  a  shows  a  maximum,  when  plotted  against 


co/ojv,  at  some  finite  value  of  u)/ujv 


Moreover  the 


calculated  values  of  both  v  and  a  for  w  >  wv  are  much 
larger  than  the  observed  values.  All  these  results 
contradict  the  predictions  of  Eqs.  (16)  and  (17)  in 
which  b  and  n  are  constants' independent" of  frequency. 
One  is  led  to  conclude  that  for  these  liquids  some 
relaxation  processes  must  be  taking  place  and  that  b 
and  n  only  describe  the  low  frequency  limits  of  these 
processes.  We  recall  that  b  and  n  were  introduced  in 
(8)  under  the  assumption  that  only  the  first  space 
derivatives  of  the  particle  velocity  u  are  kept.  For 
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higher  frequencies,  this  will  be  a  poor  approximation 
and  the  hydrodynamic al  equation  (9)  will  break  down. 

Hence  a  proper  way  of  modifying  the  stress  tensor  seems 
to  be  needed.  Since  this  is  not  easy  as  molecular 
motions  in  liquids  are  in  general  very  complicated,  one 
finds  it  more  convenient  to  assume  some  relaxing  forms 
of  the  viscosities  in  a  modified  hydrodynamical  equation, 
or  in  the  case  of  bulk  viscosity,  one  can  introduce 
thermodynamic  relaxation  theory,  as  will  be  discussed 
in  Chapter  3*  In  some  liquids  like  CS^  [17],  it  is 
sufficient  to  consider  b  to  arise  from  a  single  relaxa¬ 
tion  time;  in  others  [13,  1 61  ,  more  than  one  relaxation 

time,  or  even  a  distribution  of  relaxation  times,  is 
needed . 

In  viscous  liquids,  both  shear  and  bulk  viscosities 
are  relaxing  as  the  condition  oj  >  is  easily  satisfied. 

§2.4  Sound  Absorption  due  to  Thermal  Conductivity 

We  assume,  for  simplicity,  that  the  only  loss 
mechanism  in  a  fluid  is  due  to  thermal  conductivity;  that 
is,  there  are  no  viscosities.  The  heat  transfer  equation 
is 


*1 Q  _  m  d s  _  y 
dt  1  dt  P 

o 


(22) 


where  x  is  the  thermal  conductivity  and  s  is  the  entropy 
per  unit  mass. 
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Using  thermodynamic  relations  (see  Chapter  3) 
one  can  write 


m  9 _ _  m  /9  S  \  9  P  .  m  (9 _ S\  9  p 

o  3  t  o^3  p  p  3  t  ol3pyp  3t 


=  LL-  3_E  _  _J2  Le. 


6Bt  3  t 


P06  3  t 


(23) 


V2T  =  (| |)p  V2p  +  (||)p  V2p 


SB 


1  V2p  -  A  V2P 


po® 


(24) 


For  6p,  6p  etc.  a  e^^  (24)  gives 


i£  _  +  luicp) 

6p  =  P°  (Xk2  .  .  X 

(p  +  lu)Cv) 


(25) 


Putting  this  into  (11),  and  neglecting  viscosity. 


we  proceed  as  we  did  to  obtain 

2 

i  2  _  Po^2  xk  +iajpocV 

k  —  ^  *  p  . 

T  xk  +ia)P0Cp 


(26) 


Again  writing  k  =  -  ik^,  the  low  frequency  limit 

2 


0)  << 


c  p  v 
pMo 


(27) 


yields  for  a  and  v  the  expressions 


a  =  TT(y-l)xk 


p  V  c 
Ko  o  p 


(28) 


2  2  T3  / 

v  =  v  =  B  /p 
o  s  Ko 


(29) 
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If  viscosities  are  not  neglected,  it  can  be  shown 
easily  that  provided  (27)  and  (18)  hold,  the  sound  absorp¬ 
tions  due  to  heat  conductivity  and  viscosities  are 
additive.  That  is 


a  = 


7T  k 


p  V 
kO  O 


D> 


(30) 


33 


Chapter  3 


THERMODYNAMIC  THEORY  OF  RELAXATION 

We  have  mentioned  in  Chapter  2  that  in  certain 
liquids,  sound  absorptions ' can  be  considered  to  arise 
from  relaxation  processes  with  either  a  single  relaxa¬ 
tion  time  or  more.  It  was  also  pointed  out  that  the 
usual  hydrodynamic  equation  with  nonrelaxing  viscosities 
is  then  not  applicable.  Since  a  derivation  of  relaxing 
viscosities  from  the  detailed  molecular  motions  does 
not  seem  possible,  attempts  have  been  made- to " describe 
instead  the  general  relaxation  features  of  such  processes. 
We  shall  describe  in  this  chapter ' a  very  general  formalism 
called  the  thermodynamic  theory  of  relaxation,  which  was 
first  introduced  by  Mandelstam ■ and  Leontovich  [Yl 8J  ,  and 
then  developed  by  others,  in  particular,  very  extensively 
done  by  Meixner  [l9]  .  In  the  following,  we  shall  follow 
the  description  in  Ref.  13. 

§3.1  Formulation  of  the  Theory 

Consider  a  liquid  which  contains  a  single  ordering 
mechanism.  Such  a  system  needs  for  its  specification 
three  independent  thermodynamic  variables.  Two  of  these 
may  be  chosen  from  S,  T,p  and  P,  and  the  third  is  the 
"order”  parameter  (which  we  denote  here  by  £)  or  the 


,.:cJ  vo  .ioeJ  fin  r  '  ms  ;a  I  yd  beoubcridj  i  ts'iJtl 

i  t3n  wc  lo  driJ  nl  .  '.*>JJ 

.  iL  .  tsfl  n.  noi 


V  oim  n^Lonr  rt:  ,1  n  bn-  jbi  :.  am.hl 


34 


variable  Z  conjugate  to  For  thermodynamic  reasoning 

it  is  not  necessary  to  specify  For  example,  it  may 

refer  to  the  temperature  of  the  internal  degrees  of 
freedom  of  the  molecules  -  the  relaxation  of  £  then 
gives  rise  to  the  so  called  thermal  relaxation  process; 
for  a  chemical  reaction,  £  would  represent  the  degree 
of  reaction. 

In  all  the  thermodynamic  descriptions  throughout 
this  chapter  we  shall  refer* extensive  quantities  to  unit 
mass.  The  first  law  of  thermodynamics  may  be  written  as 

i 

6u  =  T6S  +  ^  6p  +  z6£  (1) 

po 

The  Gibb’s  free  energy  is 

6G  =  -  S6T  +  —  6p  +  Z6£  (2) 

po 

When  the  system  is  at  equilibrium  at  constant  T  and  p, 

G  is  a  minimum  so  that  we  have 

Z(T,  P,  €)  =  0  ,  (3) 

which  determines  the  equilibrium  value  of  E,  as  a  function 
of  T  and  P.  For  equilibrium  at  constant  T  and  p,  the 
Helmholtz  free  energy  is  a  minimum: 

SA  =  -SSI  +  ~  Sp  +  Z6S  =  0  ,  (4) 

po 

•f* 

The  theory  can  be  easily  extended  to  the  case  of  several 


order  parameters. 


.  .  _  .  . 


£.•  ^  1  9  -3 
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so  that 

Z(T,p,€)  =  o  .  (5) 

A  thermodynamic  quantity  measured  at" constant 
Z(=0)  corresponds  to  the  value  measured  in  the- equili¬ 
brium  state,  i.e.  at  zero  frequency.  On  the  other  hand, 
a  measurement  of  quantities  at  constant  value  of  the 
relaxation  parameter  £  corresponds  to  the  infinite 
frequency  case,  when  things  happen  at  too  fast  a  rate 
for  £  to  attain  its  equilibrium  value  £,  which  is 
determined  by  (3)  or  (5).  Such  quantities  are  sometimes 
called  'frozen'  quantities  and  are  denoted  usually  with 
a  superscript  A  quantity  without  a  superscript  will, 

for  simplicity,  correspond  to  one  measured  at  zero 
frequency . 

Since  when  £  has  its  equilibrium  value  £,  one 
may  to  a  first  approximation  assume  that  for  small 
deviation  6£=£-£,  the  system  returns  to  equilibrium 
according  to  the  equation 


d6£, 

dt 


L6Z  , 


(6) 


where  SZ  may  be  regarded  as  the  ordering  force,  and  L 
is  a  constant  or  a  slowly  varying  function  of  the  thermo¬ 
dynamic  variables  p,  T,  etc. 

If  we  write 


6Z  =  (!4) 


9  T'  p  ,  £ 


6T  +  (--) 


•9  p'T,£ 


X  _L  /3Zx 

6  p  +  U-cO 


•9  t  T ,  p 


(7) 
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and  make  use  of  (5),  Eq.  (6)  becomes 


6M 

dt 


XTV 


) 


where 


-1 

tTV 


L  ( 


a_z. 

8  S;T,p 


For  to  relax  at  constant  T  and  P,  we  have 


(8) 


(9) 


6M  =  _  &Z-&Z 

dt  TTp 

where  Is  now  to  be  determined  by  (3),  and 

T-1  =  T 

T  ~  yd  £'T,P 

Similarly  one  can  define  tqv  and  tqp  . 


(10) 


(ID 


§3.2  Some  Useful  Formulae 


We  shall  in  this  section  give  some  thermodynamic 
relations  which  will  be  useful  to  us  later.  First  we 
recall  some  of  the  general  thermodynamic  quantities  and 
relations  for  the  zero  frequency  limit. 


/9 

pol9  pjS,Z 

;  bt  = 

po(; 

=  T  (3-^) 

0^3  T'p,Z 

)  ^  p  — 

to(; 

=  _  i_(3_B.) 

Pq^9  T;P,Z 

■po^T)T 

,z  ; 

BS  _  cp 

Bs  ”  bt 

bt  cv 

p 

9  P 


9  S 


'V 


6  =  - 


o 


o  V 


5  =  (9-^) 

T  l9TJp,Z 


2  2 

3  Bt  ; 


2  9_S 

po^9  p  '  T ,  Z 


37 


3  (cP  cy)  bt 


(12) 


where  and  BT  are  the  adiabatic  and  isothermal  bulk 
modulus  respectively;  Cp  and  Cy  are  the  isobaric  and 
isochoric  specific  heats  respectively;  y  is  the  specific 
heats  ratio  and  3  is  the  coefficient  of  volume  expansion 
For  infinite  frequency  case,  we  immediately  have 
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(13) 


From  (4),  one  gets 


(3_S) 

'3  p;T,Z 


-  (9-^) 

2  l3Tjp,Z 

Po 


3B 


T 


=  (ML) 

Z;p,T  l3T;p,Z 


1_  (L2)  =  -(9-^-) 

2  v3  Z;T,p  v3  p'T,Z 

p  3 


(14) 


Similar  relations  can  be  obtained  from  the  other  Legendre 
transforms  of  (1).  Consider  now 


. 


38 


6P  (3T)p,£<5T  +  6p+  (3^T,p 


R 

=  6  VT  +  ^  6p  +  ff  T)  P  65 


(15) 


Therefore , 


Br 


OO 
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Similarly 
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07) 


(15)  also  gives 
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Consider  now 
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Therefore , 
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-  c°°  =  T  (— )  (9-I) 
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Also,  from  (19),  we  have 
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Similarly 
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It  can  be  easily  verified  that  we  also  have 
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_  i_  =  p~3 (3_£} 
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Combining  (20) ,(21)  and  (25),  we  have 
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Combining  (16),  (22)  and  (18),  we  have 
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Also,  from  (12),  (13)  and  (28),  we  have 
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We  shall  now  obtain  relations  between  the  various 


relaxation  times. 


From  (9)  and  (11),  we  have 
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Therefore 
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Combining  (31),  (32)  and  (33),  we  have 
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Adding  (39)  to  (36)  and  using  (26)  and  (32),  we  have 


^  +  TSP  dt^6p  Bs  ^  +  TSV  dt^6P  * 

(40) 

For  6p,  6P  a  eiajt,  we  have  from  (40), 


o  ,4  n  6P  „  ,  (BS“BS)la)TSV 

BS(lu)  =  6p  =  BS  +  1+ iwT  gY 

(41) 

Thus  we  have  arrived  at  a  complex  frequency-dependent 
bulk  modulus. 

B .  Isothermal  Processes 

If  we  consider  6T=0  and  do  similar  things  as  in 
A,  we  would  obtain,  similar  to  (40)  and  (4l),  the  follow- 


ing  results 

(1+tTP  dt^6p  =  B^(1  +  tTV  lt^6P 

(42) 

E  fin)  -  E  1  (BT-BT)iuTTV 

BT(im)  Bt  +  1+1wTtv 

(43) 

C .  Temperature  Induced  or  Thermal  Relaxation  Processes 

This  is  a  very  important  kind  of  relaxation  process, 
and  is  one  induced  by  a  change  of  temperature  only.  Thus 
for  such  a  process,  isothermal  changes  of  pressure  and 
volume  occur  without  relaxation  effects.  The  condition 


then  is  given  by 
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Adding  (39)  to  (36)  and  using  (26)  and  (32),  we  have 
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For  6p,  5P  a  e  ,  we  have  from  (40), 
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Thus  we  have  arrived  at  a  complex  frequency-dependent 
bulk  modulus. 

B.  Isothermal  Processes 


If  we  consider  6T=0  and  do  similar  things  as  in 
A,  we  would  obtain,  similar  to  (40)  and  (4l),  the  follow¬ 


ing  results 


( 1+T 


TP  dt 


)6p  B^(l  +  ttv  dt)6P 


00 


BT(ica)  =  Bt  + 


(  ®rp  -^rp  )  T rp-y- 


1+iooT 


(42) 


(43) 


TV 


C .  Temperature  Induced  or  Thermal  Relaxation  Processes 


This  is  a  very  important  kind  of  relaxation  process, 
and  is  one  induced  by  a  change  of  temperature  only.  Thus 
for  such  a  process,  isothermal  changes  of  pressure  and 
volume  occur  without  relaxation  effects.  The  condition 


then  is  given  by 


(<LP)  =  (<LZ)  =  0 

'9  '  iji  p  '  ip  p  vy  • 


(44) 


44 


In  (44),  since  (f"|)Tjp  =  -  (|-|-)T  p  >  one  of  the 

Bm 

two  conditions  in  (44)  will  suffice.  From  (25)  ,  this 
condition  is  also  equivalent  to 

Bt  =  B~  ,  (45) 

or  from  (43),  B^ico)  =  Bp  is  real  and  independent  of  go. 
Therefore  there  is  no  relaxation  of  the  isothermal  bulk 
modulus . 

If  a  process  is  thermally  induced,  then  from  (28), 

we  have  8  =  3°°,  which  together  with  (45),  (12)  and  (13) 

00  00 

imply  that  cp  -  Cy  =  cp  -  Cy.  These  are  the  consequences, 
not  the  condition  for  thermal  relaxation. 


D .  Pressure  Induced  Relaxation  Processes 

These  are  processes  which  are  induced  by  a  change 
of  pressure  only,  irrespective  of  whether  the~relaxat ion 
takes  place  adiabatically  or  isothermally .  Such  processes 
are  designated  by 


(?_§.)  =  (?J£)  =  o 

O  ^  X  P  O^^SP  W  * 

From  (21),  this  condition  is  also  equivalent  to 


(46) 


00 


If  we  write 


6p 


<&>P,e4s 


+ 


(?_£)  <sp 
^  P;S,£ 


(9_£) 

^9  £  ;S,P 


(4?) 


«1 


45 


Then 


(<Le.)  =  (!_£)  C— )  +  (^-2-) 

T,P  V9  S;P,£V9  C  T3P  v9£;S,P 


or 


(4_e.)  =  (i^.) 

o  P  o  £  ^  S  P  * 


(48) 


From  (25),  (26)  and  (33),  we  can  write 

2 


00 


KT  KT  ^9  ^T,P  Cp 


00  = 


KS  KS  (g-f)s5p  Cp 


(49) 


Now 


(^_P)  =  (3_P)  +  (l_P)  ftJL) 

o  S  P  a  ^  T  P  9  T y P  E,  d  E,  S  P 


=  (l4)m  o  +  (M) 


k9  C  T,P 


oo  "9  E,  yT,P 


(50) 


Using  (20)  and  (21),  we  have 


00 


V  (3_R) 

cp  v3  C  S,P 


=  fLPi  +  PoeTo  (ZA-) 
d  E,J T  P  c  v9  ^  ;T  P 


(51) 


Hence  (49)  becomes 


00 


^p  Kp 


(Le.)2  c“ 

^3  g;T,P  P 


oo 


Ks  Ks 


C  f(8-P)  +  MT(3_S)  1 2 

cpL^a  £;T,P  Cp  v9  £  '  T ,  P- 


(52) 


Eq.  (52)  provides  us  with  a  criteria  for  temperature  or 
pressure  induced  relaxation  process.  Thus  if 


•  r  i  i  ns  3  a  t  £Z')  ..  ff-  (o:>  t  ( <  !l ,  rio’i  l 
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(— ) 

Kd  e  t5p 


<< 


P  $T 
o  o 


3_Ss 

3  T,P 


or 


oo 


Kp  Kg  Kg  , 


we  have  thermal  relaxation,  and 


00 


Kp  Kpp 


or 


oo 


Bp  -  Bt  . 


(53) 


(54) 


(55) 


On  the  other  hand,  if 


(L£) 

£;T,P 


>> 


po3To  ,3  S v 
cD  ^jT,P 


we  have 


(56) 


K 


T 


00 


(57) 


This  corresponds  to  pressure  induced  relaxation. 

In  analogy  with  the  above  discussion,  we  may 
define  also  a  density  or  volume  induced  relaxation 
process  by  requiring  that 


(9_S) 

9  6  T,p 


(— ) 

S,p 


(58) 


which  differs  from  (46)  in  that  the  role  p  and  P  have 
been  exchanged.  For  this  case,  we  have 


c 


V 


00 


C 


V 


} 


(59) 


) 
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and 


(9_Z)  =  (9_Z) 

a  £  y  T ,  p  v9  £  ;  S  ,  p 


(60) 


Similar  to  (53)  and  (56),  we  have  now  the  condition 


(82) 

WT,p 


<< 


gBTTo  f3_S, 
cy  °  £  x ,  p 


(61) 


or 


BT  -  bt  <<  Bs  “  Bs 


(62) 


for  temperature  induced  processes,  and 


(8-2) 
k3  C;T,p 


>> 


gBTTo  (9_S) 


'V 


9  5'T,p 


(63) 


or 


bt  _  bt  ~  Bs  ”  Bs 


(64) 


for  volume  or  density  induced  processes.  We  may  note 
that  eqs.  (62)  and  (54)  are  essentially  equivalent, 
as  they  should  be.  However,  eqs.  (64)  and  (57)  are, 
in  general,  not  equivalent. 


§3.4  Representation  of  Bulk  and  Shear  Viscosities 

For  the  case  of  a  fluid  system  in  which  the  only 
cause  of  sound  absorption  is  a  relaxation  process  with 
a  single  relaxation  time,  the  sound  absorption  and  sound 


. 
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velocity  can  be  calculated  from  Eq.  ( 4 1 )  together  with 
the  hydrodynamical  equations  (2-3)  and  (2-5).  For 
6P  and  6p  proportional  to  e1 ^ ‘ r ^ ,  (2-3)  and  (2-5) 
give 


6P 

6p 


03 


(65) 


Eq.  (65)  is  obtained  from  Euler’s  equation, 
which  has  been  derived  for  force-free  fluid.  We  wish 
to  make  a  remark  here:  If  -r—  in  (65)  is  real  and  is 

6p 

equal  to  Bg/p,  then  eq.  (65)  cannot  describe  sound 
absorption.  To  obtain  sound  absorption,  the  Euler’s 
equation  has  to  be  replaced  by  Navier-St okes  equation 
or  other  modified'  force  density  equation.  An  alternative 
way  is  what  we  are  doing.  The  Euler’s  equation  is  used 
but  in  (65)  is  complex,  so  that  sound  absorption  is 
automatically  taken  into  account  in  the  thermodynamics 
in  terms  of  an  ordering  parameter  £ ,  which  is  what  we 
have  discussed  in  this  chapter.  The  first  way  of  doing 
things  requires  a  hydrodynamic  equation  involving  the 
bulk  and  shear  viscosities,  which  is  not  applicable 
at  high  frequencies.  It  is  interesting  to  see  how  the 
bulk  viscosity  arises  in  the  thermodynamic  theory  of 
relaxation.  Substituting  (4l)  into  (65)  ,  one  gets 

1  +  lOJT 


k2  = 


sv  w  Bs 


bs+bs1(a)Tsv 


(66) 


v 


o 
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Now  the  description  of  sound  absorption  in  terms 
of  a  bulk  viscosity  b  gives,  if  one  puts  n=0  in  (2-13), 


2 

a)  p 
+  iwb 


(67) 


Substituting 

h  <Bs-Bs>* 

b  =  - 

1  +  icot  gy 


(67) 

SV 


into 


(66), 


we  have 


(68) 


Actually  from  the  manner  b  is  introduced  in  (2-10), 
it  cannot  be  frequency  dependent ," which  is  precisely  why 
the  hydrodynamical  equation  of- Navier-Stokes  breaks  down 
at  high  frequencies.  Hence  strictly  speaking,  one 
cannot  substitute  (67)  into  (66)  to  obtain  (68).  The 
proper  way  should  be  to  compare  the  low  frequency  limit 
of  (67)  with  that  of  (66),  because  it  is  in  this  limit 
that  description  in  terms  of  the  hydrodynamical  equation 
is  meaningful.  Doing  this,  one  obtains  b  as  the  low 
frequency  limit  of  the  relaxation  process 


b  (Bg  Bg)  tsv  .  i&9) 

However,  (68)  provides  us  a  way  to  describe  a  relaxing, 
or  frequency-dependent  bulk  viscosity. 

In  a  similar  way,  one  can  describe  the  shear 
viscosity  n  as  the  low  frequency  limit  of  a  shear 
relaxation  process 


} 


n 


(70) 


,  .  V  ,t  L  SO  0  1  7  >i  ■  >!•  '  4  0 


c  )  -  -  t  -■  1 J  ( '  :  'J 


limJ  'i 


js^oo'iq  nolJ“BXBl9'i 
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00 

where  t  is  the  shear  relaxation  time  and  GQ,  GQ  are 

T")  o  o 

the  shear  modulus  at  zero  and  infinite  frequencies 

00  oo 

respectively.  In  liquids,  Gg=  0  and  Gs  =  G^ .  Hence 

(70)  becomes 


n 


OO 

G  x 


n 


(71) 


A  relaxing  shear  viscosity  then  can  be  written  as 
00 

G  T 

/ .  x  n 

nCirn.)  1+1(j0X 


(72) 


n 

As  it  is  well  known  that  the  frequency  dependence 
of  a  quantity  can  be  considered  as  resulting  from  the 
Fourier  transform  of  a  time  dependent  quantity,  one  can 
assume  that  there  exists  a  function  b'(t)  for  the  bulk 
viscosity  of  the  following  form 


b' (t)  =  (B~  -  Bs)  e_t/TSV 


(73) 


Then 
b  ( ico ) 


/  e  la)t '  b ' (t '  )dt ' 


0 


(Bs-Bs 


1  +  iWT 


)tsv 

SV 


as  in  (68) 


Similarly  the  time  dependent  function  for  shear  viscosity 
is 

n'(t)  =  G°Vt/Tn  (74) 

Thus  in  describing  a  thermodynamic  relaxation 
process,  one  can  either  employ  an  ordering  parameter 
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£,  which  adds  a  new  degree  of  freedom  to  the  thermo¬ 
dynamics  of  the  system,  or  one  can  use  (68).  If  one 
prefers  the  latter  way,  one  has  to  find  a  suitable 
way  to  introduce  the  relaxing  bulk  viscosity  into  the 
hydrodynamical  equation.  We  recall  that  for  the 
unrelaxed  case,  the  term  in  the  hydrodynamical  equa¬ 
tion  involving  b  is,  from  (2-9),  b  grad(div  u)„ 

Therefore  an  integral  representation  of  a  time 
dependent  function  of  the  following  form  is  suitable 
t 

c(t)  =  /  b’(t-t’)  grad  (div  u(t'))  dt ?  ,  (75) 

0 

where  b'(t)  is  given  by  (73)*  For,  the  frequency 
dependence  of  c(t)  can  be  obtained  by  a  Laplace 
transform 

00  00  00 

/  c(t)dt  =  /  dt  /  e”^^  bf(t-t’)gnad  div  u(t')dt' 

0  0  0 

00  .  H  00  .  J 

=  /  dt  "b ’  ( t"  )e_ia)t  /  e-^^  grad  div  u(t’)  dt ?  .  (76) 

0  0 

The  first  factor  in  (76)  gives  the  correct  frequency 
dependence  as  in  (68). 

We  shall  in  a  later  chapter  discuss  the  difference 
between  the  two  ways  of  representing  the  relaxing  bulk 
viscosity . 


* 

.  f  3  .  9  1  3  >  n  il  v.  )0"  ..  \J£: 


.  ( l  V )  vd  n  \:  .c  {  ’  d  i  »ri\. 

3  .  <J  D9nlB^cfc  d  itfij  ( d }  o  lo  &9n  thnoqsb 
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Since  the  shearing  processes  do  not  involve 
changes  in  volume  or  temperature  of  the  fluid,  a 
relaxing  shear  visocity  is  most  conveniently  described 
by  an  integral  representation  similar  to  the  one  out¬ 
lined  above  for  the  bulk  viscosity. 
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Chapter  4 

MICROSCOPIC  FORMULATION  OF  THE  SCATTERING  PROBLEM 

The  early  formulations  of  the  light  scattering 
problem,  described  in  Chapter  1,  do  not  allow  one  to 
calculate  readily  the  spectral  distribution  of  the 
scattered  light.  Recently  Kamarov  and  Fisher  [2  0]  and 
Pecora  [2l]  have  adapted  Van  Hove’s  correlation  function 
technique  (developed  in  connection  with  neutron  scatter¬ 
ing)  to  this  problem.  In  this  chapter  we  first  briefly 
describe  Van  Hove’s  theory  and  then  show  how  it  is  trans¬ 
lated  from  the  language  of  neutron  scattering  to  that  of 
light  scattering. 

§4.1  Van  Hove’s  Theory  of  Neutron  Scattering  [22] 

Consider  a  target  system  S  containing  N  atoms,  the 
position  of  the  ith  atom  being  denoted  by  r^ .  Let  the 
incident  neutron  of  mass  m  have  momentum  fik^  .  It  is 
scattered  by  S  into  a  state  with  momentum  ,  while  S 
changes  from  the  initial  state  | I>  with  energy  E^  to  the 
final  state  |F>  with  energy  E„o  The  transition  probability 

r 

per  unit  time  is 

W ( k , a) )  =  ^|  <F,k2|Hint  |l,k1>|2p(eF)6(fi(A)-Ep+EI)  ,  (1) 


where  Hint  (r)  is  the  interaction  potential;  ep=^p+ 


2m 


■ 
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2  2 
■f)  lc 

is  the  energy  of  the  final  state;  =  fiu  =  Ep  -  Ep 

is  the  change  in  energy  of  the  neutron;  6 (ftw-Ep+E^) 
expresses  conservation  of  energy  and  p(ep)  is  the  density 
per  unit  energy  range  of  the  final  state. 


Passing  into  continuum  for  the  space,  we  have 
m  V 


p(0  = 


F'  (2,)V 


(2) 


The  differential  scattering  cross-section  per  unit  solid 
angle  per  unit  energy  range  is 


d 2o  _  ^2  ,  m  x  2  2 

d^de  -  kx  27Tft2 


<F ,k2 | Hint (r ) | I,k1>| 26 (hm-Ep+EI ) 


(3) 


Now  in  general  there  is  a  statistical  distribution  of 
initial  target  states  | I> ,  each  with  a  probability  P(I), 
and  will  end  up  in  a  certain  |F>.  Therefore,  we  have 


dfc_ 

dftde 


V2(~-)(-~)2i:  P(I)  |  <F,k2|Hint(?)  |l,k  >|  2S(-nu-E  +EJ.) 
1  2-nfi  I,F 

(4) 


Usually , 

Hint  (r) 

where  a. 

i 

particle . 
<k2  |  Hint  ( 


Hint  (r)  can  be  written  in  the  form 


2  71"fi  „  j.  /  "►  \ 

-  E  a .  6 ( r-r . )  , 

m  .i  i  5 

i 


(5) 


is  an  operator  depending  on  the  spin  of  ith 
Therefore , 

-*%!,  1  ,  ,  3  +  v  27Th  r  ■>  \  i(k,-kn).r 

r)lk1>  =  v  f  d r  *  ai  ~m  6(r-ri)e  1  2 


2Trfi  „  iic.r. 

-iT—  E  a.  e  i 
mV  .  l 

l 


(6) 
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We  now  replace  the  6-function  in  (4)  by  the  integral 


VEI 

1  j.  e  i(w  ^  )t  ^  *  and  substitute  (6)  into 


2  7T"h 


(4)  to  obtain 


d2c  = 
dfide  2ufi  'k 


1  {tt~)  f  dt  e"lajt  £  £  P(I) 


I5P  i,J 


x  <  F  e 


iH  t  iH  t 

*  a±  elL?i(o)  e’”3-  1 1> 


x  <  I | a .  e-lk'rj(o)|F>  , 
J  J  '  5 


(7) 


where  we  have  used  the  equation' Hq | I>  =  E^|I>,  Hq  being 
the  unperturbed  Hamiltonian,  and  we  have  also  chosen  our 


initial  time  as  the  origin.  We  now  denote 

E  P(’I)  <T  |  ....  1 1>  by  the  ensemble  average  <. 
I 

E | F  >  <F  |  =  1,  (7)  can  be  written  as 
F 


>.  Since 


dja  =  1  (52)  s  fe~i^  <  a  e-iK-?i(o)a.  (t)eiL?j  (t)>dt  , 

dftde  2'7iTi  kn  .  .  1  J 

-L  1  3  J 

(8) 

where  r.(t)  and  a.(t)  are  now  Heisenberg  operators,  i.e., 

J  HH5  1Hofc 

^  ( t )  =  e  *  r  j  ( o )  e  ^  , 


with  similar  expression  for  a.(t) 

J 
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Consider  now  the  expression 


/d3x  d3J'  eik-(x-x,)<  6(x’-r1(o))a(x-rj(t))> 


-+  -* 


=  </  d^x'  e  ik,x  6  (x  '  -r^(o)  )fd^x  6(x-r.(t))e 


ik .  x 


> 


-ik . r . ( o )  iK.r.(t) 

=  <  e  1  e  j  > 


(9) 


Therefore,  (8)  becomes 


d!o_ 

dfide  2ttTi  vk 


1  d3^'  el5-(x-x,)/dt  e-lut 


x  <  E  a .  (o )  6  ( x  '  -r .  ( o  ))&  ( t )  6  ( x-r .  ( t ) )  > 

A  A  -L  “*■  J  J 

±  )  J 


(10) 


If  we  make  the  transformation  x-x'=r,  x=r’,  and  use  the 
even  property  of  the  6-function,  we  obtain 


d2a 


1  /k2w,3+  ,3-*-,  i(K.r-ajt) 

(t— ) /d  r  dJrf  dt  e 


d£2de  2t\Ti  k 


1 


x  <  E  a . ( o ) 6 (r+r . ( o ) -r ’ )a  .  ( t ) 6 (r f -r  .  ( t ) )  > 

i,j  1  1  J  J 


(11) 


Define  now  the  pair  correlation  function 


f(r,t)  = 


N”1 <  E  /d^rf  a . ( o ) 6 (r+r . ( o ) -r ’ ) a . ( t ) 6 (r ’ -r . ( t ) ) > 

A  A  1  J  J 

y  J 


(12) 


Therefore  (11)  becomes 


dfo_ 

dftde  2 Trfi  xk 


N  (_2)/ei(k.r-ut)  r(j  t)  a3?  dt 


(13) 
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One  can  define  another  pair  correlation  function 
G(r,t) 

G(r,t)  =  N_1  <  Z  /d^r ' 6 (r+r . ( o ) -r ’ ) 6 (r ’ -r . ( t ) ) > .  (14) 

In  general,  r(r,t)  can  be  studied  in  termB''of  the 
properties  of  G(r,t).  The  latter  is  more  familiarly 
known  as  Van  Hove’s  pair  correlation  function.  In  fact 
if  one  makes  the  approximation  a^=aj=a,  then 

f  (r  ,t )  a2G(r  ,t ) 


§4.2  Application  to  Light  Scattering 


The  following  derivation  is  due  to  Komarov  and 
Fisher  [2 o]  . 

Consider  a  plane  polarized  monochromatic  beam 
of  light  with  electric  vector 


( ft  t-K 
o  o 


(15) 


which  is  incident  on  a  fluid.  As  we  are  only  interested 
in  dipole  scattering,  the ■ interact ion  of  the  incident 
light  and  the  fluid  molecules  will  result  in  the  induc¬ 
tion  of  a  dipole  moment  density  given  by 

p(r,t)  =  a  2  (r,t)  Z  6(r-r.(t))  ,  (16) 

^  -L 

where  ag  is  the  effective  polarizability  of  each  molecule, 
and  r^(t)  is  the  position  vector  of  the  ith  molecule  at 
time  t.  The  Hertz  vector  is  then  given  by 


. 

. 
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V2n  (r,t)  -  ^2  ~^~2  n(r,t)  =  -4Trp(r,t) 

c  3  z 


(17) 


and  the  scattered  wave  has  an  electric  Sector  S»(r,t) 


given  by 


2' (r,t)  =  grad  div  ff  -  V2  S 


(18) 


In  the  wave  zone  (i.e.  for  large  R),  (18)  becomes 


S’(Rnt)  =  grad  div  n  -  -^-p  fi  . 

c"  3  z 


(19) 


Now  the  solution  of  (17)  is  readily  obtained  in  terms 
of  the  retarded  potential  as 


OO 


ff(S,t)  =  /<Nr '  /dt'  P(r'  V  '  )  6  ( t ' -t  + 

| r ' -S | 


-ii 


)  .  (20) 


—  00 


Using  (16)  and  the  integral  representation  of  the  6- 
function  in  (20),  we  have 


a 


n(S,t)  =  2^/d3r'dt 'd W 


Jo  i  ( n  t  ’  -k  .  r  f  ) 
e  o  o 


-si 


ifi' (t'-t+ 


zEl 


x  e 


E  6 (r-r . (t) ) 

•  _L 

l 


(21) 


For  large  R, 


(grad  div  -  ^p-  -^-p) 

c  3  t 


o 


ifl  1  ( t  ’  -  t  + 


-g| 


_  S,  i£J '  ( t ' 

i-(£1-U_S)e 


R  S.r ' 


t  +  —  —  t-, 

c  cR 


c  R 


R 


(22) 


I 


as  IslonevO',;  i  ofnatf  **  arid  o 
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Hence  (19)  gives 

£'(5,t)  =  a-e'2-  -V-  S)/d3?’dt,d^’2eif^ot  T_Ko*r) 

2ttc  R  R 


R  R.rf 


ifi'  (t ’  -  t  +  H-  -  D  , 
x  e  c  E6(rf-r.(t)). 

i  1 


(23) 


The  Fourier  transform  of  £f(R,t)  is  defined  by 

oo  . 

£'(fi,t)  =  /  dfi'  e1^  t  2'(S,SJ')  . 


(24) 


-OO 


Hence 

2 

£'(8,fi)  =  — (£-. - -  §0/d3r  fdt 1  e1  (^ot  Ko,r  ^ 

2ttc  R  R 


x  e 


-ift ( t f  + 


,  +r  _  i 


c  cR  y  Z  6 (r  f -r . ( t ) ) 

•  ZL 

1 


(25) 


S'tfi.fl)  =  a£p  (S,  -  R)/d3r"dt"e_1(!:2ot"+l<:o-r") 

2irc  R  R 


m(t"+  S.  +  ■■ 

x  e  Z  6  (r,f+r  .  ( t”  ) )  . 

j  J 

Therefore , 

2  4 

Q  01  N  p  p  o .  o . 

<|£(R,ft)|  >  =  — 2  "4'  2  Eo  siri  ’  dtdt" 

4tt  c  R^ 


(26) 


x  e”iDKo“  cR'),r+(^ 


X  N  1  Z  <  6  (r '  -r .  (t+t"  ))6  (?-r '  +?.  ( t"  ) )  >  , 

i,j  1  J 


(27) 
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where  we  have  written  r=r'+f"  and  t=t’-t”,  and 


4>  =  cos 


-i 

|  S,  |  R 


The  last  factor  with  the  double ' summation  is  the  Van 
Hove’s  two-particle  correlation  function  G(r,t)  given 
in  (14).  Hence  we  have,  for  the  intensity. 


1(8, fl) 


Na  fl  o  o  .  r  / 

e  T  •  2  .  r  ,  3  — 1  (.  K 

- I  sin  d)/d  r/dt  e  L  o 

n  2^2  o 
2ttc  R 


flR 

cR 


)  ,r+(fl-fl  )t] 
G(r,t)  , 


(28) 


where  we  have  made  use  of  the  relation 


I(S,£J)=lim  £<  |E(K,fi)|2>  . 


(29) 


T+oo 


We  now  define  a)=fl-fl 


k  =  rn  —  k  =  k  —  k_  , 
cR  o  o  5 


(30) 


^  O  "D 

where  k  =  — rr  is  the  wave  vector  of  the  scattered  wave. 
cR 

Therefore,  (28)  becomes 

2  4 
Na  fl 

=  I(R,ai)  =  -  e  0 

2ttc  'R 


— 2"  I0  sin(^  s(k>a))  j 


(3D 


where 

S(k,u>)  =  /  d^r/dt  e_i((A)t“k*^)  Q(r,t)  ,  (32) 

V 

and  where  we  have  written  flQ  for  fl  because  for  light 
scattering  in  fluids,  oj  «  fl  . 


-* 
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:  r  •  i'  •  £  r 


arid  lo  sa u  j.  *.  v  -\ci  <w  9T<rlw 


.  <  i  !■>  :  ■ 
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Now  for  the  light  scattering" problem,  since 
a)  <<  ,  the  time  constants  and  lengths  characterizing 

density  fluctuations  are  much  larger  than  those 
characteristic  of  molecular  scattering' processes . 
Therefore,  only  the  relatively ' long" time'  and  large 
r  limit  of  G(r,t)  is  relevant ; ' For  this" case  G(r,t) 
is  reduced  to  the  autocorrelated  density  function 


G(r,t)  =  N  d^r’<  n(r T -r ,o )n(r ’ ,t ) >  , 


(33) 


where  n(r,t)  is  the  number  density  at  r  and  t.  Also 
from  symmetry  of  isotropic  fluids. 


G(r ,t )  =  G( | r |  ,  t )  . 


(34) 


We  now  let 


n(r,t)  -  n  =  E  n(kf  ,  t)  e“lk  °r 
°  kT 


n(k,o)  =  n(k)  , 


(35) 


where  n  is  the  mean  number  density  of  the  molecules, 
o  a 

The  nQ  part  does  not  contribute  -  to * scattering .  Subs¬ 
tituting  (35)  and  (33)  into  (32),  and  making  use  of 


1  f  i  ( k '  -k"  ) .  r  3?  =  6^  ^ 

y  J  e  a  r  t- 1  j 


k '  ,k' 


(36) 


we  have 


S(k,uj)  =  N-^  f  <  n(-k)n(ic,t)>dt 


(37) 


—  00 


62 


From  (35),  the  inverse  transform  is 


n(k,t)  =  /  d^r  e^k‘r|n(r,t)  -  n 


o 


(38) 


If  we  also  define  the  Fourier  transform  in  time. 


—r 

n(k,w)  = 


/  dt '  e  n( k, t ’ )  , 


(39) 


(37)  becomes 


S(k,oj)  =  N  ■*"<  n(-k)n(k,a))> 


(40) 


The  integrated  intensity  from  (31),  (39)  and  (40)  is 

(41) 


2  4  2 

+  Not  ft  I  V  ?  ^  ^ 

I  =  /I(R,uj)du)  =  — e0° i,Q —  sin  $  <n( -k)n( k) >  . 


p2  4 
R  c 


5  -1 

In  light  scattering,  typically  k  ~  10  cm  ,  so  that 
k-'*'  >>  intermolecular  distance.  For  small  k  and  except 
near  the  critical  point  of  a  fluid. 


n^knT 

^  /tu  oBo 

<  n(-k)n(k)>  =  - 


(42) 


V  B 


T 


Substituting  (42)  into  (4l),  we  see  that  (4l)  is  identical 
with  Einstein’s  expression  (1-19)  provided  we  identify 


4ttoi  n  = 
e  o 


[pOt]  p=p 


n=no  , 


(43) 


which  is  the  same  as  (1-20). 
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Chapter  5 


SPECTRAL  DISTRIBUTION  OF  SCATTERED 
LIGHT  IN  RELAXING  LIQUIDS 

§5.1  Introduction 

In  this  chapter  we  shall  use  the  formalism 

outlined  in  Chapter  4  to  calculate  the' spectral 

distribution  of  light  scattered  by  a  relaxing  liquid. 

Such  calculations  both  for  a  non-relaxing  and  relaxing 

liquids  have  been  made*  previously  by  Mountain  £23,24]. 

In  a,  non-relaxing'*  liquid '  the  only  mechanisms 

for  the  decay  of  density  fluctuations* are ‘thermal 

conductivity  and  the  static  shear  and  bulk  viscosities. 

Mountain  £23j  showed  that  in  the  limit  of  low  frequen- 

4  2 

cies  (w(b  +  'jrl)/P0v  «  1)j  the  spectrum  consists  of 
three  Lorentzians,  one  for  the  Rayleigh  line  whose 

i 

width  depends  on  the  thermal" conduct ivity ,  and  the 
other  two  for  the  two  Brillouin  components.  The 
intensity  ratio  satisfies  the  Landau  Placzek  formula, 
as  it  should. 

In  a  relaxing  liquid  the  density  fluctuations 
can  decay  in  other  modes  than  described  above  because 
of  the  existence  of  internal  degrees  of  freedom.  For 
the  relaxing  liquid  Mountain  £24]  chose  a  model,  in 
which  there  is  essentially  a  frequency  independent 


. 
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(static)  shear  viscosity  and  a  frequency  dependent 

or  relaxing  bulk  viscosity,  corresponding  to  the 

existence  of  a  single  relaxation  time  .  This  gives 

rise  to  a  fourth  component.  Mountain's  model  is 

2 

applicable  as  long  as  r)  is  such  that  mri/pov  «  1. 

If  this  is  not  the  case,  then,  as  is  well  known,  from 
ultrasonic  data  on  liquids  of  high  viscosity  [13,16], 
shear  viscosity  must  also  be  considered  as  a  frequency 

dependent  or  relaxing  quantity.  Indeed,  if *  one 

2 

solves  Mountain's  equations  ~  for  ojn/pov  >  1,  one  finds 
that  there  are  no  Brillouin  peaks  in  the" scattered 
spectrum,  contrary  to  observation  in  liquids  like 
glycerine  [see  Refs.  35,  36]. 

In  the  present  work  we  calculate  the  spectral 
distribution  of  scattered  light  for  the  following 
model  of  a  relaxing  liquid.  We  assume  (1)  that  the 
ordinary  shear  viscosity  is  the  low  frequency  limit 
of  a  shear  relaxation  process  having  a  single  relaxa¬ 
tion  time,  and  (2)  that  the  same "is  true  for  bulk 
relaxation.  For  the  latter  we  explicitly  introduce 
an  ordering  parameter  and  use  thermodynamic  theory 


Mountain  actually  considered  a  non-relaxing  bulk 
viscosity  also;  this  merely  amounts  to  replacing  n 

i 

3 

by  n  +  fj-b  in  the  relevant  equations. 
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of  relaxation  processes  to  calculate  p(k,t),  instead 
of  using  the  time  integral  representation  of  bulk 
viscosity  as  done  by  Mountain.  We  shall  see  later 
that  our  treatment,  though  somewhat  *  lengthier ,  gives 
better  insight  into  the  type  of  bulk  relaxation 
processes  implied  in  Mountain’s  work.  Anticipating 
our  analysis  we  may  mention  that  Mountain’s  represen¬ 
tation  for  the  bulk  viscosity  is  appropriate  for  a 
density  induced  relaxation  process,  rather  than  the 
thermal  relaxation  process  as  stated  by  Mountain 
[24]  +  . 

In  Sec.  5.2  the  exact  expression  for  the 
spectral  intensity  of  light  scattered~by " our- model 
of  the  liquid  is  derived.  In  Sec.  5.3  ,  this  spectral 
distribution  is  approximately  expressed  as  a  sum  of 
five  Lorentzians  plus  two  other  terms  which  are  usually 


In  a  recent  paper,  which  came  to' hand  after  the 
work  described  in  this  chapter  had  been  submitted  for 
publication  [25,26J,  Mountain  \21~\  has  treated  the 
relaxing  bulk  viscosity  by  introducing  a  thermodynamic 
ordering  parameter  just  as  we  do  in  the  present  work. 
Our  formulation,  however,  is  different  from  his.  Our 
work  differs  from  Mountain’s  also  in  that  we  consider 
a  relaxing  shear  viscosity  rather  than  a  static  shear 
viscosity. 
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small.  The  expressions  for  'the 'heights~and  widths  of 
the  Lorentzians  take  particularly  simple  forms  for 
certain  limiting  values  of  the  shear  and  bulk  relaxa¬ 
tion  times;  these  cases  are  discussed  in  subsequent 
sections  of  this  chapter. 

§5.2  Equations  for  Calculating  n(S,t) 
or  p(£,t)  and  S(k,oj) 

We  have  seen  in  Chapter  4  that" the ' expression 
for  the  spectral  intensity  S(k,w)  involves  the  cal¬ 
culation  of  n(5,t)  or  p(lc,t):  this  calculation 
proceeds  along  the  same  general  lines  as  in  the  work 
of  Kadanoff  and  Martin  [281  and  Mountain  [24].  It 
is  based  on  the  suggestion;  due  to 'Landau  and  Placzek 
JjLl]  ,  that  for  a  non-relaxing  liquid  the  decay  of 
density  fluctuations  is  determined  by  the  linearized 
hydrodynamic  equations  and  the 'heat  diffusion  equation. 
In  our  work  we  have  to  start  with  a  generalized  form 
of  the  former  equations  to  incorporate  a  relaxing 
shear  viscosity,  and  to  combine' it  with  the' thermo¬ 
dynamic  equations  appropriate  for  a  fluid  containing 
an  ordering  or  relaxing  parameter.  The  procedure  is 
to  first  calculate  p(k,t)  in  terms  of  an  initial 
fluctuation  p(k)  and  then  to  take  the  average  over 
initial  states  indicated  by  <.  .  .  . 
for  S(k,oj)  . 


>  in  the  expression 
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The  relevant  equations  for  the  calculation  of 
p(5,t)  for  our  model  are 


jMr  +  p  div  u  =  0 
3  t  Ko 


(1) 


(3) 


x  V2  T 


(4) 


!“£■  -  -  L  6  Z  , 


(5) 


where  the  last  equation  is  from' ( 3-6 ) .  We  shall  take 
p,  T  and  Z  to  be  the  independent - variables .  This  set 
has  the  advantage  of  simplifying  our  calculations  since 
the  thermodynamic  fluctuations  in  density  are  indepen¬ 
dent  of  those  in  T  and  Z. 


Eq.  (5)  then  can  be  rewritten  as 


where 


(7) 


We  shall  now  make  plentiful  use  of  the  thermodynamics 
relations  which  we  derived  in  Chapter  3*  Writing 


VP  =  -p  Vp  +  6BtVT  +  ( 


o 


(8) 
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and  eliminating  u  between  (1)  and  (2),  we  have 


3  2P 

9  t2 


B 


T  „2  .  „„  „2m  .  /9  PN  n2n  .  4 


t 


V  p  +  6B  V  T+(~)t  V  Z  + 

P0  T  9  Z  T ,  p  3po  Q 


x  V2(|fr)  dt’ 


/  n'  (t-t' ) 


(9) 


The  equation  of  heat  flow  (4)  can  be  written  as 


X 


P  T 
Mo  o 


2 

V  T  = 


_1 

2 


6BT  ft  + 


_V  il  (3_S}  ,3_Z . 

T  8  t  d  Z  T  p  ^9  t ' 


(10) 


The  procedure  now  is  to  obtain  the'Fourier  (space) 
and  Laplace  (time)  transforms  of  (6),  (9)  and  (10).  Then 
the  resulting  equations  can  be  solved  for  p(K,S),  which 
is  the  Fourier-Laplace  transform  of'Sp(r,t).  Thus  p(k,t) 
or  n(£,t)  can  be  obtained  as  the  inverse  Laplace  transform 
of  p(£,S).  Accordingly,  we  define 


p(K,S)  =  V"1  /  dt  /  d^r|jD(r,t)-p  lelk,r  St ,  Re  S  >  0  . 

0  V  ° 


(11) 

p(ic)  =  V-1  /  [p(r,0)-p  ]  elk,r  d3r  (12) 

V  ° 

and  similar  formulae  for  T(k,S),  T(k)  and  Z(k,S),  Z(k). 

Now  we  note  that  in  the  expression  for  the 
scattered  spectrum  given  by  (4-37),  n(k,t)  occurs  in 
<  n( -k)n( k, t ) > ,  where  the  symbol  <....>  denotes  statis- 

I 

tical  average,  which  is  equivalent  to  time  averaging. 
Owing  to  the  arbitrariness  of  the  time  origin,  we  have 
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j  '  ol  nota.ie-  Jxe  srfcf  nl  3&ti3  scJon  9w  woM 


69 


<|>(t)  =  <  n(  -£)n(£ ,  t )  >  =<  n(-£  ,t ?  )n(£  ,t+t  ’  )  >  .  (13) 

Since  the  liquid  is,  on  the  average,  isotropic,  there- 
fore  <  n(-k)n(k,t)>  = <  n(k)n(-k,t ) > .  It  is  clear  then 
that  by  a  change  of  time  variables  in  (13),  we  have 


4> ( t )  =  <J> ( — t )  .  (14) 

The  even  property  of  <j>(t)  also  requires  that  =  0 

-4.  ~y. 

at  t  =  0.  Anticipating  that  n(k,t)  or  p(k,t)  is  to  be 
used  inside  <  >,  one  can  require  p(£,t)  to  satisfy 


p(k,t )  =  p ( k 


3  p  ( k ,  t )' 


9  t 


t  =  0  0  5 


(15) 


which  is  the  usual  boundary  condition  used  in  such 
calculations  £29_]  .  Hence  one  may  obtain  the  following 
results 

V_1/d3?/e1^,r_St  :*-£•  dt  =  Sp(k,S)  -  p(£) 

V_1/d3r/elk,r_St  3— §■  dt  =  S2p(k,S)  -  Sp  (£) 

3 

V-1/dt/eik'r’"St  V2pd3r  =  -k2p(k,S)  ,  (16) 


with  similar  results  for  T  and  Z. 

We  now  multiply  (6),  (9)  and  (10)  by  V  1e^k’r  ^ 
and  integrate.  One  obtains  the  following  three  equa¬ 


tions  : 
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(18) 


T(k,S)  -  S(^)TjpZ(k,S) 
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(19) 


Eliminating  T(k,S)  and  Z(k,S)  and  making  use  of  thermo¬ 


dynamic  relations  in  Chapter  3>  we  have 


p(ic,S)  =  (vok)_1  p(ic)  ,  (20) 

where  we  have  neglected  terms  containing  T(k)  and  Z(k) 
because  from  fluctuation  theory, <  p(-k)T(k)>  = 

<  p(-k)Z(k)>  =  0  .  Hence 

<  p(-ic)p(k*S)>  =  <p(-£)p(£)>  v-  (21) 

or 


<  n( -&)n(5, S ) >  =  <n(-ic)n(lc) >a(k,S) , 


(22) 
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where 


a(k,S)  =  ( v^k)”1  ||^y 
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Here  v  is  the  sound  velocity  at  infinite  frequency. 
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Now  from  (22) 9 


a(k,S)  =  <"(.-?>( £>S)>  .  f(y) 
<  n(-K)n(ic)>  vo,k  ^ 


Therefore 


a (k  ,t )  - 


<  n(-£)n(K)> 


oo 


-St 


a(k,S)  =  /  a(k,t)  e  dt  , 
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odc.c)  =  y.(-S)n(S»u>)> 
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=  /  e  ia)t  a(k,t)dt  . 
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Since  from  (27),  we  have 


a(k,t)  =  a(k,  -t)  , 
therefore 


(27) 

(28) 


(29) 


(30) 


a(k,a))  =  2  Re  a(k,  S=ia)) 


=  2  Re 


t 


V  kG(^~~) 
o  v  k 


(31) 


4* 

Strictly  speaking,  one  should  write  (31)  as 

tp  t  i^"^ <5  ^ 

a(  k,oj)  =  lim  2  Re  - ,  6  >  0  ,  (32) 

6->0  v  kG(^r-) 

o  v  k 
o 

which  satisfies  the  requirement  on  S  that  Re  S  >  0 .  (32) 

has  to  be  used,  for  example,  in  the  case  of  an  idealized 
fluid  having  no  sound  absorption.  (The  spectral  distribu¬ 
tion  for  this  case  consists  of  just  three  6-functions.) 


> 


t  .  -  !>M  .  = 

-  >tw  / 


— 


:  )  si.;  W  blv  ;  z  mo  1 3 : 1  *'  3©q3  :  J  c f3 

. 

1  Z  r  >  ter1ef?JU  s»i  irfi  ee.n  niriw 

T  >  sic.  tan 30  9^- bo  alrli  10I  noli 
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and 

S(k,a>)  =  V2N~1  <  n(-k)n(k)>a(k,a))  ,  (33) 


where 

<  n(-Ic)n(ic)> 


n2kDT 
o  B  o 


V  B 


T 


(34) 


from  Appendix 

We  shall  for  simplicity  consider  the  spectral 
distribution  to  be  given  by  a(k,w)  of  Eq.  (31).  The 
total  integrated  intensity  is 


OO  00  00 

/a(k,w)du)  =  /  dt  a(k,t)  /  e-iaJt  da) 

—  00  _oo  _oo 


00 

=  27r/a(k,t  )6(t  )dt  =  27ra(k,o)  =  2  tt  ,  (35) 

_00 

where  we  make  use  of  a(k,o)  =  1  from  (31) 


§5.3  Approximate  Determination  of  the 

Spectral  Intensity 

Eq .  (31)  is  the  exact  expression  for  our  model. 

Since  F(y)  and  G(y)  are  quite  complicated,  it  is  not  easy 
to  visualize  from  (31)  the  variation  of  a(k,oj)  with  a)  and 
other  parameters  defined  in  F(y)  and  G(y).  A  more  con¬ 
venient,  though  necessarily  approximate,  form  of  cr(k,u)) 
can  be  obtained  which  separates  out  the  various  components. 
The  method  is  as  follows: 

The  inverse  Laplace  transform,  a(k,t),  of  a(k,S) 
is  first  obtained  by  writing  a(k,S)  as  partial  fractions. 


« 


> 


D  C  J  O  C  i'B  IrlB  V  ^rfct  (x  j  )  ..10*1  X  9.x  X  i  -  'If  od 

. 
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Then  a(k,w)  can  be  evaluated  from  a(k,t).  We  have 

oo 

a(k,t)  =  /a(k,t f ) 6 ( t-t ' )  dt ’ 

—  OO 


OO  00 


=  5-  /  dt’  /  du) '  a  ( k ,  t f  ) 

cl  TT  * 


—  oo 


—  oo 


00  #  00  #  oo 

=  /  dw'e^  ^{/a(k5t’)e  ^  dt 1  +  /cr(k,t f  )e^w  ^  dt f  } 


—  00 


0 


0 


00 

=  ^  f  da)’  e^w  ^  Re  a(k,iu)f)  . 


(36) 


—  oo 


In  evaluating  (36)  by  contour  integration,  one  needs  to 
know  the  poles  of  a(k,S),  hence  the  roots  of  G(y)=0  and 

i 

this  can  be  done  only  approximately. 

We  shall  assume  that  a(=  — ^ — )  «  1.  This  is 

OO 

P  V  cTr 
Ko  o  V 

true  for  almost  all  nonmetallic  liquids  for  the  value  of 

k  (~  10  )  relevant  to  the  problem  of  light  scattering. 

-9  /  \ 

For  example,  for  glycerine  a  -  5  x  10  k,  and  from  (1-32), 

5-1  * 

k  is  at  most  3  *  10  cm  for  scattering  of  visible  light  . 


Throughout  this  paper  we  consider  this  case  (a  «  1) 
only.  It  may  be  mentioned  that  if  a  >  1,  then  Eq.  (19) 
no  longer  adequately  represents  the  energy  transport  and 
one  should  add  on  the  left-hand  side  of  it  a  term  t^-^- 
where  t  is  a  relaxation  time  characteristic  of  phonon-phonon 
collisions.  For  solids  at  low  temperatures  or  for  liquid 
helium  this  term  may  have  to  be  taken  into  consideration; 
see,  M.  Chester,  Phys.  Rev.  131 ,  2013  (1963). 


p  T— 
po  dt 


' 


1 '  f  rr  r  -  :  noo  *  v  1 ; 
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Let  us  further  assume,  for  the  moment,  that 
ax^  <<  1,  aXp  «  1.  This  covers  most  of  the  interesting 
relaxation  region  x.^,  x 2  <  1.  One  then  readily  sees 
from  (25)  that  there  is  one  root • -y  of  G(y)  =  0,  which 
is  much  smaller  than  unity  and  is  approximately  given 

by 


aCV  =  xk 
YC;  "  povocp 


y„  «  l 


O 


(37) 


To  determine  the  other  roots  we  try  to  factorize  G(y) 
in  the  form 


=  (y+e+^-)  (y+0-^-)  (y+q1+q2^'y+qi’q2')  ^y+yo^ 

12  o  o 

where  v,0,  and  q p  are  as  yet  unknowns,  and  where  v 

can  be  interpreted  as  the  velocity  of  the  two  propagating 

vQ 

modes  and  vQk0  their  ■  absorption  per- second,-  2tt(^-)0  is 
the  absorption  per  cycle.  This  will  be  obvious  from  the 
expression  for  a(k,t)  which  we  shall  obtain  later. 

Now  0  is  usually  small  compared  with  unity  as 

7 

absorption  data  from  ultrasonic  i  show.  Hence  in  the 

2 

following  we  shall  neglect  0  compared  to  unity. 

Expanding  (38),  we  obtain 


^  8  01  ci  •  V 
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p 

-  y5+yZl(2q1  +  20+yo)+y3(q^-q2+X2-+40q1+yo(2q1  +  20)) 
^  ^  vo 

+y2  [20(q^-qp)  +  2q-,~  +  y  (q?-q?+ 


L1  ^2 
2 


L1  2  ’  JoVMl  ^2'  2 

v  v 

o  o 


2  2 

+y[(q^-Q2)^7  +yo(20(q2-qp)+2q.1  )] +y^  ( q?-q2  )^o  .  (39) 

vo 


■1  ^2  H1  2^+yo^ql  q2'  2  * 
v  ,  v 

o  o 


On  the  other  hand,  from  (25) 


X1X2 


■x,xn  2 
L  d  2  v 

o 


+  +a(—  +  — — ~ — )~ ' 


x. 


rB 


00 


00 


+  y 


i|  00 

_S  1_  Vf  +  ac 
lBs  x2  Bsxd  x2xdcy 


,  00  li  00  . 

y  ot(B,p+-jG  )■ 


B 


S 


CVXd 


+  y 


00  i|  00  00 

Brj-jOt  B  ,j, 4" ot  c  y 

+  — -  + 


L.xdx2  BSx2 


B, 


cVXdJ 


;1 


oo 


ac 


V 


cvYX2xd 


(40) 


where  we  have  used  (3-34)  to  define 


00 


'V 


x  ,  =  v  k  t mTT  =  —  X-, 
d  o  TV  cv  1 


(41) 


Comparing  now  terms  Independent  of  y  and  the  coefficients 

of  y^  In  (39)  and  (40),  we  have 

2 

2  2  _  (  s-1  vo 

Q  2_  ^2  ■  x2Xd  2 


(42) 


2Q1  =  x"1  +  x-1  -  26  +  c(l  - 


-1  ,  -1 
"  x2  +  xd 


(43) 


20 
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Hence 


2q 


2  "  L'a2 


[(Xo1  +  xd  1  -  20 )?-  4(x?xd)  1 

V 


(44) 


2 

Next  equating  coefficients  of  y  in  (39)  and  (40),  making 
use  of  (42)  and  (43),  and  retaining  only  the  linear  terms 
in  a,  we  have 


I,  __  0  rs  (B0+  TrG  ) 

29  =  r?2xd(r-}  fe2(7-)'[(F-)_-  — §-  3'  3 

O  O  O  S 

00  00 

+  x,(^-)2r(M2-  %  +  a(h^  r-  -  (2- )2j 

dvv  '  L>y  y  R  J  v  o__  Lv  Vv  J  —1 


/V  \2  rv v  s 2 

+  ax2xd  {~ >  C— ) 


CV  "Y 

V  ^00  00  v 

(BT+3^ 


V 


O 


o 


B 


S 


3>- 


(45) 


The  only  unknown  left  now  is  v/v  .  This  we  determine 
approximately  by  comparing  the  coefficients  of  y  in  (39) 
and  (40).  Eliminating  q^  and  q^  from.  (42)  and  (43), 

neglecting  terms  containing  a  and  terms  proportional  to 

2  2 
0  ,  one  obtains  for  v  the  equation 

,  00  4  n  00  \ 

f  v  v8  -  n2,  r  v  \6r  2  2,  /  N  2  BS+3°  - 

( (X0X,j  +  ( — )  IXo+X  ,+XnX -,-(X0X  ,  ) 5 — = — — 

v  2  d  v  u  2  d  2  d  2  d  B0  -1 

00  S 


-  (  )  (x2+x  )|>2(Bs  +  3  G  )/V  xdBS/BS-l 

o 


+  C^~)2  &?xd(Bg  +  |  G°°)/Bq  +  1]  -1  =  0 
vo 


(46) 


We  are  interested  in  that  solution  of  (46)  which  gives 
2  2 

v  =  vQ  when  xp=xd=0o  A  particularly  simple,  though 


. 
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2 

approximate,  expression  for  v  is  obtained  if  one  sets 

(^—)B  =  1  +  A,  and  assumes  that  A  <<^1 0  Then^retaining 
o  ^G 

only  the  linear  terms  in  A  and  in  -  -D —  and  ~  one  may 

Bs  Bs 

readily  obtain 
2 


V  -Id. 

— p-  -  1  +  A 
vo 


4  00  2  oo  2 

!  +  3^_  X2  +  BS“BS  Xd 


BS  1+Xp 


BS  1+Xj 
d 


(47) 


Eq0  (47)  gives  correct  values  for  the  velocities  in  the 


various  limits,  namely  v  -*■  v  for  xp, 

_00  4^00  0  (- 


X 


0  and 


v  -*  v 


00 


b; 


S 

for  x,,Xi  +  v  ( — )  2  for 
2 5  d  5  p 

d  Ko 

4  G00  ^ 

x0  ■>  0,  x,  -*■  oo  and  v  (•=•  — )  2  for  x,  0,  xn  ■* 
2  5  d  ip  d  5  2 


00  , 


For  intermediate  values  of  xp  and  x^ ,  it  is  only 

approximate .  If  on  the  other  hand" we- determine  the 

velocity  without  neglect ing" the ■ terms  proportional  to 
2 

0  ,  we  find,  on  ignoring ‘ terms' proportional  to  a  and 
retaining  terms  only  up  to  quadratic  powers  in  xp  and 


x ,  that 

d 

vf.  ~  !  +  BS~BS  (1  1  BS~BS,  2  +  3°  (1  .if!) 

v2  "  BS  "  7  BS  d  BS  3  BS 

O 


oo  (B  -B  ) 

2  G_  v  S 

3  Bs  Bs 


X2Xd  + 


o  o  c  o  o  c  o 


(48) 


Finally,  to  complete  the  term  by  term  comparison  of 
(39)  and  (40),  we  note  that  the  coefficients  of  y  in 
them  are  equal  if  terms  proportional  to  a  are  neglected. 
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Having  determined  the  roots  of  G(y)=0  approxi¬ 
mately,  we  are  now  in  a  position  to  calculate  a(k,t)„ 
We  first  write  a(k,S)  as  partial  fractions 


a(k,S) 


2 

Z 


a. 


m 


aB+ig 


aB-1^ 


n  S+S 
m=0  m 


S+F+ivk  S+T-ivk 


(49) 


where 


r  =  v  k  0  , 

o  5 


S  =  v  ky 


J£i£ 


2 


o 


O  "  o  pc 


o  P 


S1  v0k  ^ql+q2^J  "2  vo“  VM1  M2 


S0  =  v  k  (q. 


q0) 


(50) 


where  a  ,  o] }Og  and  g  are  real  constants  (independent  of 
S),  which  are  to  be  determined .  Now  using  (36)  and  (49), 


we  have 


a(k,t)  =  E  a  e  ^rrr  +  2e  •  ^  '  (aRcosvkt  +  g  sin  vk|t|) 
m=0  m  a 


(51) 


In  the  above  equation,  the  interpretations  we  have  given 
to  v  and  r  become  clear . 

In  the  calculation  of  o(k,t)  by  contour  integration 
the  singularities  which  are  connected ~ with' the • root s  of 
G(y)=0  all  lie  in  the  upper  complex  plane .  The  quantities 

a  ,  °B  and  g  are  Slven 

-S 

p(: — 


=  lim  ( S  +  S  ) 
S+-S 

m 


v  k 
o 


v  kGC—  f) 
o  •  v  k 
o 


(52) 


0  . 
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on  =  Re  lim 


B 


( s+  r+ivk ) 


s^-r-ivk 


Tn/-r-iVkN 

F(~“k~-) 

_ o 

i n / -r-ivkv 
o  v  k  ' 


F( 


=  Im  lira  ,  (S+r+ivk) 
S^-r-ivk 


o 

-T-ivk 

v  k 
o 


) 


i n / -r-ivk \ 
v  kG  ( -- — .  - — ) 
o  v  k 
o 


(53) 


(54) 


These  expressions  then  give 


°o  = 


1  -  - 
Y 


a 


2  F(q 15^2)5  o 2  f ( q  ^ i  Q  2  ^  > 


aB 


albl+a2b2 

2  2 
bl+b2 


£l  -j^b  p  — a  p  b  -j 

2  2~~ 

bl+b2 


(55) 


where 


f(q-|  ,q  ?)  =  '  ((1-Xp(q  +qp) )  [1  -  i  +(i  -  ~)(q1+qp)x  ] 


B, 


OO 


■1,H2 


2VM1  M2 


Y  •  y  Bg  ^2  d- 


OO 


+  (q1+q2)(l-xd(q  +q  )  [9  +q  -x  (q  +q  )2-  j  §-  } 

S 


■C2q2^i+q2^xd^2  Cf  ql+q2^ 


2 


20 (qi +q? )  +  ^  * 

V 


V_jr  q  -1 

o 


(56) 


2  2  00- 

1  .  i  _  v_  +  x  x  r^-)2bv ~-Y  Ml 

Y  2  d  V  l-  v2  Yj 

o  o 


OO 


YJ  +  e£xd(3f  -  fe) 

v  v  S’ 

o  o 


+  xn( 


3v  ^  G  1  y~k  ^  2  p^V^2^/V\2  1"|  r  n  <-r  \ 

2'—  3  bT  "1+7)]+e  x2xd  B— }  “6(~}  "  7-  (57) 

v  ~  S  '  ^  ^ 


o 


o  o 

2 


a  =  )!_  (x  (is  _  1  _  vfn  +  y  M  _  I  _  Zl  +  £  G_s 

a2  v  1  d  Bs  y  2 '  2  y  v2  3  Bq  ■ 

vo  o 


o  2 

2  00  _ 

+  20Q-I  +  XpXg(y  +  ~~V~  "pT~ - )1  30  "  (xo+X^  )  } 


V 


2  d 


o 


(58) 


,  _  -n  /V\2p/0\2/l  \  v 

bl  x2xd  v  Dv  x5x,  -  2 

o  2  d  v 


-  'tqqel 


(59) 
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x0x  , 
2  d 


1 


(60) 


Prom  Eq,  (29)  and  Eq,  (51),  one  then  has 


2  a  2S 

a(k,aj)  -  E  — ~ — p  +  a 
m=0  a)  +S 

m 


2r  _ _  +  2T 

B  Ku3+vk)24T2  (a3-vk)2  +  r2 


(61) 


§5«^  General  Discussion 


Eq,  ( 6 1 )  is  an  approximate  formula  for  the  spectral 
distribution  a(k,uj).  It  is  correct;  as  long' as  the  sound 


2 


absorption  is  small  enough  so  that  0 


<<  1  Is  justified 


and  that  x?}  x^  are  not  too  large  so  that  ax9,axd  «  1  holdc 
The  first  three  terms' in  (6l)  each  has  a  peak  at 
a)=0,  i  <,  e  „  at  the  incident  light  frequency0  Thus  the 
central  or  Rayleigh  component  for  our  model  consists  in 
general  of  the  sum  of  three  LorentzianSc  The  width  of 
the  first  of  these  is  determined  by  the" thermal  conduc¬ 
tivity,  whereas  the  widths  of  the  other  two  depend  on  the 
relaxation  ,times  and  t ^ ,  in  a  way  showing  that  the 
shear  and  bulk  relaxation  processes  are  mixed  together. 

The  next  two  terms  in  (6l)  have  a  peak  at  oj=-vk  and  a)=+vk 
respectively  and  represent  Brillouin  components.  The 
width  f  of  each  of  them  is  the  sound  absorption  per  second. 


* 
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Thus  if  the  last  two  terms  in  ( 6 1 )  could  be  ignored, 
the  spectral  distribution  'would  consist  of  a  sum  of 
five  Lorentzians,  three  peaked  at  the  incident- 
frequency  and  the  other  two  at  w-±vk  respectively. 

We  observe  that  the  last  two  terms  in  ( 6 1 ) 
contribute  nothing  to  the  integrated  intensity.  Also, 
as  may  be  readily  verified,  their  effect  on  intensity 
distribution  may  be  neglected ■ so  long  as  0(=r/vk)  <<  1. 
If,  however,  0^1,  they  af feet ‘ markedly  the  position 
and  shape  of  the  Brillouin  peaks.  It  may  be  mentioned 
that  these  asymmetric ' terms* would  occur  also  in  the 
corresponding  expressions- in  the  works  of  previous 
authors  [24,  28],  although  they  did  not  write  them 

i 

* 

explicitly  .  The  origin  of  these  terms  can  be  traced 
to  the  requirement  that' any  solution' for  a(k,t)  must 
satisfy  the  boundary  conditions  similar  to  (19),  viz. 

a(k,t)  =  a(k,-t)  and  '-\  =  0  .  (62) 

3t  j  t  =  0 

The  second  of  these  conditions  can  be  satisfied  only 
if  g  is  not  approximated  to  zero  in  eqns.  (49)  and  (51). 

j 

Using  (62)  into  (51)  one  obtains  a  consistency  relation 
for  the  solution 


The  existence  of  these  terms  has  been  recently  noted 
also  by  Montrose  et  al.  [30]. 


(63) 


2 

2gvk  =  E  am  Sm  +  2anT 
n  m  m  B 

m=0 

The  importance  of  the  asymmetric  terms  in  the  spectral 

i 

distribution  is  probably  best  illustrated  by  consider¬ 
ing  a  model  in  which  the  effect  of  thermal  conductivity 
can  be  ignored  and  the  only  cause  of  sound  absorption 
is  due  to  non-relaxing  viscosities.  One  can  obtain  a 
simple  exact  expression  for  a(k,m)  for  this  case  which 
we  shall  discuss  at  the  next  section. 

For  the  total  integrated  intensity,  we  have 
00 

from  (35)  that  /  a(,m)  dm  =  2tt.  If  we  perform  this 

—  00 

integration  on  Eq,  (6l),  we  obtain  the  relation 


ao  +  °1  +  °2  +  2aB  =  1  • 


The 

R  = 


intensity  ratio 


■R 


c  +an+a0 
o  1  2 


21, 


2o 


B 


is  then  given  by 


~  °B 


(64) 


(65) 


This  quantity  is  in  general  very  complicated  and  differs 
from  the  value  of  y-1  of  Landau  Placzek. 

We  shall  now  study  some  special  cases  for  which 
the  various  complicated  expressions* simplify  somewhat 
or  very  much.  Moreover  we  shall  show  that  for  some 

j 

cases  our  results  can  be  reduced  to  those  obtained  by 


previous  workers. 


.ti  )  'lev  I  Bier,  "j  3  nJ  \  1  fneup  alriT 


.  '  r  ■  "  ••  u  •  \  ->J J-*  I  . I  1  :  2  :  V 
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§5°5  Special  Cases 

( a )  Unrelaxed  viscosities  as  the  only  loss  mechanisms 

This  case  is  ob’tained  if  we  put  a=0,  x^,Xp  -*  0, 
but  (Bc-  B0)x,  =  v  kb.  G  x0=  v  kp  are  finite.  Then  for 

O  O  Q  O  cl  O 

1  4  2 

small  dispersion,  r  ■>  p —  (—  n  +  b)k  ,  a  =o^  =  op  =  0, 

~  po  p 

aD  4  g  hT/v  k0  However,  one  can  easily  carry  out 
an  exact  analysis  for  this  case.  The  factorization  of 
G(y)  can  be  done  without  approximation .  Prom  (23),  (24) 
and  (25)  j  we  have 


c(k,S) 


1  F(y) 

vQk  gT77 


F(y)  =  y  +  a 


3 


G(y)=y+ay+l 
4 


a  = 


k(b  +  j  n) 

P  v, 

^o  k 


(66) 


Then 

G(y) 

with 


(y  +  e  +  —■) ( y  +  e  - 

O  O 


v_ 

v 


=  (1  -  if-  a2)h 


o 


0  -  ^a , 


or  T  =  4 


4  ,  2 

(b  +  T  p)k 


(67) 


o 


From  these  one  determines  easily  that 


a  =a-i=a0=  0  , 
o  1  2  5 


a  = 

UB  2 


-  x. 


g 


r/vk 


(68) 
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All  these  are  exact  results.,  We  see  that  g  is  not  zero, 
so  that  the  Brillouin  peaks  do  not  occur  at  co=±vk» 

Their  position  can  be  obtained  by  finding  the  maxima 
of  the  spectral  distribution 


a(k,co) 


x 


'2 


f  2  r  +  2  r  1 

L(a)+v.k)2  +  r2  (w-vk)2  +  r2J 


r 


vk 


co  +  vk 


co  -  vk 


*-(  co+vk ) 2  +  f2  (co-vk)2  +  f2- 


(69) 


Prom  (31),  this  is  also  equal  to 
a(k,co)  = 


2a  v3  k3 
o 


2  2  2,2  /  2  2,2x2 

a  co  v  k  +  ( co  -v  k  ) 
o  o 


9  o  ( k,co)  r,  •  ,  •  jr. 

Then  — -^-2 — -  =  0  gives  two  maxima  at 
9  co  ° 


(70) 


CO  =  CO 


B 


+vok  [l  -  3sa2]" 


(71) 


and  a  minimum  at  co=0„  So  in  general,  coB^vQko  We  also 
see  in  Eq .  (71)  that  if  the  viscosities  are  so  large 


that  a  >  2  or 


( b  +  2-n  )  k 


P  v 
Ho  o 


>  /  2 


(72) 


then  there  are  no  Brillouin  peaks. 


(b)  x2  <<  1,  ax^  <<  1 

This  case  is  relevant  to  the  study  of  a  relaxing 
bulk  viscosity  and  we  should  then  be  in  a  position  to 
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compare  our  model  with  that " of " Mountain T s 0 

In  this  case,  since  a,  «  1,  we  have  also  ax p  <<  1 
(50)  and  (55)  then  reduce  to 

2 


v 


o 


I  T  5 

n 


2  2  3 
v  xd 


0  P0CP 

aQ  =  1  -  i  ,  a1  ~  0  ( Xg )  ~  0  , 

(B  -B  )k  2  v-4  p 

°2  {—0 - (V  -  +  vok  (1  -  7^} 

M0  V  V  T  Y 

o  d 


x  (v2k2  + 


V 


T2 

v  Td 


o  }-l 


2aB  1  °o  al  °2 


(73) 


From  (45),  we  have 


00 


2  B„  2 

20  “  ^2  Xd  <BT  -  V 

v  S  v 

o  o 


+ 


00 

ac„  2 
V  ,  v 

'•  Q 


•v  v 


o 


(74) 


We  have  neglected  any  static  viscosities  in  our  model. 

If  one  wants  to  include  the  static  shear  viscosity  also  , 


c 

one  can  do  this  by  taking  Xp  *  0  while  G 


and 


00 


G  Xp  =  jn  vQk0  Thus  (45)  becomes 


The  frequency  independent  part  of  the  bulk  viscosity, 
b  ,  as  considered  also  by  Mountain,  can  be  included  easily 
by  replacing  n  in  (75)  by  n  +  b  . 


'»  o  3  rid  J 
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00 


2  Bc 

29  ,  ^xd(gi 

vo  S 


00 

2  acTr  2 

I_)  +  _Z(V_ 

2 J  c  v  2 

v  V  v 

o  p 


1 

Y 


4  .  2 

+  ^ - 5-(l 

P  V° 

o  o 


axd } 


(75) 


For  the  special  case  0,  it  can  be  verified 

easily  that  Op  0  and  2gb  =  ^  „  (76) 

Hence  there  is  just  one  mode  for  the  central  component, 
and  the  intensity  ratio 


R 


Y  -  1 


(77) 


corresponds  to  the  Landau-Placzek  value.  The  sound 
absorption  per  cycle  is  from  (75)  given  by 


or-  r<  » °°?  2 a-]  .  2  4  ,  / 

2f  =  Dv  -  vQ)Jk  Td  +  3  p —  +  (Y  - 


00 


00 


where  vf  =  (BQ./p  )  2  is  the  velocity 

o  O 

frequencies  such  that  Xp  <<  1  and 

For  the  case  x,  :>>  1,  we  have 

d  5 


a0  -  —  (1  - 
2  Y 


v' 


2 
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00 
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2aB  1  °o  a2 


2 

,  v" 
1  o 


B 


T 

00 


Y  V 1 “2  Bs 


,  2 

1)  ,  (78) 

DC  5  vi/ 

po  P 

of  sound  waves  at 
>>  1  0 


(79) 


Therefore 


R 


Y 


2 

f00 

“2~ 

v 


o 

Eq.  (80)  was  first 
temperature -induced 


1 

obtained  by 
relaxation 


Rytov  Til]  . 
process,  B^ 


(80) 

For  a  purely 

oo  r 

=  Bt  [see 


. 


r 

to  vbcm  a  o  1  L  t  9*3  arid1  or  ■ 
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Seco  3-3c | ,  and  Eq.  (80)  may  be  written  in  the  alter¬ 
native  form 


R  = 


00 

Y 


1 


(81) 


Let  us  now  compare  our  results  with  those  obtained 
by  Mountain  [24].  The  expressions  (6l)  with  g=0  (which 
Mountain  assumed)  and  SQ,  S^,  a  ,  a1  etc.  given  by  (73) 
agrees  with  the  expressions  (43)  in  Mountain 1 s  work 


provided  we  identify  his  t  with  our  5  Tsv°  ^-n  Sec  .3-2, 
we  have  discussed  and  given  the  transformations’  of  the 
various  types  of  relaxation  times;  We  wish  to  point  out 
that  it  is  important  to  remember  which  relaxation  time 
is  determined  from  the  light  scattering  data,  if  it  is 
to  be  compared  with  that  obtained  from  ultrasonic  data. 

Now  Mountain  in  his  work  represented  the ' relaxing 
bulk  viscosity  by  a  time  integral.  As  we  remarked  before, 
this  representation  does  not ■ take • into  *  account  properly 
the  interaction  between  effects  of  heat  flow  and  relaxa¬ 
tion.  In  fact,  if  ax^  is  not  small  compared ' with  unity, 
then  our  expression  for  a(k,w)  would  be  different  from 
his.  This  is  best  illustrated  by  neglecting  in  our  exact 
expressions  (24)  and  (25)  for  F(y)  and  G(y)  the  terms  due 

00  V 

to  shear  processes  (i.e.  G  =0,  x?=0)  and  comparing  the 
a(k,S)  so  obtained  with  the  corresponding  expression  for 
a(k,S)  obtained  from  Mountain’s  paper  [Eqs.  (28)  and  (29) 


.If  :  .  it  fi  yd  fJ  io:  v  y  £ud 

1  £&  9d  M l  v  (u  l’A)x>  I'-'l  no.  iS9riqxe  ijjo  nsrid 
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of  Refo  24  with  bQ  =  oj .  We  find  that  the  two  expressions 
differ  from  each  other  wherever  the  product  ax^  occurs. 
The  two  become  identical  only  if 

00 

Cy  =  Cy  ,  (82) 

or  equivalently. 


BS  “  BS  =  BT  ”  BT  5  (83) 

that  is,  for  a  density  induced  relaxation  process. 

We  shall  not  give  here  the ' expression  for  a( k,w) 
valid  when  xp  «  1  and  ax^  >  1,  because  it  is  only  in 
liquids  of  high  viscosity  that  the  second  condition 
can  usually  be  satisfied  for  values  of  k  relevant  to 
light  scattering,  and  in  these  liquids  then  ax^  >  1 
also  [l3,  16J  .  For  the  case  G°°=  0  (i.e.'no  shear 
viscosity)  and  ax^  >>  1,  one  may  obtain  o(k,u))  from  the 
expressions  given  under  case  (d)  below  by  putting  G  =0„ 


(c)  x2>xd  >:>  ax?,axd  <<  1 


This  represents  the  high  frequency  limit  of  the 
expression  ( 6 1 )  for  a(k,ai),  for  it  is  by  assuming  the 
second  condition  that  ( 6 1 )  was  derived. 


For  this  case,  Eqs.  (50),  (55)-(60)  give 
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o  d  T) 


°1  =  f(Q1»q2)  -  [2q2(qi+q2)VQkS,Tdtr^] 
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x  {[l-vokTri(q1+q2)]  [(Bs-Bt)/(pov“  ) 
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-  (q1+q2)v0icTd  _  2 
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Q 2  =  fU15  -q 2)» 
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a  =  1  -  — 
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2aB 
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00  4  00 
Bq-Bq 

2r  =  - - 2 -  +  (YW-D 


Bq  .  2 

00  -  *  S  yk 


(84) 


00 


00 


p  V  T ,  P  V  T 

o  d  o  n 


00  o  c°° 
p  V  P0^p 


Of  course,  v  in  the  formula  ( 6 1 )  is  now  v  .  The  inten¬ 
sity  ratio  is 


R 


ao+al+a2 


2a 


B 


00 

1  ,  ,  v  v  2  , 

2^-  -  1  =  Y(— )  -  1  , 

B  o 


(85) 


which  has  the  same  form  as  (80),  except  that 

4  00 

.  . .  ..-in.  -  x, 

oo  oo  .  oo 

vf  =(BQ/p  )  is  replaced  by  v  = 

o  O 


B^+^rG 


o 


The  expression  for  ,  ap,  S^,  also  take  simple 
forms  if  one  of  the  relaxation  times  is  much  greater  than 


the  other. 
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For  t >>  t 
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Sp  *  — 

2  „oo 


BsTd 
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(86) 


On  the  other  hand,  for  j  >>  x^, 

4  o° 

BS 

S-.  ^  - r. — -  ,  a,  -  — = - T^— 

1  /  .  H  ^  oo  \  _  5  1  ,  „  .  4  -  °0 


(BsT  )Tn 


y(Bg+^G  ) 
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4  00 

BS+  3° 

pov“"Td 


Bm (Bn  -  B„) 

2  po^^vK5 


(87) 


(d)  xpixd  >>  1>  axQ,axH  >>  1 


2 5  d 


We  still  assume,  however,  that  a  <<1.  For  this 

case  the  results  obtained  In  ( 6 1 )  etc.  no  longer  hold, 

and  one  has  to  factorize  G(y)  again.  It  can  be  verified 

?  . 

easily  that  the  small  root  -y  of  the  equation  G(y)  =  0 
is  given  by 


OO  4  OO  00  4  ~  00 

»  BT+3^  Brp+^G 

y  «  - —  a  =  - — 

U  pooL  _  ,_oo^- 
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00 
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P  C,rV 
Ko  V  o 


(88) 


instead  of  yQ  given  by  (37)°  Then  writing  G(y),  analo¬ 
gously  to  (38),  in  the  form 


( b) 
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=  (y+0+t^Hy+e_  (y+q|+q2^y+qi_q2^y+y~)  (Q9) 

12  o  o 


o 


and  following  the  procedure  similar  to  that  In  obtain¬ 
ing  (6l),  we  have 

bm 


l|  p(Brp+^-G  )  Dpp  -1 
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(w+v°°k)2+r2  (a)-v”k)2+r2  I , 


(91) 


where  x-,  =  Xmy  ,  and  we  have  neglected  g,  which  Is  of 
the  order  a;  T  Is  again  given  by  (84),  and 


sl  =  v^ky^  ,  Sj  =  vok^ql+q2 )  >  Sp  =  v0k(qi“q2^ 
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00 

The  Intensity  ratio  Is  y(~)~  -  1  as  in  (85).  For  the 

vo 

case  xd  >>  x^  , 


r,'w  •  (sP-  •>>?  •  j 
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and  for  x  >>  x,  . 

n  d  3 
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(94) 


A  comparison  of  (91)-(94)  with  the  expression  for 
a(k,aj)  for  the  case  (c)  shows  that  both  the  intensity 
distribution  and  the  integrated,  intensity  in  the  Brillouin 
components  are  the  same  for  case  (d)  as  for  case  (c). 

There  is,  however,  a  redistribution  of  intensity  amongst 
the  three  Lorentzians  which  contribute  to  the  central 
peak  as  one  goes  from  the  case  axp,ax^  «  1  to 
aXp,axd  >>  1.  However,  since  the' total  integrated 
intensity  of  the  Rayleigh  component  remains  the  same, 
it  would  probably  be  hard  experimentally  to  distinguish 
the  two  cases.  Another  remark  we  wish  to  make  is  that 
the  intensity  distributions  in  the  Lorentzians  associated 
with  the  relaxation  processes  are  simply  related  for  the 
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two  cases  (c)  and  (d).  For  example,  Eqs.  (93)  and 
(94)  may  be  obtained  from  the  corresponding  Eqs.  (86) 
and  (87)  by  replacing  in  the  latter,  B<,  BT,  B~  ■>  B^ 
and  T^(=T3y)  showing  that  the  relaxation 

modes  are  now  decaying  isothermally ,  in  accordance 
with  the  condition  ax^ax^  >>  1.  The  Brillouin  com¬ 
ponents,  on  the  other  hand,  remain  unaffected,  since 
we  have  still  retained  the  condition  a  <<  1. 

§5.6  Intensity  Distribution  calculated  from 
Thermodynamic  Fluctuation  Theory 


The  total  integrated  intensity  for  the  scattered 

light  from  liquids  by  density  fluctuations 'has  been  shown 

2 

in  Chapter  1  to  be  proportional  to<  6p  >.  From  thermo¬ 
dynamic  theory  of  fluctuations  Qsee  Appendix!  ,  we  have 


0  pJk„T 
„  x  2^  Ho  B  o 

<  6p  >  = 


B 


(95) 


T 


This  value  is  independent  of  the  possible  kinds  of  loss 
mechanisms  which  are  responsible  for  sound " absorption 
in  the  liquids.  The  flatter  only  have  to*  do  with  the 
distribution  of  intensity  amongst  the  different  compo¬ 
nents  of  the  spectrum.  For  a  model  of  fluid  for  which 
the  only  source  of  sound  absorption  is  thermal  conduc¬ 
tivity,  the  distribution  of  the  integrated  intensity 
between  the  different  components  can  be  calculated  using 


i‘  '10  VJ+O.r/ 


e.i  u:o  i3i  >qtnco  Jnf>‘  9‘i1JL£>  ©rid  ne9wd9d 
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the  theory  of  thermodynamic  fluctuations.  This  was 
done  in  Chapter  1,  Sec.  1-3D.  1  In  this  case,  the 
central  nonoscillatory  Rayleigh  line  has  a  value  of 
intensity  proportional  to  1-  ^  and  the  oscillatory 
Brillouin  components  have  together  a  value  propor¬ 
tional  to  —  . 

7 

When  there  is  an  additional  cause  of  sound 
absorption  arising  from  a  relaxing ■ bulk  viscosity, 
if  one  represents  it  by  a  time ' integral, ' it  is  not 
built  into  the  thermodynamics' of " the  system;  '  hence 
one  cannot  use  the  thermodynamic  theory  of ■ fluctuations 
to  investigate  the  spectral  composition  of  the  scattered 

i 

light.  In  our  rriodel  we  '  have  used  the  thermodynamic 
ordering  parameter  instead.  We  have  also  developed  in 
the  Appendix  a.  thermodynamic  fluctuation  theory  to 
include  an  ordering  process.  Hence  we  can  use  such 
a  fluctuation  theory  to  calculate  the  spectral  composition. 
However,  since  the  shear  relaxation  process  is  not 
described  in  the  thermodynamic  theory,  our  calculation 

oo 

will  correspond  to  the  case  G  =  x^=  0.  We  shall  now 
give  such  a  calculation. 

First,  taking  S,  P  and  Z  and  independent  variables. 


we  have 


)  sn  cj  U  L:  evnl  oJ 

093  0  9vsri  svi  lebom  tjo  n I 
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From  Appendix 
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Therefore  using  thermodynamic  relations  derived  in  3-2, 
it  can  be  verified  easily  that 
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(99) 


(100) 


(101) 


The  sum  of  (98)  to  (101)  is  unity,  in  accordance  with 
Eq .  ( 64 ) . 

If  we  take  S,  P  and  £  as  independent  variables. 


we  have 
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+  P  S  ,  P  ^3  S'  p,5  <SSS^> 


=  <6 p  >'s  +<6 p2>p  +  «Sp2>^  +  <6 p 2 > g ^  (102) 

2 

where  we  have  used  <  >’  to  denote  average  of  6p  using 

a  different  set  of  variables  i.e.  S,P  and  E, .  From 
Appendix,  we  have 
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In  the  last  equation,  we  have  made  use  of  the  relations 
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and  (<L£)  =  (<LE)  _  (<L£)  f  —  ) 

o  S  ^  P  Z  o  S  ^  P  £*  d  £  S  P  u  Z  ^  S  P  0 


Again  it  can  be  easily  verified  that  the  sum  of  (104)- 
(107)  is  unity. 

In  Brillouin  scattering  we  are  interested  in 
separating  out  the  various  contributions  into  oscilla¬ 
tory  and  nonoscillatory  modes,  that  is,  we  want  to 
write 


2  2  2 
<  6p  >  =  «5p  >  +  <6 p  > 

K  K  nonosc.  K  osc 


(108) 


We  shall  now  study  Eqs.  (96)  and  (102)  in  the  light  of 
such  a  separation. 

For  the  study  of  not  very  high  frequency  case, 
i.e.  for  aX‘£  «  1,  Eq.  (96)  is  the  more  appropriate  one 
to  consider,  because  the  relaxation  process  is  not 

i 

coupled  to  the  entropy  mode,  so  that  the  contribution 

to  the  total  integrated  intensity  normalized  to  unity 

2  1 

from  the  entropy  mode  is  given  by  <6p  >g  =  1  -  —  . 

[From  now  on  we  shall  for  convenience  assume  that 

p  P^k-nT 

<  Sp  >  is  normalized,  i.e.  °w  ■  is  put  equal  to  unity. J 

.Drp 

In  the  low  frequency  limit,  i.e.  for  vokT^  <<  1,  energy 
transfer  can  easily  occur  so  that  the  relaxation  process 
is  completely  coupled  to  the  oscillatory  mode.  Therefore 


<  6P2>osc  E  2oB  ■  «<>%  +  +  <6p2>PZ 


y  •  (109) 


<  q8>  *  <  q6  > 
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These  correspond  to  the  Brillouin  components  and  the 

results  agree  with  those  of  Eq.  (78).  For  high 

frequency  limit  (while  still  retaining  ax^  <<  1),  (109) 

is  no  longer  true  because  energy  transfer  cannot  occur 

2 

easily,  and  <6p  >  has  to  be  evaluated  at  constant 

09  K  osc 

£  and  S.  However,  because  of  ax^  «  1,  we  still  have 
2  l 

<  6p  >g  =  1  -  —  as  contribution  from  the  entropy  mode 

The  oscillatory  component  can  be  extracted  from  the 
2  2  2 

sum  <6p  >p  +  <6p  >2  +  <6p  >p^  as  follows: 

We  write  6 Z  as 


6Z  ^S^p,cfis  +  (3P)S,S<5P  +  65  * 


(110) 


and  then  take  6S=6£=0  so  that  6Z  =  (|-p-)g  ^6P  .  Hence 
we  have 

<6z2>  =  =  h 

po 

<6P6Z>  =  (|-|)s  «SP2>=  -  U-  (f|)s  p<«P2>-  (HD 

po 

Using  (3-26),  (97)  and  (111),  we  have,  in  the  high 
frequency  limit. 


r  2 

<  6  p  > 


osc 


£<6p2>p  +  <6  p  2  >  7  +  <6  p  2  >  p  7~] 


PZJS,5 


—  =  -  —  =  2pr 
°°  y  R°°  B 

Bs  Bs 


(112) 


To  conserve  the  total  intensity,  we  get  an  amount  of 

1  BS 

contribution  equal  to  —(1  -  — 1 — )  for  the  nonoscillatory 

Y  Bs 
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mode  due  to  the  relaxation  process.  These  results  agree 

with  those  of  Eq.  (79). 

For  very  high  frequency  (i.e.  ax^  >>  1),  the 

entropy  mode  is  coupled  to  the  relaxation  process,  hence 

its  contribution  to  the  intensity  is  no  longer  given  by 
2  1 

<  6p  >c  =  1  -  —  .  The  contribution  from  the  oscillatory 
K  S  y  17 

mode  is  still  given  by  (112),  as  expected.  There  is, 
however,  a  redistribution  of  intensity  between  the 
entropy  mode  and  the  relaxation  mode.  In  this  case,  it 
is  more  convenient  for  us  to  consider  Eq.  (102)  instead 
of  Eq.  (96).  Putting  6S=6£=0  in  that  equation,  we  have 


2  p  1  1 

2 cD  =  <6 p  >  =  <6p  >„  =  -  — 

B  M  osc  M  P  Y  R°° 


which  is  the  same  as  (112). 

The  contribution  from  the  entropy  mode  at  this 


high  frequency  limit  can  be  obtained  easily  by  taking 

From 


2  »  P  '  2  f 

6£  =  0  in  the  sum  <6p  >g  +  <6p  +  <6p  . 


fluctuation  theory  [jsee  Appendix]  ,  this  is  equivalent 

2  00 

to  replacing  the  value  of  <6p  >  =  kgCp  by  kgCp.  Doing 
this,  it  can  be  obtained  from  (104)  that 


<  6 


,  OO  . 

(y  -1)B 


T 


(113) 


■ 

.(  z£  orxa  art"  s  1o  .  v 

. 


:•  '  I)  rc  .  T  (: 
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This  value  is  then  a 


o 


Again  to  conserve  intensity. 


the  relaxation  mode  has  a  contribution  to  the  intensity 

1  BS  (v00— l)Bm  Brn 

equal  to  1  -  —  —  -  - - -  =  1  -  —  .  The  results 


V  oo 

Y  B, 


B, 


00 


Br 


00 


agree  with  those  of  (92)  with  G°°  =  0 . 
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Chapter  6 

VELOCITY  DISPERSION  FOR  A  SINGLE  RELAXATION 

PROCESS1* 

As  we  discussed  in  Chapter  1,  the  Brillouin 
components  in  the  spectrum  of  light  scattered  from 
liquids  can  be  interpreted  in  terms  of  a  Doppler  shift 
effect,  due  to  the  motion  of  the  two  sound  waves  moving 
in  opposite  directions.  The  shift  in  the  *  frequency ,  w, 
is  then  related  to  the  sound  velocity  v,  according  to 
(1-32)  and  (1-38),  by 

v  =  w/k  .  (1) 

However,  this  is  only  true  in  the  ideal  case  when  there 
is  no  absorption  of  sound  [see  (1-49)],  in  which  case 
the  hypersonic  waves  are  plane  waves  of  undiminished 
amplitude.  In  general  loss  mechanisms  are  always  present 
in  liquids  and  Eq.  (1)  will  have  to  be  modified.  It  seems 
that  this  has  not  been  generally  recognized  except  until 
only  very  recently  [30,26],  in  the  sense  that  the  positions 
of  the  experimentally  determined  Brillouin  peaks  are  taken 
to  give  the  sound  velocity.  This  has  not  generally  given 
rise  to  difficulties  because  sound  absorption  is  small 
in  most  cases  and  the  resolving  power  of  detection  devices 
always  gives  some  uncertainty  in  the  measured  values. 

i 

This  work  was  published  in  Bhatia  and  Tong,  Can.  J.Phys. 
17,  No.  3,  361  (1969). 


'*  ;  t‘  i  .  'I  t:  1>;.  •. 
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From  the  theory  and  calculations  which  we 
developed  in  the  last  chapter,  it  becomes  obvious 
that  Eq.  (1)  is  equivaleny  to  neglecting  the  g  terms 

p 

in  a(k,o>),  which  is  small  only  if  — ^  <<  1. 

In  terms  of  the  frequency  shift  of  the  Brillouin 
components,  which  we  denote  by  Wg ,  one  can  certainly 
define  a  ’sound  velocity’,  Vg,  by 


VB  ^B/k  * 


(2) 


The  quantities  v  and  Vg  are,  therefore,  not  the  same 
thing  in  general. 

In  a  recent  paper  by  Mountain  and  Litovitz  [32], 
negative  velocity  dispersion • in  Brillouin  scattering 
was  discussed.  By  using  a  perturbation  procedure  about 
k  ~  0,  they  stated  the  result  that  the  velocity  disper¬ 
sion  graph  (i.e.  v  versus  k)  for  a  single  relaxation 
process  will  show  a  negative  slope  for 


00 


V 


>  3 


(3) 


v 


In  this  chapter,  we  wish  to  discuss  the  fact  with  its 
consequences  that  v  and  Vg  satisfy  different  dispersion 
relations  so  that  when  one  talks  about  velocity  disper¬ 
sion,  one  has  to  mention  explicitly  which  quantity  one 
is  talking  about.  In  fact,  from  Eq.  (5-^8),  one  obtains. 


' 

.. 


nl  o  niupJ  Ixt-’  ni  no! 3*7  b  /Jiool9V  svl^Bsan 


c  ' 
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for  a  single  relaxation  process  (I.e.  putting  G°°  =  0), 
the  dispersion  relation 


v 


v 


00 

Bo“Bn 

i  +  & 


1  Bs  Bs-,  2 

-J  x 


T\ 


d 


(4) 


o 


which  yields  a  condition  for  negative  dispersion 


oo2 

>  5  ,  (5) 

V 

o 

and  this  condition  is  different "from 'that  in  (3). 

There  is  yet  another  possibility  of  defining  a 
sound  velocity.  As  is  well-known,  the  velocity  dis¬ 
persion  of  sound  waves  in  a  medium  depends  on  whether 
the  absorption  is  spatial  or  temporal  [33j34j,  corres¬ 
ponding  to  whether  one  makes  k- or  w  complex.  The 
definitions  of  the  velocities  v  and  - v'  for  the 
temporally  and  spatially,  absorbed  waves  are  thus, 
respectively , 


v  =  Re  co/k,  k  real,  co  complex  (6) 

and 

v'  =  u)[Re  k]-1  ,  k  complex,  03  real^  .  (7) 


The  difference  between  v  and  v’  is  particularly  signi¬ 
ficant  if  the  absorption  is  due  to  nonrelaxing  viscosities 
alone,  the  dispersion  being  negative  for  v  and  positive 
for  v’,  as  can  be  readily  verified.  The  spatial  absorption 


t 


Ultrasonic  measurements  usually  refer  to  v1. 


■  I.)  y  \t  j  )o  ;  9v  s»rl^  nwOirt-SIw  .1  eA 

.  j  n t  Vi  a1  1  -  o  i:oi8T  q 
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is  relevant  to  progressive  waves  whereas  the  temporal 
absorption  refers  to  the- decay  of  standing  waves  with 
time  in  a  medium.  The  latter  is  appropriate  to 
Brillouin  scattering  of  light' in  liquids,  and  Eq.(6) 
defines  the  same  velocity  as  that  in  Eq.  (4). 

We  shall  now  study  the  different  dispersion 
relations  which  arise  owing  to  the  different  ways  of 
defining  the  sound  velocity.  For  simplicity  we  shall 
consider  a  single  relaxation  process.  In  particular, 
we  shall  show  that  the  conditions  in  (3)  and  (5)  are 
exact  results. 


§6.1  Velocity  dispersion  for  v 


The  relevant  dispersion  equation  for  determining 
v  is,  from  (3-41)  or  (3-66), 


3  2  BS  2  iBs  2 

tuj  -  im  -  —  k  u)T  ,  +  -  k  =  0  . 

P  dp 

Ho  Mo 


(8) 


Eq .  (8)  can  also  be  obtained  from  (5-25).  We  recall 
that  the  imaginary  part  of  the  complex  roots  in  G(y)=0 
gives  the  velocity  of  the  temporally  absorbed  hyper¬ 
sonic  waves.  For  a  single  relaxation  process,  the 
dispersion  equation  (5-25)  reduces  to 

oo 

3  1  2  BS  1 

y  +  x  y  +  y  +  x  = 


0 


(9) 


.  O  3  i  US  9*1  $08X9 
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where  x=vQkTd.  Putting  y  =  ,  (9)  reduces  to  (8). 

We  assume  that  is  a  noninfinite ■ positive  number  so 

oo 

that  if  ■>  0,  -*  0  also.  Writing  now  a)=a)^  +  ia)2 

in  (8)  and  eliminating  a)2 ,  one  obtains  for  v  the 
equation 


( z-z  )  ( 4z-z  ) 2  + 


oo  ' 


X2  [8z2-6zzoo+^(  z2+i8zoo-27  )]  +z-l 


=  0  , 

(10) 


where  we  have  written 


z 


(ID 


Eq .  (10)  can  be  more  simply  obtained  if  one  defines 
y=y’  -  to  transform  (9)  to 


OO 


’3  ♦ 

DS  3x  27x- 


B 


OO 


3xB 


s  +  i]  =  0 

XJ 


(12) 


IT  ^ 

and  then  write  y’  =  y^  +  iy2,  eliminate  y^  and  define 
2 


^-5-  =  y 2  -  z  •  One  obtains  an  alternative  form  of  (10)  as 


B  -  z„  + 

3x 


3x 


27x- 


3x  x 


+  i]2 

Y-i 


(13) 


Hence , 


9  z 
9  x 


x^[4z2-3zZa)+i(z^+l8zoo-27  0+z-l 


{- 


3x- 


_ 1 _ 

(4z-zoo+Jt)(4z-3z00+  ^ 
3x  x 


} 


(14) 


-»V 


’ 


■ 
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a  2z 

a  x2 


-  IH  +  -TT  {2E')z‘;-3zooz+|(z^+18zoo-27)] 

3x 


1/  2 


+  I  If  +  8xz  5T  -  3^x  |f}/{(4z-z(B+-±j)(4z-3zc„+ij)} 

3x  x 


ft  0^-38.+^  (4  ^ 

X 


»-  •  i?Hl*  &  - 


'  {(  4z-zoo+^-)  (  4z  -  32^+  -i)}"1  .  (15) 

3x  x 

Prom  Eqs.  (10),  (13),  (1*0  and  (15),  we  can  draw  the 
following  conclusions: 


(1)  For  x  ->  0,  (10)  gives 


v 

v 


2  =  z  ~  1  -  If  x2(z«>  -  D(z„  -5)^1 
o 


a  z 
a  x 


-igx(z  -  l)(z  -  5)  -►  0  for  all  z 

'  OO  00  '  < 


11m  a  z 
x+0  -  2 

a  x 


=  -h(z  -  l)(z  -5)  . 
^  X  00  '  v  00 


Hence  lim  a  z  <. 
x+°  3  x2  5 


0  if 


> 

z  = 
< 


(11)  For  x  oo  ,  (9)  gives 


z  -+  z 


or  v  -►  v 


OO 


00 


a_z 
a  x 


o  2  3 

2z  x 


(zm  +  3)(z„  -  1)  0  . 


(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


)]> 


< D-rr+,  1  ~  A  ~  A  :+  A  - 


S  He  *10*1  0  (?  -  s)(I  -  s  )xjf-  «-  — 
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Since  from  general  thermodynamic  considerations, 
00 

Bs  ^  Bs  *  we  have 

Zoo  1 


(22) 


Hence  for  large  x 

If  >  0  •  (23) 


(iii)  The  condition  for  nontrivial  stationary  points 
3  z 

(i.e.—  =  0)  can  be  obtained  by  combining  (10)  and  ( 1 4 ) 

o  X 

to  give 


X4 - £-=-! - 

(Hz-z.j2(z-z  J 


(24) 


Hence  there  are  no  stationary  points  for  z  if  1  <  z  <  z  . 

9  z 

That  is,  for  this  range  of  values  of  z,  ^ —  is  never 

o  X 

negative . 

Furthermore,  from  Eq.  (10),  because  the  equation 
2 

is  quadratic  in  x  ,  one  can  see  that  for  every  value  of 
z  there  can  be  at  most  two  values  of  x,  since  x  >  0. 

One  may  may  conclude  then  that  negative  disper¬ 
sion  (decrease  of  z  with  increasing  x)  must  start  from 
the  origin,  and  this  occurs  when  z^  >  5 .  Figure  1 
sketches  the  variation  of  z  with  x  for  different  values 

of  z  . 

00 

For  x  <<  1,  one  has  the  expansion  given  by  Eq.(l6) 
which  is  identical  with  Eq.  (4). 


~  t  s  \  c 

.  rsa  r. 


'»  so ' i  >  .  d  t  (01 )  .  R  10 Tl  ? i  .rr 


t  s?  :  ::.6nl  a!.  9bul0:ioo  ^£i .  ysn  anO 

anj.efi*.  ion!  i  U:w  s  lo  a*  -  naob)  no  ’ 
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2  2 

Fig.  1.  Schematic  variation  of  z(=v  /v  )  with  x  on  -log 

log  plot  for  different  values  of  z  .  (a):  l  <z  <5  • 

CO  CO  3 

(b):  5<zoo<9;  (c):  z  >9.  For  z00=9J  the  minimum 

00  —i' 

in  curve  (b)  occurs  at  x=(27)  2  and  z  .  =0.  As 

mm 

shown  in  curve  (c)  for  z  >9,  there  is  a  range 
of  values  of  x  for  which  eq.  (10)  has  no  real 
positive  roots. 
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§6.2  Velocity  dispersion  for  vf 


This  velocity  is  defined  by  Eq.  (7)»  The 
relevant  equation  to  consider  is  also  Eq.  (8). 
Writing  k=k1-ikp,  one  can  eliminate- kp  to  obtain 

an  equation  in  kn  and  w.  This  equation  can  be 

1  2  r00 

v  ’ 

rewritten,  by  defining  z'  =  — p—  ,  zto  =  g—  ,  to  give 

v  °° 

o 


z,3x2  {(z  -1)2+  4z3x2}  +  4z  x2z,2{l  +  z  —  z^x2} 


00 


00 


00 


00  00 


+  4z’(l  -  2z2x2)  -4=0. 

'  /v\  ' 


(25) 


Prom  (25),  it  can  be  seen  easily  that  for  x  0,  z’  -*  1, 
or 

lim  v*  v  (26) 

x+0  ° 


For  small  x. 

Since  z  >  1 
00 

For  x  -*■  z 


,2 


z’  =  =  1  +  5-  x2(3zoo+  1)  (z„o-  1)  • 

V 


3  z 


3  x 


is  always  nonnegative  near  the 


v 


’2  B 


00 


00 


V 


B, 


=  z 


(27) 

origin . 


or 

v  ’  =  v 

00  00 

For  large  x 
,  3  z  1 

and  5—  is 

d  X 

since  (25) 


(28) 

>  z'  “  z»  -  — TT  (zco+3Xz„-l)  (29) 

4z  x 

00 

positive  for  all  values  of  z^  >  1.  Again 

2 

is  quadratic  in  x  ,  by  similar  argument  as 


9 


4 

- 


Ill 


we  had  for  z  in  (l-(iii)),  we  conclude  that  zf  is  a 
monotonically  increasing  function  of  x  for  all  x  and 
for  all  z^.  Figure  2  sketches  the  variation  of  zr 
with  x. 


§6.3  Velocity  dispersion  for  v 


B 


To  calculate  the  velocity  vg  defined  by  (2)  we 
need  the  expression  for  the  spectral  distribution.  This 
is  obtained  from  (5-24),  (5-25)  and  (5-31)  by  taking  the 
limits  a=0,  G°°=X2  =  0.  We  get 


00 


a( k,to)  =  2  Re 


lim 


2  BS 

y  x+y+(^  -  1)  x 


V  00 


V  k  .  „  -  Bn 

V=“ 7T  y3x+.y2+^-xy  +  1 


y  v  k 
o 


B, 


(  Z  -  1  )  X 

'  oo 


V  k  r  2  n  2 

° 


2,  2 
v  k 
o 


+ 


2  2 
(0  x 

2,  2 
v  k 
o 


[■ 


z 


(0 


2  -i 


00 


2,  2 
v  k  -1 
o 


(30) 


Hence 

9_a  ; 

9  (a) 


4(zoo  -  1)  tox 


v3k3 

o 


(1  - 


CO 


2,  2 
v  k 
o 


2  2 

)2+  w_2L  (Z  _ 
;  2,  2  1  °° 
Vok 


ca 


n  2 


2,  2 
v  k 
o 


) 


x  {x^(z  - 

'  oo 


(0 


2,  2 
v  k 
o 


)  -  2(1  - 


0) 


2,  2 

Vok 


)  - 


o  2  2 

2x  u) 

2,  2 
v  k 
o 


( z 
'  00 


03 


2,  2 
v  k 
o 


)}. 


(3D 


The  maxima  and  minima  in  the  a  versus  to  curve  are  given 
by 


fcijJ  -t;  >  >  su  t((i  -I  ,  ni  2  *iol  I  rt  9W 


'u) 

; 


■ 
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2  2 

Pig.  2.  Schematic  variation  of  zT(=vr  /v^)  with  x. 

zr  is  a  monotonically  increasing  function 
of  x  for  all  x  and  for  all  z  . 
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9  g(k,p)) 

9  03 


0 


(32) 


The  solutions  of  (32)  are 


03  =  0 


v2k2 


03 


2  =  (°)  l2z^  -  1  t  A**]  , 


3x' 


oh  p  2 

where  A  =  z  x  -  4z  x  +6x  +1 

00  00 

it  can  be  verified  easily  that 


By  evaluating 


(33) 


9  2a 

9  032 


9  O 


9  03 


l6(z  -l)x 


00 


(A)  “>03 


47T 

vok 


2  , 
03+  1 


-  2  2 
3x  o)± 

2,  2 
v  k 
o 


-  2z  x  +  1} 


00 


16(z<o_1)x  “+  i, 

+  in  A  2 

v  k 
o 


(34) 


The  discriminant  of  A  is  12(z  -l)(z  -3).  If  z  <3, 

00  00  00  3 

A  is  always  positive.  If  z^  >  3*  there  is  always  a 


range  of  x  for  which  A  is  negative  so  that  03+  and  03 
are  complex.  One  can  also  verify  that 


9  2q 

9  03 2 


03  -►0 


-(z2x2  -  2) 
x  00  ' 


(35) 


From  the  above,  we  can  say  the  following: 


(i) 


At  03=0,  a(k,oi)  has  a  maximum  or  a  minimum  according 


to  whether  z  x  -2  >  0  or  <  0. 


(ii)  Provided  03+  is  real,  nonzero  and  distinct  from  03  , 

2  -1 

a(k,o))  has  a  maximum  at  oj  .  For  x  >  (2z  )  ,  03, 

T  00  + 

2  —  1 

satisfies  these  conditions  if  A  >  0.  For  x  <  (2z  )  , 

N  oo '  ? 


I 


’  '  ..  K  '  ;•  '  f  - 


f<«i  J£  mumixBra  s  sari  (u>e>Oo 


. 1 1  •  ri<  d  ’  (  >  . .  \  i©  '  ,::i  1 ~ 
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2  2 

the  corresponding  requirement  is  that  A  >  (2zoox  -1)  , 

2  2 

which  is  equivalent  to  the  inequality  z^x  -2  <0,  i.e. 


a  minimum  at  u)=0  according  to  ( 35 ) .  For  z^  <  3,  a(k,oo) 
has  a  maximum  at  oj+(^0)  for  all  x.  For  z^  >  3,  however, 
there  is  a  range  of  values  of  x  for  which  go+  is  complex, 
(iii)  Provided  co_  is  real,  nonzero  and  distinct  from 
a)  ,  a(k,oj)  has  a  minimum  at  u)  .  From  (33),  it  can  be 
seen  that  these  conditions  on  go  are  satisfied  only  if 
(a)  A  >  0  and  (b)  x  is  greater  than  both  (2zoo)  and 


Taking  together  (i),  (ii)  and  (iii)  imply  that 
if  a(k,Go)  has  a  maximum  at  w,  (?*0),  then  it  has  a  minimum 
either  at  co=go  or  at  go=0;  in  the  former  case  it  has 
another  maximum  at  co=0. 

(iv)  For  x  0, 


2 


->  1 


(36) 


For  large  x 


00 


( z  -1  ) 

N  oo 


00 


(37) 


and 


>  0 


(38) 


o 


(w«>Oo  <£  >  s  .  (££)  ocf  -ni'DtoooB  0=u>  nwm.tn.tm  b 


i  30  £  ■  il  1  IK  rtO*:  X  l0  J9U1JE  '10  9^nBfI  Ja  at  9T3ri>1 

■  :  <  SO  ,1  '.)  0  .  ^  t£  T  T  in  B  2Bd  »J ,  >  )  0  t  .u> 

w  no  artolJ.  bnoo  ea9ti&  JerfJ  n 998 


U  it  ,,  ;mt  (iJ  )  Li  6  (it)  ,i  '  j  isritfos©*  snltfsT 

i  !  ■?.  :  C  s  [  S  8Bfi  ( U)  (  >0  O  li 

v.  "i /i  :  ni  ;  »co  tff*  >  u>®u)  cM  ii9ri;fx9 


■  ^  1  ■ 
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Also , 


0) 


2.  2 
v.  k 
o 


z  -3 


Z  X 
00 


(39) 


(v) 

ti 

+ 

3 

0)  if  A=0 . 
0 

If  z  < 
00 

3, 

this  is  never 

satisfied . 

For  z 

00 

>  3, 

u)+  and  u> 

are  equal 

to 

each  other  for 

two 

values 

of  X 

,  which  are 

roots  of 

the 

equation  A=0. 

At 

these  points  it  can  be  easily  verified  that 
8  to .  8  to 


?  2 

except  for  z  =3.  At  x  =  2/z  ,  to  =0,  to.^0  for  z  <4 
and  to+=0,  to__  is  complex  for  zto  >  4 .  If  zoo  =  4,  then 
(o+=(o  =0  at  this  point. 

The  foregoing  discussions  show  that  we  can  identify 

to+  with  Wg,  and  enable  us  to  make  schematic  plots  of 
2  2  2 

z  '=  to^/v  k  versus  x.  Such  plots  for  certain  values 

JD  u  O 

of  zto  lying  in  certain  ranges  are  given  in  Figure  3.  We 
see  that  if  zb  <3,  Zg  shows  negative  dispersion  which 
starts  at  the  origin.  Also  there  is  a  range  of  values 
of  x  for  which  there  are  no  Brillouin  peaks. 

The  difference  between  v  and  vg  is  particularly 
marked  if  zm  ^ 3 .  Unfortunately  no  liquid  is  known  at 
the  present  time  for  which  zm  >  3.  However,  even  for 


3  ■  ■  n  8 j.  a :  J  .  '  >  'I  .  '  V 


C\  3  riO'l  .i  c  lo;>9  of  '  f:  r'  e -ir 


.  £  -  s  *iol  ^q90xs 


r  o  -v  worts  snaiasuoaifo  aalogsio'l  9riT 

e 


>J  v\a0)  =  aS 


11 6 


2  2 

Fig.  3*  Schematic  variation  of  zD(=vD/v  )  with  x  for 

D  DO 

different  values  of  z  .  (a):  l<z  <3; 

00  CO  ' 

(b):  zm=3i  (c):3<2qo<4.  For  z ^  4  ,  the  branch 
showing  negative  dispersion  in  (c)  meets  the 
abscissa  at  x^/z^;  the  branch  showing  positive 
dispersion  in  (c)  remains  qualitatively  unaltered. 
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z  <  3>  the  difference  between  v  and  vD  (  and  v’)  can 
be  significant  in  the  region  x  ~  1.  Table  I  gives  a 
comparison  of  the  variation  of  v,  vg  and  v'  with  x 
for  =  2.24,  which  is  the  value  appropriate  to 
glycerol.  It  is  only  when  z  -  1  <<  1  or  sound  absorp- 

i 

tion  per  cycle  is  much  smaller  than  unity  that  v’-v-Vg. 


Table  I.  Values  of  Vg,  v  and  v’  for  several  values 
of  v0kTgyj  assuming-z^  =  2.24 


X 

<W2 

(v/v0)2 

(v’/v0)2 

0 

1 

1 

1 

0.2 

1.020 

1.037 

1.096 

0.5 

1.170 

1.327 

1.506 

1.0 

1.743 

1.928 

1.960 

1.5 

2.003 

2.098 

2.105 

2.5 

2.153 

2.188 

2.189 

OO 

2.24 

•^r 

OJ 

• 

C\J 

C\1 

• 

OJ 

-q'lOe  6  bnwoa  io  I  >>  X  -  s  nsriw  ylno  al  ctl 


eulBV  Xb*i9V93  lo'i  1  v  bnfi  v  K  :v  lo  89jjIbV 
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Chapter  7 

SPECTRAL  DISTRIBUTION  FOR  LIQUIDS 
HAVING  DISTRIBUTION  OF  RELAXATION  TIMES 

In  Chapter  5  we  calculatedthe  spectral  distribu¬ 
tion  of  light  scattered  by  a  relaxing  liquid  in  which 
both  the  bulk  and  shear  relaxation  processes  were  each 
characterized  by  a  single  relaxation  time.  There  are 
many  liquids,  however,  in  which  there  may  be  more  than 
one  relaxation  time  for  the  bulk  relaxation  -  for 
example,  for  thermal  relaxation  in  polyatomic  liquids. 

A  specially  interesting  case  is  that  of  liquids  of 
high  viscosity  in  which  both' the  shear  and  bulk  relaxa¬ 
tion  processes  have  to  be  characterized - by ■ a  distribution 
of  relaxation  times  QL3,l6j.  In  this  chapter  we  generalize 
the  treatment  of  Chapter  5  to  cover  these  cases. 

§7*1  Basic  Equations 

Following  the  same  procedure  as  in  Chapter  5,  the 
modified  linearized  hydrodynamic  equations  are 


(1) 


-VP  +  i  /  n’(t-t’)  grad  div  u(t')  dt ’ 
5  0 


(2) 


and  the  energy  transport  equation  is 


t  c kXAsn  %o  woi^uaiR  ?aia  omivah 


*  '  3  1  .  '  S  i.S  1  1  3.1(. B,’iO  1 


"•I  r  '  i-hjp±  J 


I  lanY^oib^rf  b^'  *i  ;-->rr 1  boJt’l  .  n 
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9_S 
9  t 


(3) 


where  we  assume  the  following: 

(i)  There  is  a  distribution  of  shear  relaxation  times 
g^(i)  such  that 

oo 

/  g_(T)  dT  =  1  (4) 

0  n 

nf(t)  =  G°°/  g  (t)  e~t//T  dx  .  (5) 

0  n 

(ii)  There  are  n  ordering  processes  each  being  charac¬ 
terized  by  an  ordering  parameter' corresponding  to  an 
ordering  force  Z^.  The  relaxation  equations  are 


9  £.  n 

•Tj T —  =  -  Z  L .  .  6Z  ,  i  =  1 ,  .  .  .  n  .  ( 6 ) 


Thus  (1),  (2),  (3)  and  (6)  are  to  be  solved 
together  to  give  p(k,t)  for  the  calculation  of  a(k,oo). 
For  independent  variables  it  is  convenient  to  choose 
S,  p  and  Z,  where  Z  denotes  Z, ,  Z0  ....  Z  .  This  will 
be  apparent  later  when  we  discussed  comparison  with 
experimental  data. 

Making  use  of  thermodynamic  relations  derived 
in  §3.2,  we  can  write 


2 

V^P 


Iq^^t  2  BS  2  n  2 

— - -  V  S  +  —  V  p  +  2  Pi  V  Z. 

Ctt  p  .  t  Z  .  1 

V  Mo  1=1  1 


2 

V  T 


T  0  T  3Bm  0  n 

—  V2S  +  — —  V2p  +  Z  T' 


'V 


P  C,r 
K0  V 


1-1  Zi 


2 

V  Z 


(7) 

(8) 


.  Tb 


36 


.stfj&b  IjeJnsml'ieqxs 
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where  P*  =  (^-£— ) 


t'  =  d^-) 


-j_  9  SjpjZ^  5  8  Z^  Sap,Z^  , 


and  Z.  denotes  Z-.  .  Zos  ...  Z.  ...  Z .  ,  ,  ,  ...  Z  . 
i  l5  29  i-l5  i+l5  n 

Eq.  (7-6)  can  be  written  as 


9£.j  ao  9^4  r, 

( — i)  (9_e.)  +  ( 1)  9_£  = 

'9  S  p,Z  9  t  9  p  '  S ,  Z  9  t 


9  €± 

?^9Z7^p  ,s3z! 

J  J  9  J 


x 


9  t 


+  L  .  ."]  6  Z  . 

i  iJ  i 


ij 


(9) 


It  can  be  shown  by  thermodynamics  that 


9^1 
( — -) 

l9S  Jp,Z 


=  -  T 


Z  . 

l 


9  q 

( — -) 

9  p  ;S,Z 


o 


Z. 

l 


(10) 


If  we  define 


-1  9 

*ij  ■  fzT  S ,  p ,  Z  ; 

J  J 


-1 

T  .  . 

IJ 


L  .  .  cb  .  . 
IJ  IJ 


(ID 


then  eq.  (7-9)  can  be  written  as 


7  l-l  +  K  K  |4  =  z  *7.  (.4-  +  4-)  6z4 


Z±  9  t 


2  Z  .  9  t 
P  i 
Ko 


ij  9  t 


ij 


(12) 


Eq.  (1)  and  Eq .  (2)  can  be  combined  to  give 


9Jl£  = 

9  t2 


4-  /n’ (t-t' )v2(|-2t)  dt ' 

3p0  0  3  t 


2 

VP  + 


(13) 


4 


*  W  ti“  -  tiT.  ;  :_s.q.e  tv 


SB  i  tl I  '  f  O  (C-')  . 


i 


y  to, 1.00  otf  mo  (S)  .p,l  brr  (.,)  ,pS! 


S+  c 
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§7-2  Calculation  of  a(k,w) 


As  before,  the  Laplace  transform  a(k,S)  has  to 
be  first  obtained.  This  is  done  by  multiplying  (3), 
(12)  and  (13)  by  V-1  eik.r-St  and  integrating,  using 

(5-15)  and  (5-16).  The  results  are  the  following  n+2 


equations 


T  6B„k2  p  T  S  T  k2 

p(k,S)  { — 2^ — }  +  3(5,3)  (  °.  0  +  -2 — } 


P  C,r 
Ko  V 


X 


'V 


2  1  ->  P  T 

+  E  K  T7  Z.(k,S)  =  -^-2-  S(k) 

j  zj  J  x 


(I'D 


p(k,S)  ’  {- 


SPv.  ,  (1+St. . ) 

— }+  s(5,s)  {s  t7}  -  z  -z — z.(5,s) 

J 


O 


_  d> .  .  t  .  . 

J  1J  1J 


Z  .  (k) 


Pz. 

=  p(k)  {■ — i-}+  S(tc)  {Tp  -  z  -4 — 

p„  1  J  ¥ij 


■,  1=1, ...n  (15) 


o 


p ( k , S )  {S2  +  -I—  +  G°°k2S  /  ^ 


g  ( t  )  t  di 


P  3  P 
*o  Ko 


1+Sx 


T,  3Bmk  Q 

+  s(5,s)  {  0  1 — }  +  i  k2  p'  z  (5,s) 

cv  j  h  J 


=  p(tc)  (S  +  ji-  G“k2  /  ^ 


g  ( T ) T  dr 


3p 


O 


1+St 


■}. 


(16) 


In  general,  Eqs.  (l4)-(l6)  can  be  solved  for 
the  n  +  2  unknowns  p(^,S),  S(1c,S)  and  Z^(tc,S).  For 
the  light  scattering  problem,  we  are  only  interested 
in  p(k,S),  as  only  <p ( -I* ) p (]£ , S ) >  contributes  to  the 


O 

{“~cT 


'  (  J  +  —  r°— >;  '  8 ,  )0  q 


V°oQ 


t  t  .  i5  > 


v° 


nap  (c  r)-(A.  )  ,  a  ,  hi"  .  >s  r.L 
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spectral  Intensity.  Also  from  fluctuation  theory 
[Appendix],  we  have  <p(-k)Z^(k)>  =  0  so  that  one 

-f 

can  neglect  all  the  Z^(k)  terms. 

Thus  we  can  write  our  solution  as 


p(k,S)  = 


A 


B 


where 


n+2 . 


A 


and 


p  T  . 
S(k) 


A12  p ( k ) 


^2  Pz  + 

P  1 
Ho 

T’  S(k) 

Z1 


B 


(17) 


are  2  determinants  of  order 


p  T  S  T  k‘ 
o  Q  I  Q 

X  cv 


VBTk< 

CV 


ST, 


'1 


2  ’  2  ? 
k  T„  ....  k  Tz 

L1  n 


k2P*  ....  k2p’ 
Z1  Zn 


(1+St11)  ( 1+St  ) 

In 


hlTll  Win 


4  Pj  p(k)+ 

pg  \ 

S(k) 


I 


S  T 


n 


(l+STni) 

I  •  •  •  • 

^nlTnl 


where  A^2 


[s  + 


4  G°°k2 

3  “ 


(t)t 


(1+St  ) 

_ _ nn 

6  T 

Ynn  nn 


1+St 


(18) 


. yJtan9 Jnl  iBiJosqa 

..  0  r  U. 
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B 


W‘ 

2 

P  cTr 
Ko  V 


B 


12 


S_  p- 
2  rZ, 
Po  1 


S_  p’ 

2  rZ 
po  n 


2  2  ’  2  » 
p  T  S  T  k  k  T„  .  k  T7 

-2^-  +  -2—  Zl  Zn 

X  °v 

rr®6BTk  k  PZ1  .  kZ?'w 

S  t’  -(1+St11)  -( 1+St11 ) 

^ll7!!  ^llTln 


S  T, 


n 


(1+Sti;l) 


4> 


nlTnl 


- ( 1+St  ) 
_ _ nn 

<b  t 
Ynn  nn 


(19) 


,2  ,  BSk  ,4  G”k2S  ,  sn(T)T  dT 
where  B12  S  +  Pq  +  3  Pq  1  1+St 


Anticipating  that  p(k,S)  is  to  be  considered  only 

$B 

in  <  p(-k)p(k,S)> ,  and  that  <  p(-k)S(k)>  =  -  — 2^-<p(k)p(-k)> 


o 


[see  Appendix],  we  can  replace  S(k)  in  (17)  by  - 


3B 


T 


p (k)  . 


We  now  try  to  simplify 


A 


and 


B 


by  applying 


some  elementary  operations  on  them. 


We  multiply  the  1st  column  of 


A 


2.  2 
P  k 

by  ■■■°  , ’  and 
p(k) 


the  2nd  column  by  k  /S.  This  will  increase  the  determinant 

2  4  ->■ 

by  the  factor  pQk  /Sp(k).  Then  we  add  to  the  1st  column 


.r  ^il  no  an  i  Jsrraqo  *.•  chv  oJ.  omoe 
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3B,p  times  the  2nd  column.  This  will  not  change  the  value 
of  the  determinant.  Finally  we  interchange  the  1st  and 
the  2nd  columns.  This  will  change  the  sign  of  the 
determinant . 


For  the  determinant 


B 


,  we  shall  multiply  the 


2  2  2 
1st  column  by  pQk  /S  and  the  2nd  column  by  k  /S,  and 


then  interchange  them. 

The  results  of  such  operations  yield 


A 


Sp(k) 

2  .  4 


,  2  p  T  S 
k  ^ 


S  ^ 


o  o 

-r  ~ 

X 

2 

T  K 

+  "3— } 

CV 


W 


Sc 


V 


2  » 

k  TZ 


2  » 

k  TZ 


n 


W 


Sc 


V 


Pok2[Ai2+ 


2  2  2 
To6^BTk^ 

SpoCV 


■] 


2  » 

k  PZ 


1 


k2p’ 


n 


k'Tz 


k2pZ 


■(1+ST11> 

<I)11T11 


- ( 1+STni ^ 

^nlTnl 


k2T’ 


n 


k2p’ 


n 


-(1+STln) 

^lnTln 


•  •  •  • 


- ( 1+St  ) 
_ _ nn 

(b  t 
Ynn  nn 


(20) 


)1  - ;  snoJ  Je  *qo  riou  ?.  lo  a^Iu  sen  ori'. 


'  r- 


5 
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B 


274 

k 


i  2  p  T  S 
k  /  Ko  o 

S  1  X 
T  k2 

+  ^— } 

CV 


W 


Sc 


V 


k2T' 


k'Tz 


n 


W 


Sc 


V 


2  2 

p^k 

o 


B 


12 


k2pZ 


k2p' 


n 


k2T' 


k2pz 


■  (1+Sti:l) 
*11T11 


-(1+Slnl> 

<t’nlTnl 


k2pz 


n 


k2p' 


n 


-d+STln) 

^lnTln 


■  ( 1+St  ) 
_ _ nn 

(p  T 

nn  nn 


(21) 


We  recall  that  d>..=d>...  Further  it  is  reasonable  to 

ij  Ji 

assume  that  L..=L..,  then  (20)  and  (21)  are  symmetrical, 

i  J  J  * 

which  simplifies  the  expressions  considerably. 

Consider 

6^i  "hfS^PjZ  6S  +  (9"7“)S,Z  6p+  ?(aT7)p,S,z'  6Zj  5 

J  J  J 


i=l 


n  (22) 


This  can  be  written  as 


[6£]  =  ss 


IA 
9  S 


+  6p 


9_£ 

9  P 


3«i 

9  Z  . 


[6Z] 


(23) 


<nT8+I)- 


■q  necfi.  r  ■  d  hbo  airlT 
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where  (jS^J  , 


i.e.  [$C\  = 


p_L] 

[Ls.1 

3  S. 

> 

L9  pJ 

and  ffiZ]  are  column  vectors. 


f6^l\ 

<5^ 


\66  / 


etc. ,  and 


3_ii 

3  Z  . 
j 


is  an  n  x  n 


symmetric  matrix.  Hence  there  exists  a  non  singular 
similarity  matrix  d  such  that  D 


3  Z  . 
J 


n-l  . 

D  is 


diagonal.  If  we  now  define 


PCI  =  d[6S] 


[SZ]  =  D[6Z]  , 

then  we  have  the  result  that 


(24) 


I 

r3hi 

roJ 

C_J. 

3  Z. 

L 

6  .  . 
ij 


(25) 


So,  without  loss  of  generality,  one  can  assume 
that  depends  only  on  p,  S  and  the  corresponding  6Z^ 


t 


Or, 


3_q 

3  t 


-hi  6Zi 


(26) 


This  kind  of  relaxation  where  each  parameter  can 
be  assumed  to  relax  independently  of  the  others  is  usually 
called  relaxation  in  parallels.  [13,  l6j 


,10 


12? 


with  the  understanding  that  any  t  defined  from  (2 


is  at  constant  S  and  p,  i.e.  x^y. 

We  shall  now  write  <j>^  =  0^,  x^  e  and 

again  apply  elementary  operations  on  A  and 


6) 


B 


to  evaluate  these  determinants,  remembering  that  now 


( 1+St .  . ) 

0  .  .  T  .  . 
1J 


is  diagonal.  Let  and  denote  the 


ith  column  and  ith  row  respectively.  We  now  apply  the 
following : 

^iTi  2  » 

(i)  Cp  +  £  i+St .  k  .  c2+iJ  followed  by  similar 
operations  on  Rp . 

0.x.  p 

(ii)  ci  +  E  k  _  c2  +  iJ  followed  by  similar 


operations  on  R^ 


A 


The  determinants 


A 


and 


B 


then  become 


(1+Sx. ) 


+  - -  .  -  *  n  „  '  "  i 


-r%—  pOc)  p(s)  (-d“  &  -t— — 

p5k°T  i  i 

Ko  o 


B 


(1+Sx. ) 


— £T% —  G(s)(-I)n  n  ,  ■ 
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Hence  Eq.  (17)  becomes 
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F(S)  =  S2  {1+  Xk 
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We  next  write  down  some  thermodynamic • relations , 
for  the  case  of  multiple  relaxation,  which  are  generali¬ 
zations  of  eqns.  (17),  (24)  and  (27)  of  Chapter  3.  One 
may  verify  that 
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00 
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If  we  now  go  to  the  continuum  limit,  so  that 
there  are  many  relaxation  times  in  a  small  range 
between  and  +  dx^,  we  can  define  three  distribu¬ 
tion  functions  as  follows: 
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s. 
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which  are  each  normalized  to  unity,  i.e. 
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We  can  then  replace  in  (29)  and  (30) 


z  Po1  h  pz2  f(Ti}  -  (Bs-Bs}  f  gd(T)f  (T)dT 
1  1 
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where  f(x)  is  any  function  of  x.  It  may  be  verified 
from  (3*0  that  only  two  of  the  three  distribution 
functions  are  independent,  since 


g^(x)  a  /gd(T)g’(x)  .  (36) 

With  the  substitution  of  (35)  and  (36)  into  (29)  and 
(30),  one  sees  that  the  expression  for  the  spectral 
distribution  will  in  general  contain  two  distribution 
functions  for  bulk  relaxation. 

We  next  observe  that  if 


TZ  =  constant  (independent  of  i) 

i  1 

=  A,  say  ,  (37) 

then  it  follows  from  (35)  that 

gd(i)  =  g^U)  =  6^(t)  *  (38) 

so  that  only  one  distribution  function  is  sufficient  for 
the  description  of  bulk  relaxation  processes.  When  (37) 
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Is  satisfied  we  shall  say  that  all  the  individual  bulk 
relaxation  processes  in  the  fluid  are  of  the  same  type. 

It  is  interesting  to  note  that  condition  (37)  is 
satisfied  if  all  the  individual  relaxation  processes 
are  either  (1)  density  induced  or  (2)  pressure  induced 
or  (3)  temperature  induced.  The  definition  of  these 
three  types  of  relaxation  processes  was  given  in  Chapter 
3,  and  one  may  verify  that  for  these  cases  A=0,  Top~13Cp1, 
respectively. 

In  principle,  there  is,  of  course,  no  reason  why 
in  a  given  liquid  all  the  relaxation  processes  be  of  the 
same  type.  However  it  is  known  that  in  highly  associa¬ 
ted  liquids  like  water  and  glycerine,  the  relaxation 
processes  are  predominantly  pressure  induced,  whereas  in 
other  liquids  (Kneser  type)  they  are  primarily  temperature 
induced  [13,  16] .  Since  use  of  (37)  simplifies  the 
analysis  considerably,  we  shall  calculate  here  spectral 
distribution  for  this  case  only. 

§7*3  Approximate  Determination  of  Spectral 

Density 

In  the  preceding  section  we  obtained  an  expression 
for  a(k,S)  =  F(S)/G(S)  for  the  case  where  the  relaxing 
bulk  viscosity  is  characterized  by  a  discrete  set  of 
relaxation  times.  Also  we  gave  a  prescription  for  going 
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to  the  limit  where  these  relaxation  times  lie  in  a 
continuous  range.  When  the  relaxation- times  form  a 
discrete  set,  they  are  usually  only  few  in  number. 
In  Chapter  5  we  have  already  discussed  the  case  of 
two  relaxation  times,  and  for  this' reason  we  shall 
not  discuss  the  discrete  case  further- here " . 


For  the  continuum  case ,  expressions' ( 29 )  and 
(30)  can  be  written  as,  on  using  eqns.  (35)  and  (38), 
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We  shall  now  proceed  to  factorize  G(S)  approximately. 
It  Is  convenient  to  define  the  following  dimensionless 
quantities : 
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Hence  Eq.  (7-27)  becomes 
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The  low  and  high  frequency  limits  are  then  defined 


OO 


by  v  k  t  -*■  0  and  v  k  t  -*  00  respectively,  where  t  =  /  xg(x)dT 
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is  the  mean  relaxation  time.  This  definition  of  the 
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oo 


has  the  value  Mq  and  M  respectively  If  one  takes  the 
low  and  high  frequency  limits. 

f 

We  shall  in  the  following  analysis  assume  that 

ax,  and  ax  «  1,  where  x ,  =  v  k  t , ;  x  =vki 
d  rj  9  dod^nor) 

Writing  (44)  as 
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we  see  that  the  small  root  is 
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Therefore,  G(y)  =  (y  +  yQ)j ~j2+  (c2+£2)y  +  c1  +  £^,(48) 
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Expanding  (48),  and  comparing  coefficients,  we 


get 


a 

"yo  yA 

00  ’ 


B  0  — B  o 

a  •  r  S  S  ,  x 
yA  1  Bq  J 


)  00 

l+yx  gd(x)dx  +  3  h  s  Try?  V*)dx} 

(49) 


Define  now 


26  =  c2  +  e2 


d  =  ci  +  ei 


(50) 


Therefore,  G(y)  =  (y  +  yQ)(y  +  26y  +  d) 


5 


136 


oo 

®5~®S  x 
26  =  ■ /  x 


gd(x)dx  +  £  §-  /  1+j^  Sr,(x)dx 


00 


B, 


l+yx 


+  —  (i  -  — )  -  9lLa — V.  r  yx —  g  (v)dx 

A  V  Y;  A  J  l+yx  )a 


(52) 


OO 


J  _  1  ,  Oi  ^  G  /-.  1  \  /• 

d  -  1  +  A  3  (1  -  Y>  f 


I  00 

irfr  «n<x)dx  +  ^ir  “  i  k  * 


f  g  ( x )dx  f  -rz*--  g,(x)dx  +  ^  “  /  ttttt  g^(x)dx 


l+yx  &n 


l+yx  °d 


A  y  l+yx  °d 


B  -B  B  B 

+  “L— r - g — -I  f 


X 


l+yx  °d 


gj  (x)djc  . 


(53) 


We  now  write 
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Comparing  (54)  with  (43),  we  get 
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Now  consider  Eq.  (54).  The  first  term  corresponds 
to  a  central  component  arising  from  thermal  conductivity. 
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just  as  in  the  case  of  single  relaxation  times .  The 
second  term  contains  other  central  components  and  the 
Brillouin  components.  In  general  it  is  very  difficult 
to  separate  out  these  components.  We  shall  use  an 
approximate  method,  taking  as  our  guideline  that  the 
results  should  reduce  to  those  that  we  previously 
deduced  for  the  case  of  single  relaxation  times.  We 
first  assume  that  the  sound  absorption  is  small 
( T/vk  <<  1)  so  that  the  Brillouin  components' and 
central  components  are  well  separated.  Now  from  (5*0, 
consider 


y2  +  26y  +  d 


(56) 


where  we  recall  that  6  and  d  are  functions  of  y.  Relevant 

to  our  problem,  y  is  eventually  to  be  considered  complex 

and  is  equal  to  — ^  .  Hence  6  and  d  are  also  complex. 
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For  the  description  of  the  Brillouin  components,  we  require 
a  form  of 
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which  is  to  be  extracted  from  (7-55). 
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where  26  =  26^  +  262 


d  =  d^  +  d2 
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(59) 


d^,  d2,  6-,  being  real,  and  d2,  beinS  pure  imaginary. 
We  shall  now  replace  some  of  the  y’s  in  (58) 
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where  and  g  are  still  to  be  determined .  We  shall  find 
them  by  comparing  (64)  with  the  second  term  of  (54). 


Writing 


b  =  b^  +  b2 

c  =  c-^  +  c2 

c2  =  y  c2  ,  (65) 

where  b^,  c-^,  c2  are  real  and  b2,  c2  are  pure  imaginary, 
we  identify 


2aB  Cbi  +  c2^y'=v/vo 

2a  0  +  2g  f-  =  [c  +  b2  y]  ,  /  .  (66) 
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Therefore , 

2aR  =  -  {1 
B  y 


VBs 

Bn 


2  i  °° 

/  ■  -X-  ■  ■  K  (x)dX  _  1  5_  / 
J  .  ,  2  sdl  ;  3  B„  J 

1  +  zx  S 


x 


l+zx‘ 


2  gn(x)dx} 


(67) 
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2g  - 


00 

Bs-Bs 


4  G 


_£  {26+2e[-H3^  /  gd(x)dx  +  f  57  ' 


X 


1  +  ZX 


l  +  zx‘ 


X  g^(x)dx] } 


(68) 


Hence  the  spectral  distribution  is  given  by 

2f  2f 


2a  S 

a(k,aj)  =  2°-  +  a 


a)  +S 


t 


B  l(oo+vk)2  +  r2 


(aj-vk)  2  +  f2- 


+  2g 


oo+vk _ 

(oo+vk)  2  +  r2 


oa-vk 

(ai-vk)  2  +  r2 


+  central  components 


due  to  relaxation  processes. 


(69) 


In  the  above,  we  have  written 


S  =  v  k  y  =  , 

o  o  Jo  p  c p  3 

O  r 


o 


The  intensity  ratio  is,  as 
1 


2a 


-  1. 


B 


1  -  i  5  r  =  vok  0.  (70) 

before,  given  by  R=  1^/21^= 


The  central  components  arising  from  the  relaxation 
processes  are  in  general  quite  mixed  together  and  the 
resultant  distribution  is  in  general  not  Lorentzian. 
However,  by  virtue  of  conservation  of  intensity,  their 


'> 


combined  contribution  is  equal  to 


■  1  -  a0  -  2°B 


(71) 


One  can  obtain  a  very  approximate  expression  for 
such  a  resultant  distribution.  Assuming  that'  it  is 

approximately  Lorentzian,  one  can  write  down  its  con- 

o  * 

tribution  to  a(k,S)  as  •gpgT  •  •  Neglecting  the  contribution 
to  the  intensity  at  oj=0  from  Brillouin  components  we  have 


—  =  lim 


by  +  c  1 


y+0  2,OJf  , ,  v  k 
J  y  +26y+d  o 


00 

1  ■  rBS“BS  -  4  G°°  -  1 

Y  {  B„  Td  +  3  B„  Tn}- 


(72) 


Hence 

00 

qf  .  /  rBS“BS  -  ,  4  G°°  - 

s  y/  {  Td  +  3  B0  Tn}  • 


(73) 


The  contribution  to  a(k,oj)  is  then  given  by 


2a  f  S  1 


2  2 

a)  +Sf  ' 


(74) 


It  can  be  verified  easily  that  in  both  low  and  high 
frequency  limits,  (6l)  and  (67)  both  reduce  in  the  case  of 
single  relaxation  times  to  the  expressions  for  2f  and  2 aB 
given  in  §5.5.  The  expression  for  v  /v  in  (5-47)  and 
(5-48)  can  also  be  obtained  as  such  limits  of  (63). 
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Chapter  8 


DISCUSSION  AND  SUMMARY 

In  the  preceding  chapters  of  this  thesis  we  have 
given  a  formulation  of  the  problem  of  calculating  the 
spectral  distribution  of  a  liquid  which  possesses  both 
a  relaxing  bulk  and  a  relaxing  shear  viscosity.  The 
former  is  ascribed  to  the  relaxation  of  ordering  para¬ 
meters  and  the  thermodynamic  theory  of  relaxation 
processes  has  been  applied  to  the  calculation  of  the 
density  correlation  function,  and  hence  of  the  spectral 
intensity.  It  was  shown  that  such  a  treatment  enables 
one  to  determine  the  extent  to  which  the  integral  repre¬ 
sentation  for  the  relaxing  bulk  viscosity  is  appropriate. 

As  for  the  shear  viscosity,  it  must  be  taken  as 
the  low  frequency  limit  of  a  shear  relaxation  process  if 
the  occurrence  of  Brillouin  components  is  to  be  explained 
in  liquids  of  high  viscosity.  In  general,  the  use  of  a 
relaxing  shear  viscosity  is  necessary  whenever  mri/B^  £  1. 
This  condition  is  easily  met  in  liquids  of  high  viscosity 
for  the  values  of  w  involved  in  light  scattering  experi¬ 
ments.  Parenthetically  it  may  be  mentioned  that  in 
neutron  scattering  this  condition  is  met  in  other  liquids 
also,  and  it  is  well  known  that  the  neutron  scattering 


. 


r 
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data  requires  that  high  frequency  shear  waves  propagate 
in  liquids  without  much  absorption,  indicating  the 
necessity  of  using  a  relaxing  shear  viscosity. 

In  Chapter  5  we  have  calculated  the  spectral 
intensity  for  the  case  where  both  the  bulk  and  shear 
relaxations  are  each  characterized  by  a  single  relaxa¬ 
tion  time  for  many  liquids  (e.g.,  the  so  called  normal 
liquids).  However,  one  usually  has  x^^v^t^)  <<  1, 

whereas  x^  can  have  values  such  that  x^  £  1.  As  already 
mentioned  in  Chapter  5,  this  special  case  of  our  results 
has  been  discussed  by  previous  workers  in  some  detail 
[24  ,  27,  323.  We  shall  therefore  confine  our  comments 
here  mainly  to  liquids  of ■ high' viscosity . 

The  viscous  liquids  have  the  property  that  their 
shear  and  bulk  viscosities,  and' hence *t  and  t , ,  are  of 
the  same  order  of  magnitude  at  -  a  given  temperature,  and 
that  both  increase  sharply  as  the  temperature  decreases . 

It  is,  therefore,  possible  by  varying  the  temperature  to 
obtain  the  spectral  intensity  of  scattered  Tight  for 
different  values  of  x^  and  x^  covering  the  range 

x2,xd  <<  ^  x2,xd  >>  ax2,axd  <<  1  to  ax2,axd  >>  1. 

The  expected  qualitative  behaviour  of  the' spectral  dis- 

I 

tribution  as  one  goes  from  low  frequency  (high  temperature) 
to  high  frequency  (low  temperature) limit  is  as  follows: 


« 
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(1)  In  the  low  frequency  limit,  x^  <<  1,  xp  <<  1, 
the  central  components  due  to  the  relaxing  modes  have 
very  large  half-widths  but  small  total  intensity,  so 
that  they  form  a  weak  background.  The  central  line  is 
predominantly  due  to  thermal  conductivity,  and  is  very 
sharp.  The  Brillouin  components  are  also  very  sharp 
and  somewhat  taller  than  the  central  components. 

(ii)  As  temperature  decreases  so  that  xp,xd  ~  1* 

0(=r/vok)  is  not  very  small  compared  with' unity.  The 
Brillouin  components  then  become  broader  and  shorter. 
Further  because  the  g  terms  (see  5-61) 'have  now  become 
significant,  the  shape  of  the  Brillouin  components  will 
deviate  from  the  Lorentzian  shape,  with  the' maxima 
displaced  towards  the  centre  (u)=0)  of  the  spectrum. 

The  central  components  due  to  the  relaxing  modes  now 
have  smaller  half  widths  and  greater  heights. 

(iii)  As  temperature  decreases  further  so  thatxp,x^  >>  1, 
the  Brillouin  components  become- very  short  compared  with 
the  central  line  which  is  now  predominantly  due  to  con¬ 
tributions  from  the  relaxing  modes. 

The  above  behaviour  is  in  qualitative  agreement 
with  the  measurements  of  Knaap ,  Gornall  and  Stoichef'f 
[36J  and  Pinnow,  Candau  et  al.  [35"]  on  glycerine. 


. 
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For  a  quantitative  comparison  one  has  to  take 
into  account  the  fact,  which  is  well  known  from  ultra¬ 
sonic  work,  namely,  that  both  the  bulk  and  shear 
relaxations  have  each  to  be  characterized  by  a  multi¬ 
plicity  of  closely  spaced  relaxation  times.  In 
ultrasonic  work  the  only  relevant  quantity  for  bulk 
relaxation  is  the  complex  modulus  Bg ( ito )  ,  so  that  it 
is  sufficient  to  interpret  ultrasonic  data  in  terms 
of  one  distribution  function  for  bulk  relaxation  and 
one  for  shear.  These  are  the  distribution  functions 

g,  (t)  and  g  (t)  introduced  in  Chapter  7«  In  contrast 
—  n 

to  this  we  saw  in  Chapter  7  that" the  -  expression' for 
spectral  intensity,  in  general,  contains' two' distribution 
functions  g^(x)  andg^(x);  only  when 'all 'the  bulk  relaxa 
processes  are  of  the  same  type  does' the  expression  for 
the  spectral  intensity  contain  one  distribution  function 
(gd(T) ) . 

Montrose,  Solovyev  and  Litovitz'  [30]  have  empi¬ 
rically  generalized  Mountain’’ S  ' expression"  for  the  case 

of  single  relaxation  time  to  cover  viscous  liquids. 

+ 

Their  expression  thus  contains  one  distribution  function 


Like  Mountain’s  expression  in  [24],  it  is  thus 
strictly  valid  if  all  bulk  relaxation  processes  in  the 
liquid  are  density  induced. 
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for  bulk  relaxation  and  one  for  shear  relaxation, 

Pinnow,  Candau,  La  Macchia  and  Litovitz  [35]  have 
used  this  expression  to  interpret  their  light  scatter¬ 
ing  data  in  glycerine.  They  findthat  g^(i)  which 
fits  the  light  scattering  data  is  substantially 
different  from  that  required  to  fit  the  ultrasonic 
data,  whereas  as  for  shear  relaxation,  the  two  are 
practically  the  same. 

Our  work  suggests  two  possible  explanations 
for  the  discrepancy  ref erred' to ' above : 

(1)  all  the  relaxation  processes  are  of  the  same 
type  but  they  are  not  density  induced; 

(2)  all  the  processes  are  not  of  the  same  type. 

This  is  supported  by  the  fact  that  it  is  experimentally 

00 

found  for  glycerine  that  $  so  that  all  the  bulk 

relaxation  processes,  taken  as  a  whole,  are  not  purely 
density  induced  [see,  for  example.  Ref.  (13),  p.  254]  . 

The  above  discussion  and  summary  refers  to  the 
work  described  in  Chapters  5  and  7*  In  Chapter  6,  we 
drew  attention  that  one  should,  in  general,  distinguish 
between  the  velocity  v  of  the  temporally  absorbed  sound 
waves  and  the  velocity  vg  measured  by  the  relation 
vR  =  wR/k,  where  coR  is  the  shift  in  the  frequency  of  the 
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Brillouin  components.  The  difference  between  Vg  and  v 
is  significant  whenever  the  sound  absorption  per  cycle 
is  not  small  compared  with  unity.  A  numerical  example 
shows  that  the  difference  could  be  as  much  as  10  to  1 5l 
in  relaxing  liquids. 
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APPENDIX 


FLUCTUATIONS  OF  THERMODYNAMIC  QUANTITIES 


A.l  Method  of  Calculation 


Let  a  physical  quantity  x  deviate  from  its  mean 
value  x  by  a  quantity  6x.  If  the  probability  distri¬ 
bution  for  such  a  fluctuation  is  of  the  form 


a)  ( 6x ) 


cc 


exp  {- 


b  ( 6x ) 


2  k 


B 


(1) 


where  b  and  kg  are  constants  (for  later  purpose,  kg  will 
be  Boltzmann  constant),  then  the  distribution  is  Gaussian, 
and  the  mean  square  fluctuation  is 


<  6x2>  =  kg/b  ,  (2) 

If  there  are  simultaneous  deviations  of  several 
physical  quantities  x^,  ,  xn,  and  if  the  probabi¬ 

lity  distribution  is 


a)(6x1  ...6xn)  a  exp  {-  Z  Bik(Sxi<Sxk^  * 

B  i  ,k 

then  it  can  be  shown  that  [29] 


(3) 


<  6xi6xk>  =  kg  Blk  ,  (4) 

where  B-^  is  the  inverse  matrix  of  B  =  | | B 

Let  us  now  consider  a  volume  element  V  of  fluid 

o 

embedded  in  a  large  fluid  system.  We  assume  that  VQ  is 
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large  enough  to  allow  thermodynamic  description.  Although 
the  fluid  in  VQ  is  at  equilibrium,  it  will  at  any  instant 

i 

of  time  experience  fluctuations  of  the  various  thermo¬ 
dynamic  quantities.  The  probability  of  a  fluctuation  is 
related  to  the  minimum  work  required  to  carry  out  rever¬ 
sibly  such  a  fluctuation  by  pSQ 


OJ 


-Rmi-  ^/k^T 
a  e  mm'  B  o 


(5) 


Rmin 


6U  -  T6S  +  P6V  -  Z6 £  , 


(6) 


where  6U,  6S  and  are  the  changes  in  energy,  entropy 
and  the  ordering  parameter ■ of ' the  fluid  in  V  .  Writing 
U=U(S,V,£)  and  assuming  that  the  fluctuations  are  small, 
we  have  from  ( 6 ) 


R  .  -hr 
min  L 


2  2  2  2 
9  U  .  a  U  „r2  .  3  U  „2,  „  9  U 


6S  + 


9  S‘ 


9  y2  6V  +  9  ^2  +  29  S9  V 


6S  6V 


+  2  9-tff  68  6^+  29^T  6V 


h\jS  <S(as)v^  +  6V  <S(av)s  ^ 


+  6(9 


=  h  (6S  6 T  -  6P  6V  +  SZ  ) 


(7) 


If  we  substitute  (7)  into  (6),  choose  three  independent 
variables  and  expand  the  rest  in  terms  of  them,  the 
probability  distribution  becomes  Gaussian  and  the  results 
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of  (3)  and  (4)  can  be  used  to  calculate  the  mean  square 
fluctuations  of  the  various  quantities. 

A . 2  Mean  Square  Fluctuations  of  the  various  thermodynamic 

Quantities 


Eq.  (7)  can  be  written  as 


2  R  .  =  6S  6T  +  -  6P  6p  +  6Z 

mm  p  K 


(8) 


Taking  p,T  and  Z  as  Independent  variables  and  making  use 
of  ( 3-14 ) ,  we  have 


CV  2  ^o^T  2 
2  R  .  *  =44  6T^  + 

min  T  2 

0 


6p‘ +  <B>t,p6z2+'  2<H>T,pfiT4z 


In  (9)j  6p  is  independent  of  the  others 
( 2 ) ,  we  have 

,  ,.2.  _  pokBTo 
<6P  >  " 


(9) 

Hence,  using 


(10) 


<  6p6T>=  0 


(11) 


<  6p6Z>=  0 


(12) 


The  other  fluctuation  quantities  in  (9)  can  be  evaluated 
by  defining 


o  _  D  _  _ 1  ,3_£n 

11  m2  s  d22  T  ^9  Z  T,p 
T  O 
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Hence 


B|  =  -4r  { 


i  (i±)  . 

t2  Tq  ^  ZJT,p 
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(9_S)2  } 

V9  Z;T,pJ 


T  (1-S-) m  ~  c7t  >  using  (3-22) 
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using  (3-34) 
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<  6 Z  6T>  =  - 


kB  B12 
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k;  T 

kB  io  ,9  Ss 

co  d  6  '  T ,  p 

CV 


(15) 


The  mean  square  fluctuation  of  number  density  can  be 
obtained  from  (10)  as 


.  ,  2.  no  kB  To 

"  6n  "  =  V  Bt 
o  T 


(16) 


The  mean  square  fluctuation  of  the  number  of  particles 


is 


<  6N2> 


n  kD  T  V 
o  B  o  o 

Bm 


(17) 


If  we  now  take  S,  P  and  Z  as  independent  variables. 


we  have 


2  R 


mm 


h  ss2+  fe  6p2+  <H>s.p  6Z"+  5P5Z 


(18) 
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Prom  (18),  we  immediately  have 


<  6S  >  =  kg  Cp 


(19) 


The  determinant  required  for  the  evaluation  of  the 
fluctuation  in  other  variables  is 
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<  6S  6P>  =  0 


(22) 


<  SS  6Z>  =  0  . 


(23) 


If  we  now  take  S,  P  and  £  as  independent  variables. 


we  have 


. 

A-6 


2  Rmin  “  ^  632  +  ^  6p2+  <K>S,P  6«2+  2#S,P 

°P  S  '  '  (24) 


The  value  of  the  determinant  for  S  and  £  is 


T 

1  •  r  O  /3  Z  x 

rn2  ^  oo  l3 

io  CP 


(U)  2  J 
v9  £  '  S  ,  P 


1  (14) 
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,  using  (3-24). 


Hence 
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^  > - = -  (§T)s,P 

cp 


=  kB  Tq  (|^)T  p  ,  using  (3-33) 


<  6S6£>  = 


,  r9_T\ 

B  CPl>9  Z;S,P 


<  6P6£>  =  0  . 


(25) 

(26) 
(27) 


For  p}  T  and  £  to  be  independent  variables,  we  have 

OO 
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2  Rmin  2  BT  6p  +  T  6T  +  (3£)TjP<$S  +  p  ^^T,p6p6^ 

P°  °  °  (28) 


The  value  of  the  determinant  for  p  and  £  is 
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,  using  (3-16) 
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Hence 


<  6p6£>  = 


p  kD  T  _  ^ 

Ho  B  o  /9_5_n 

Bm  l9  Z  ;  T ,  p 


T 


<  ST6£>  =  0 


(29) 

(30) 


Finally  for  p,  S  and  Z  as  independent  variables. 


2  Rmin  “  HC  6S2+  ^  BS  «P2+  2(It)S,Z  5p6S  + 

P°  (3D 


'V 


Using  (3-12)  to  evaluate  the  determinant  for  S,  p,  we 
have 


<  6p6S>  =  - 
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Hence 


<  6p6S>  =  -pokBTQ6 


(32) 


<  6 p6Z>  =  0 


(33) 


If  we  now  write 

6P  6T  =  6T  +  ^^T^p5^6T  +  ^“^Tap  6 5 


and  make  use  of  (11),  (13)  and  (30),  we  have 
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W  p,  ou  ,  m 

6P  6T>  =  (1^-)  <  6  T2>  =  - T-  B  -° 

3  Ty  p  ~ 


'V 


For 


6P6p  6P^3T^p,Z  6T  +  (3p^T,Z  6p  +  (3Z^T,p  6Zl  5 
and  make  use  of  (11),  (12)  and  (10  to  obtain 


^  r,  o  p  kD  T 

.D  x  v  _  /3  Px  ~  2^  _  Ko  B  o 

6P  6p  3  p  T ,  Z  6p 


V 


o 


For 


6Z  6^3T-P,£  6T  +  ^3P^T,^P  +  ^3  £^T,P  6^5 


we  make  use  of  (25),  (27)  and  (30)  to  obtain 


(3*0 


(35) 


<  6Z  65>  =  kg  Tq 


(36) 


Similarly,  we  have 


<  fiS  6T>  =  (||)p  5<  6T2>  =  kB  Tq 


(37) 


It  is  convenient  to  summarize  the  results  as  in  Table  2 


A . 3  Fluctuations  of  Fourier  components  of  thermodynamic 

Quantities 


For  this  purpose  we  write  the  minimum  work  as 

Rmin  =  h  /  PQ  (6U  “  Vs  +  Po6V  "  Z6?)  dV  »  (38) 

Vo 

where  the  extensive  variables  in  the  integrand  refer  to 
values  per  unit  mass  and  the  integrant  then  represents 
Rmin  per  un ^  v°lume«  By  the  same  procedure  as  before, 
(8)  becomes 


37  Sri  9  i 


>  , 
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Table  2.  Mean  square  fluctuations  in  various  thermodynamic 
quantitieso  (To  obtain,  say,<  6P6p>,  read  the 
intersection  of  corresponding  row  and  column. ) 


Mean 

fluctu¬ 

ations 

6S 

ST 

6P 

6p 

sz 

6S 

kBCP 

kDT 

B  o 

0 

-p  kDT 
Ko  B  o 

0 

~kBCP  X 

x  B 

(il) 
d  Z  '  S ,  P 

6  T 

knT 

B  o 

2 

k-T 

B  o 

OO  00  2 

3  BTkBTQ 

0 

2 

kDT 

B  o  v 

0 

00 

cv 

00 

cv 

00 

cv 

V3  S  T,p 

6  P 

0 

00  OO  P 

6  BTkgTo 

kBToBS 

pokBTo 

OO 

kBToBSy. 

0 

■ 

00 

cv 

V 

0 

Vo 

Po 

C  S ,  P 

6p 

-p  knT 
o  B  o 

0 

pokBTo 

P^k^T 

Ko  B  o 

0 

PokBTo  .. 

x  6 

Vo 

V  Bm 
o  T 

bt 

(3  5) 
ZyT,p 

6Z 

0 

2 

-knT^ 

B  o  v 

“kBToBS^ 

0 

kBTo  X 

ko  T 

B  o 

00 

°v 

^0 

(3_Z) 

3  5JS,p 

(IS) 

V?'T,e 

(— ) 
£;S,P 

<$£ 

-kB°P  X 

0 

0 

PokBToy 

kDT 

B  o 

kBTo  X 

(11) 
v3  Z'S,P 

bt  " 

^9~Z  •'’TjP 

(IS.) 

v3  Z'T,p 

o  e 


A-10 


Rmin  s  35  p0  f  (<SS6T+-i  6P6p  +  6Z6£)  dV 

V  p 

o  o 


(39) 


If  we  take  p,  T  and  Z  as  independent  variables,  (39) 
becomes 


p  1  r  rCv  rm2  1  V  2  ,3^  ,72x  0,9S, 

R  -  -* "  ^  Po  {T  l-T  6T  36p  3  Z  Tjp  Z  +  2  9  Z  T,p 

Vo  p 

o  o 


mm 


6TSZ]  dV 


(40) 


We  now  define 

p(k)  =  — /  6p(r,o)  e”ik,r  d~  r  etc. 


o  V 


o 


Therefore 


(41) 


6p  =  Z  p(k)  e 
k 


lk  c  r 


etc 


(42) 


Therefore 

/6 p2dV  =  Z  /  p ( 5 ) p ( k ’ )  elkc r  elk?°r  dV 

k  ,k 1 


=  V  Z  p ( k) p ( -k ) 


o 


k 


(43) 


Similarly , 


/6T2dV  =  V  Z  T(k)T(-£) 

0  k 

/6T6ZdV  =  V  Z  T ( k )  Z  ( - k )  =  V  Z  T(-6)Z(5) 
0  k  0  k 


(44) 

(45) 


Therefore 


R  .  -  V  Z  r 

min  o  «  { 


k>  0  L  T  _  T  (  k  J  T  (  k  •  +  2  p  '  k  ■  p  ^  ^  ^  +po(3  Z'T,p 


o 


o 


X  Z(k)Z(-k)  +  2  P0(|-|)tjP  T(t)Z(-1c)} 


(46) 


. 


(5  (?(;:  "  Z  7  -  <\ 


' 
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Now  from  (42),  we  have 

P*(k)  =  p(-k)  etc  c  (47) 

Let  us  define 


p(k)  =  P1(k)  +  ipp(k)  etc .  (48) 

Therefore 


P ( -k )  =  p1(k)  -  ip2(k) 

I  p(k)p(-k)  =  Z  (p2  +  p2) 
k  k  1  d 


(49) 

(50) 


From  (45),  we  have 

Z  (T  Z  +  TpZp)  +  i  Z  (T9Z  -  T  Zp) 
k  k  1  1 

» 

=  Z  (T1Z1  +  T2Z2)  -  i  Z  (T2Z1  -  T1Z2) 


Hence  Z  (T9Z1  -  T-.Z,,)  =  0  . 
k  1  1 

Therefore  Z  T(k)Z(-k)  =  Z  (T-.Z,  +  TpZp)  . 

k  k  1  1  ^  ^ 


Hence  (46)  becomes 


R  .  -  V  Z  Z 

min  o  .  n  .  ,  n 

k>0  l ,  =1 , 2 


-poCV  m2 


Br 


T  Ti  + 

0  Po 


z2 

o'a  z'T.p  i 


+  2po(M)T,p  TiZi} 


(51) 

(52) 


(53) 


From  (53)  >  we  can  conclude  that  quantities  with  the 
subscripts  1  and  2  are  statistically  independent.  That 
is 


<  P1P2> 


^T2>  =  <  T^Z2>  =  0  etc. 


(54) 
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Also,  since  the  i=l  case  involve  the  same 
as  the  i=2,  we  have 
2  2 

<  p^>  =<  pp>  ,  <  T-^Z^>  =<  TpZp>  etc. 


coefficient  s 

(55) 


Using  (55)  and  (56),  we  have 

<  p(k)p(-k)>  =  2  <  p^  > 

<  p(k) p(k)>  =  0 

<  T  ( k )  Z  (  -  k )  >  =  2  <  T1Z1> 

<  T( k) Z (k) >  =  0 

Making  use  of  (9)-15)  and  (56),  we  have 

1  fpokBTo1_  pokBTo 
"  poVo  BT  VoBT 

=<  Z ( k ) p ( - k )>  =  0 


00 

P  V  CTT 
Ko  o  V 

k  T 

B  o  (dZs 
pQVo  ^9  C;S,p 

k  T 

KBo  ,9  S  x 

-  ^£;T,p’ 

CV 


<  p  ( k )  p  (  -  k )  > 

<  T(k) p(-k) > 

<  T(k)T(-k) > 

<  Z ( k ) Z ( -k ) > 

<  Z  ( k )  T  ( -  k )  > 


(56) 


(57) 


(58) 


One  can  obtain  all  the  other  fluctuations  quantities  from 
Table  2  by  putting  Vq  to  be  pQ^  there  and  then  divide  by 
PoVQ.  The  results  are  summarized  in  Table  3 . 


Vi  i  *  •  ,  (b;  ‘  its  (£c  )  3nlaU 


0  - 
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Table  3°  Mean  square  fluctuations  of  the  Fourier  transforms 
of  the  various  thermodynamic  quantities .  (To  read, 
see  Table  2 . ) 


Fluctu¬ 
ations  of 
Fourier 
Compo¬ 
nent  s 

s(i?) 

T(S) 

P(k) 

p(k) 

Z(k) 

S(k) 

S(-k) 

kBcp 

k-Q  T 

B  o 

0 

-kBTo6 

0 

"kBcP 

P  v 

O  0 

‘P  v 

O  0 

V 

0 

(3T) 

v3  Z;S,P 

T  ( -k ) 

knT 

B  o 

kDT 

B  o 

00  00  O 

e  Wo 

0 

-k„T2 

B  o  v 

0 

P  v 

O  0 

00 

P  V  cT, 
o  o  V 

00 

P  V  c,r 
o  o  V 

00 

poVv 

(3S) 

^3  E  }  T ,  p 

p(£) 

0 

00  00  O 

3  BTkBTQ 

kBToBS 

P  kDT 
o  B  o 

”kBToBS. 

0 

p oVv 

V 

0 

V 

0 

2 

p  V 

Ko  0 

(9  p) 
v3  £;S,P 

p(-S) 

-w 

0 

P  kDT 
o  B  o 

pokBTo 

0 

k„T 

B  o  „ 

Vo 

V 

o 

VoBT 

bt 

(11) 

^3  Z  y  T  9  p 

Z(-S) 

0 

-k  T2 

B  o 

-kBToBS.. 

0 

knT 

B  o 

knT 

B  o 

oo 

poVocV 

P2V 

0  o 

P  v 
o  o 

P  v 

Ko  o 

(il) 

o  5  '  T  ,  p 

9  £yS9P 

(3  Z) 

^3  C  ;  S  9  p 

5(-f) 

-kBcp  ^ 

0 

0 

k„T 

B  o  „ 

k^T 

B  o 

kDT 

B  o 

p  v2 

o  o 

V  Bm 
o  T 

P  v 
o  0 

P  v 
o  o 

(IT) 

^3  Z  '  S  ,  P 

■ 

(3_1) 

'3  Z'T,p 

(3_i) 

k3  Z  'T ,  P 

^  L. Vi 
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When  a  system  is  close  to  the  critical  point,  the 
quantities  kt  and  Cp,  and  hence  the  fluctuations  in  p 
and  S  become  anomalously  large.  Therefore  the  condition 
for  statistical  independence  of  neighbouring  subvolumes 
is  more  difficult  to  satisfy.  For  small  fluctuations, 
the  homogeneity  of  the  system  is  not  too  much  destroyed 
so  that  one  can  neglect  the  spatial  derivatives  of  the 
inhomogeneity  distribution.  In  the  language  of  Fourier 
transform,  this  is  equivalent  to  saying  that  the  pre¬ 
dominant  components  [jsee  Eq.  (42)]  have  small  k.  If 
the  inhomogeneity  becomes  pronounced  one  cannot  neglect 
the  spatial  derivatives  of  the  fluctuations,  and  must 
include  them  as  corrections.  For  an  isotropic  medium, 

the  first  order  correction  to  the  calculation  of  mean 

2 

square ' fluctuat ion  of  p  must  be  of  the  form  (V6p) 

J]29]  c  The  other  possible  form  of  V"6p  when  integrated 

represents  an  irrelevant  surface  -  effect ,  while  a  form 
2 

like  6p V  6p  can  be  transformed  into  a  term  involving 
2 

(VSp)  plus,  another  irrelevant  surface  term.  Because 

2 

of  the  presence  of  the  terms  (V6p)  ,  the  evaluation  of 
2 

<6 p  >  is  difficult  and  depends  on  the  shape  of  the  region 
also.  The  calculation  of  the  fluctuations  of  the  Fourier 
components  is  easily  done  and  is  of  importance  for  the 
light  scattering  problem.  Since  the  fluctuations  of  p  and 


u  r  :j  br '  .  • : V 

[  -  "  i  o  s  )  .  i  1  ■: . 

i  ~  <-  ■'  *  no  OJU  1  S'  jpa 

4  '  ■> 

t  .  '  1 
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those  of  T  and  Z  are  statistically  independent,  we  shall 
choose  these  as  independent  variables  and  can  in  the 
following  omit  the  T  and  Z  description,  or  assume  that 
they  are  constant.  In  this  case,  we  have  from  (40), 

Rmin  %  /Of  <SP2  +  dV  •  (59) 


where  b  is  a  positive  constant  in  order  that  the  free 

energy  F  is  a  minimum  at  p=p  ,  where  AF  =  R  .  .  Using 

o  mm 

(42)  and  the  procedures  in  obtaining  (57)?  we  have 

p2kRT 

<  p(k)p(-k)>  =  -°  B  °?  ,  (60) 

V(BT+bkd) 

where  V  is  the  scattering  volume. 

2 

In  Brillouin  scattering,  k  can  be  obtained  from 


(1-32). 


A- 4  Fluctuations  of  Fourier  components  for  Multiple 

Ordering  Processes 

If  there  are  present  in  the  liquid  n  ordering 
processes  described  by  ordering  parameter  £^,  . 
corresponding  to  ordering  forces  Z^  .  .  . . . Z  ,  then  Eq. 
(38)  is  generalized  to 

n 

Rmin  =  h  f  po(<SU  "  To6S  +  Po6V  "  Z  Zi  6?i)  dV  0  (6l) 

i=l 

If  we  take  p,  T,  Z-^  ....  Z  as  independent  variables, 
(6l)  becomes 


, 

. 


O?  .  I  i :  e  ■■ 
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+ 


2  E  ( 
i 


8  S  ^ 

9  Z. 


l 


T,  p 


6T6Z.1  dV 

iJ 


o 


(62) 


If  we  define  p(k),  T(k)  and  Z^(k)  etc.  as  in  (4l), 
and  follow  the  same  procedure  as  before,  we  get 


<  p ( k) p ( -k) >  = 


<  T(k)T(-k)>  = 


p2kuT 
o  B  o 

V  Bt 

kuT^ 

B  o 

oo 

P  Vc T7. 
o  V 


<  z.(5)Zl(-k)>  = 

Ko  1  5 


<Z.(k)Z.(-k)>  =  0  for  i  ^  j 
^  J 

<  Zi(k)p(-S)>  =  0 


<  T(k)p(-k)  >  =0  etc. 


(63) 


The  other  fluctuation  quantities  can  be  obtained  easily » 


. 
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List  of  Important  Symbols 


c 

c 

c 


P 

V 


2ib 

I 


n(r ,t ) 
T 


bulk  viscosity 

adiabatic  bulk  modulus 

isothermal  bulk  modulus 

velocity  of  light  in  vacuum 

isobaric  specific  heat 

isochoric  specific  heat 

shear  modulus  at  infinite  frequency 

integrated  intensity  of  the  Brillouin  components 

integrated  intensity  of  the  Rayleigh  component 

\  *  ■ 

wave  number  of  hypersonic  sound  wave 
Boltzmann  constant 

number  density  of  the  scattering  medium 
temperature 

velocity  of  temporally  absorbed  sound  wave 
velocity  of  spatially  absorbed  sound  wave 
sound  velocity  at  zero  frequency 


v 

v 


B 


sound  velocity  at  infinite  frequency 

sound  velocity  measured  from  the  Doppler  shifted 

Brillouin  components 


< 

a 

o 

> 

II 

k 

TSV 

X1 

:=  V 

O 

k 

tTV 

x0 

-  V 

k 

T 

2 

0 

n 

X 

-  V 

k 

T 

n 

0 

n 

w"  X 

. 


order  force 


This  symbol  has  been  used  to  denote  different 
quantities,  namely  the  polarizability 
[Chap,  lj ,  the  sound  absorption  per  wave¬ 
length  [Chap. 2]  and  a  dimensionless  quantity 
involving  the  thermal  conductivity  [Chaps.  5 3  7] 
average  polarizability  of  the  molecules  in  the 
scattering  medium 
coefficient  of  volume  expansion 
sound  absorption  per  second 
dielectric  constant 
shear  viscosity 
specific  heat  ratio 
wave  number  of  the  incident  light 
wave  number  of  the  scattered  light 
adiabatic  compressibility 
isothermal  compressibility 

thermal  conductivity 
wavelength  of  the  incident  light 

angular  frequency  of  the  hypersonic  sound  waves 


—  X 


r 


p(r,t ) 


T 

n 

0 


angular  frequency  of  the  incident  light 
angular  frequency  of  the  scattered  light 
density  of  the  scattering  medium 
shear  relaxation  time 

T/v  k.  It  also  denotes  the  scattering 
o 

angle . 


